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PREFACE 
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This  contract  involves  the  development  of  techniques  to  process 
multi -temporal  remote  sensing  data  for  purposes  of  extraction  of  hydro- 
graphic  information.  The  techniques  and  processing  software  developed 
under  this  contract  were  based  on  multi -date  analysis  of  a  set  of  previ¬ 
ously  processed  Landsat  scenes  covering  the  Bahamas  study  region. 
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SUMMARY 


In  order  to  enhance  its  Digital  Image  Processing  System  (DIPS)  cap¬ 
ability  the  Defense  Mapping  Agency  (DMA)  has  supported  the  development 
of  a  multi -temporal  procedure  (MTP).  This  procedure  has  been  designed 
with  maximum  flexibility  to  allow  the  operator  to  apply  the  software  to 
a  variety  of  hydrographic  applications  and  remote  sensing  data  sources. 
A  goal  at  DMA  is  to  develop  processing  technology  for  passive  remote 
sensors  which  minimizes  the  need  for  ship  supported  surface  truth  meas¬ 
urements.  With  the  present  capability  the  DIPS  cannot  be  used  to  remove 
unwanted  noise  and  effects  which  can  influence  the  depths  predicted  from 
satellite  sensor  data.  Multi -temporal  processing  provides  a  means  to 
diminish  noise  and  separate  the  effects  due  to  temporal  phenomenon  such 
as  turbidity,  haze,  and  surface  slicks.  The  multi -temporal  software 
developed  for  the  DIPS  allows  the  operator  to  perform  basically  nine 
separate  operations.  These  include  display  functions  for  loading, 
viewing,  and  combining  multi-date  subscenes  which  have  been  previously 
co-registered.  Options  are  also  provided  for  image  smoothing  and 
polygon  subarea  selection.  The  polygons  can  be  examined  to  determine 
depth  statistics  and  depth  differences  for  selected  dates.  Further  a 
SCATTERPLOT  and  REGRESSION  option  allows  the  operator  to  investigate  the 
relationships  between  predicted  depths  for  several  dates  and  adjust,  if 
necessary,  the  water  depth  equation  parameters.  With  the  parameters 
adjusted  between  scene  date*  the  operator  can  recalculate  the  water 
depths  for  each  date  and  weight  average  the  results  to  eliminate 
unwanted  noise.  The  resulting  predicted  depths  can  then  be  used  as  a 
basis  for  additional  parameter  adjustments  in  an  effort  to  further 
resolve  date-to-date  differences  in  predicted  depth.  Once  the  operator 
is  satisfied  with  the  result,  a  relationship  can  be  established  between 
the  original  depth  prediction  for  each  scene  date  and  the  final  best 
predicted  depth  for  the  calibration  polygon(s).  The  APPLY  function  can 


then  be  used  to  modify  each  of  the  scene  dates  and  obtain  a  best  esti¬ 
mated  depth  for  the  entire  subarea.  The  flexibility  in  the  options  of 
the  MTP  software  allows  the  definition  of  many  separate  procedures  to 
extract  a  best  estimated  water  depth  from  temporal  data. 

A  multi -temporal  data  set  was  constructed  from  six  previously  pro¬ 
cessed  Landsat  scenes  covering  portions  of  the  Bahamas  Photobathymetric 
Calibration  Area.  The  six  scenes  were  brought  into  registration  using 
common  ground  control  positions  and  a  series  of  affine  transformations. 
Once  these  images  had  been  satisfactorily  registered,  a  semi-rigid 
Landsat  imaging  model  was  used  to  locate  pixels  corresponding  to  SAI 
ship  transect  depths  in  the  scene  [1].  Errors  in  the  registration 
process  were  found  to  be  on  the  order  of  two  pixels  in  each  direction, 
while  the  errors  associated  in  the  location  of  ship  data  were  within  two 
pixels.  Four  data  sets  were  assembled  for  calibration  areas  3A,  3B,  3C, 
and  3D  (see  Figure  1,  page  10).  Each  of  these  data  sets  consisted  of 
the  average  ship-measured  depth  over  each  Landsat  pixel  along  with  the 
six  individual  predicted  depths  as  derived  from  the  Landsat  signal 
levels.  These  data  were  analyzed  to  gain  insight  into  the  characteris¬ 
tics  of  multi -temporal  data  and  as  test  cases  for  calibration  and 
evaluation  of  suggested  procedures. 

Large  offsets  (0.5-5.0m)  in  mean  water  depth  were  observed  between 
the  Landsat  predicted  depths  and  those  provided  from  ship  measurements. 
Of  data  from  the  six  available  dates,  data  from  three  were  found  to  be 
sufficiently  noisy  to  caution  their  use  in  any  multi -temporal  analysis. 
The  multi-date  algorithm  showed  improvement  over  the  best  single  date 
results  for  the  case  where  ship  survey  data  were  utilized  and  also  for 
the  case  where  a  best  depth  estimate  was  formulated  based  only  on  the 
Landsat  signal  values.  Observed  improvements  were  found  to  be  comparable 
to  that  predicted  from  rms  noise  reduction.  It  was  concluded  that  the 
algorithm  effectiveness  may  be  improved  by  implementing  pre-water  depth 
processing  procedures  design  to  remove  systematic  noise  components. 
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INTRODUCTION 

The  Defense  Mapping  Agency  has  responsibility  for  issuing  bathy¬ 
metric  charts  for  all  areas  of  the  world  outside  the  United  States. 
Costs  of  revising  charts  using  conventional  ship  survey  methods  have 
increased  sharply  in  recent  years.  As  a  result  DMA  is  no  longer  able  to 
meet  the  demands  for  accurate  charts  and  has  been  increasing  its  effort 
to  develop  a  Satel 1 ite/Ai rborne  remote  sensing  survey  capability  which 
would  allow  charts  to  be  updated  in  a  more  efficient  manner.  ERIM  has 
participated  in  numerous  studies  supported  by  both  DMA  and  NASA  to  de¬ 
velop  and  refine  techniques  that  predict  water  depth  based  upon  Landsat 
radiometric  parameters  [2,3].  In  each  of  these  previous  studies  water 
depth  algorithm  development  was  accomplished  by  relating  measured  depths 
to  Landsat  radiometric  parameters.  Altogether  these  studies  provide  a 
significant  data  base  for  the  evaluation  of  remote  sensing  techniques. 
Since  many  of  these  studies  were  conducted  in  the  Bahamas,  DMA  desig¬ 
nated  this  region  as  the  Bahamas  Photobathymetric  Calibration  Area.  The 
waters  in  this  region  are  exceptionally  clear  and  exhibit  a  wide  variety 
of  flora  and  bottom  types  [4]. 

Based  on  the  algorithms  developed  in  the  above  studies,  DMA  has 
supported  the  development  of  software  to  be  run  on  its  Digital  Image 
Processing  System  (DIPS).  With  the  present  capability,  however,  DMA  has 
no  way  to  remove  unwanted  noise  and  effects  which  can  easily  influence 
predicted  depth  calculations  derived  from  Landsat  data.  Processing  of 
bathymetric/hydrographic  data  images  can  require  detailed  analysis  in 
each  of  two  or  three  spectral  bands.  When  it  is  necessary  to  separate 
time-varying  phenomenon  as  turbidity,  surface  slicks,  clouds  and  haze 
from  bathymetric  features,  multi-date  imagery  is  required.  Given  that 
multi-date  imagery  is  frequently  vailable  for  purposes  of  image  selec¬ 
tion  and  to  identify  tt  ,  ..VI  atures  it  is  reasonable  to  consider  de¬ 
veloping  algorithms  which  can  exploit  co-registered  multi-date  imagery. 
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Work  under  this  contract  consists  of  development  of  a  multi-temporal 
processing  procedure  for  application  on  the  DIPS.  The  present  DIPS 
system  is  still  in  the  development  stages  and,  therefore,  processing 
algorithms  developed  for  the  DIPS  must  be  flexible  and  sufficiently 
generalized  so  as  to  alter  application  for  use  with  a  variety  of  data 
sources  and  applications.  Once  procedures  have  been  proven  through 
repeated  use,  then  processing  software  on  the  DIPS  can  be  modified  to 
reduce  present  time  consuming  constraints.  In  this  effort,  work  was 
directed  toward  utilizing  multi -date  Landsat  coverage  to  develop  a 
procedure  which  will  minimize  the  influence  of  noise  and  other  effects 
such  as  varying  bottom  reflectance  and  water  clarity  ’n  order  to  provide 
a  best  estimate  of  water  depth. 

2.1  BACKGROUND 

The  Defense  Mapping  Agency  at  its  Hydrographic/Topographic  Center 
(HTC)  has  the  capability  to  process  Landsat  MSS  data  to  produce  water 
depth  maps.  The  algorithms  used  are  a  single  channel  algorithm  based  on 
digital  values  in  band  MSS4  (green)  and  a  two  channel  algorithm  based  on 
the  ratio  of  digital  values  in  bands  MSS4  and  MSSS  (red)  [5].  The 
Landsat  estimates  of  water  depth  contain  errors  caused  by  changes  in 
water  clarity,  tidal  state,  bottom  reflectance,  surface  reflected 
energy,  atmospheric  effects,  and  sensor  noise.  The  algorithms  require 
estimates  of  water  clarity  (K,  the  irradiance  attenuation  coefficient) 
and  bottom  reflectance.  These  parameters  are  entered  as  constants  in 
the  program.  If  they  are  estimated  incorrectly,  depth  errors  will 
result  in  the  processed  data. 

Random  noise  effects  can  be  reduced  if  more  than  one  Landsat  data 
set  can  be  analyzed.  Before  the  beneficial  effects  of  noise  reduction 
(through  averaging  the  depth  estimates  made  on  two  or  more  scenes)  can 
be  realized,  the  systematic  errors  between  scenes  must  be  reduced  to  low 
levels.  This  reduction  can  be  accomplished  by  adjusting  the  parameters 
used  to  process  the  various  MSS  data  :  °ts  to  minimize,  in  a  least 
squares  sense,  the  differences  in  water  depths  computed  from  the  scenes 
being  analyzed. 


The  parameter  adjustment  process  requires  precisely  registered 
scenes  of  Landsat  data,  so  that  pixel -by-pixel  comparisons  can  be  made. 
The  technology  exists  to  accurately  correct  scenes  if  a  few  (5-10  per 
scene)  ground  control  points  are  available  and  the  spacecraft  attitude 
is  known.  The  latter  information  has  been  available  in  the  X-format 
tapes  provided  by  EROS  Data  Center. 

Two  cases  of  parameter  adjustment  can  be  distinguished;  a  case 
where  a  few  known  depth  points  are  available,  and  a  case  where  no 
ancillary  depth  information  is  available.  In  the  first  case,  water 
depth  estimates  from  a  reference  scene  are  first  corrected  to  the  known 
data  by  adjustment  of  the  algorithm  parameters  to  minimize  the  differ¬ 
ence  between  the  estimates  of  depth  (from  the  Landsat  data)  and  the 
actual  depths.  After  the  reference  scene  has  been  adjusted,  each  of  the 
additional  scenes  can  be  brought  into  correspondence  with  the  reference 
scene  by  a  similar  parameter  adjustment  procedure.  At  each  step  of  the 
parameter  adjustment  procedure,  the  resulting  revised  parameters  should 
be  checked  to  assure  that  the  adjustments  are  reasonable.  If  unreason¬ 
able  adjustments  arise  from  the  least  squares  procedure  it  is  an  indica¬ 
tion  that  something  may  be  wrong  with  the  data  set  being  analyzed.  In 
the  case  where  no  ancillary  bathymetric  data  are  available,  the  esti¬ 
mates  of  water  depth  from  the  various  Landsat  scenes  can  be  brought  into 
correspondence  by  adjusting  algorithm  parameters  to  minimize  the  mean 
square  depth  difference  between  the  scenes.  But  because  of  uncertain 
tidal  state  and  solar  irradiance  and  bottom  reflectance  effects,  the 
average  computed  depths  may  be  biased  with  respect  to  true  depths. 

2.2  OBJECTIVES  OF  THE  PRESENT  STUDY 

The  present  study  had  four  objectives.  First,  the  mathematical 
details  of  the  parameter  adjustment  procedure  were  developed.  Second,  a 
multi -temporal  data  set  was  assembled  from  six  previously  processed 
Landsat  scenes  covering  the  Bahamas  Photobathymetric  Calibration  Area. 
Third,  the  procedure  was  applied  to  the  composite  Bahamas  data  set  and 
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the  results  analyzed.  Fourth,  software  was  written  to  perform  the 
required  analyses  within  the  context  of  the  DIPS  at  DMA/HTC.  In  the 
remainder  of  this  report  the  results  of  the  four  efforts  are  discussed. 

2.3  EXISTING  DIPS  CAPABILITY 

The  DIPS  provides  a  basic  operating  capability  to  process  Landsat 
MSS  data  and  other  sources  of  remote  sensing  imagery  to  detect  and 
position  unknown  navigational  hazards  or  update  charts  which  poorly 
describe  hazard  features.  The  DIPS  provides  real  time  interactive 
display  and  manipulation  capabilities  that  allow  the  operator  to  process 
one  or  more  Landsat  bands  for  purposes  of  extracting  hydrographic 
information  in  the  form  of  predicted  water  depths  or  location  of 
specific  bottom  features.  When  fully  operational  DMA/HTC  plans  to  use 
the  DIPS  to  support  the  following  hydrographic  work: 

(1)  Evaluation  of  hydrographic  charts  for  accuracy. 

(2)  Updating  and  chart  revision. 

(3)  Provide  regular  inputs  to  Notice  to  Mariners  reports. 

(4)  Confirm  and  position  doubtful  dangers. 

(5)  Provide  planning  inputs  to  shallow  water  hydrographic  ship 

surveys. 

(6)  Provide  a  monitoring  function  for  unstable  navigational 

hazards. 

Presently  the  DIPS  image  processing  and  analysis  functions  are 
limited  to  the  spatial  units  of  a  single  display  image  (512  x  512). 
Each  display  image  (subarea)  can  be  transformed  into  geographic 
coordinates  with  the  aid  of  operator  selected  ground  control  points. 
The  image  warp  function  can  be  used  to  transform  geographic  coordinates 
of  ship  soundings  to  image  line  and  point  coordinates.  Signal  levels  of 
these  latter  pixels  can  be  used  to  perform  a  linear  regression  analysis 
yielding  a  relationship  between  water  depth  and  signal  level. 
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The  linear  equation  is  a  logarithmic  transformation  of: 


w  _  u  +  y 
v  “  s  o 


where 


V  =  Landsat  signal  count 
z  =  water  depth  (m) 

K  =  irradiance  attenuation  coefficient 

V$  =  Average  deep  water  signal 

VQ  =  Average  V-V$  signal  for  zero  depth 


where:  E0 

p 

G 

T 


solar  irradiance  at  ocean  surface 
bottom  reflectance 
scanner  sensitivity  constant 
atmospheric  transmittance 


(1) 


(2) 


The  logarithmic  equation  has  the  linear  form: 

Y  =  A  +  BZ  (3) 


where  Y  =  ln(V  -  V  ) 

A  =  ln(VQ) 

B  =  -2K 


Equation  3  above  can  also  be  used  to  express  water  depth  in  terms  of  a 
two  band  ratio  (MSS4  and  MSS5). 

Presently  there  is  no  capability  on  the  DIPS  to  extend  depth  pre¬ 
dictions  derived  on  one  subarea  to  an  adjacent  one  or  to  mosaic  proces¬ 
sed  subareas. 


3.0 


MULTI -TEMPORAL  PROCEDURE  DEVELOPMENT 

The  multi -temporal  processing  procedure  developed  for  the  DIPS  is 
one  which  utilizes  multi-date  Landsat  data  to  remove  uncertainty  in  the 
depth  calculation  input  parameters  which  are  either  assumed  or  measured 
at  selected  point  locations  in  the  scene.  In  this  regard  the  calcula¬ 
tion  of  water  depths  from  Landsat  data  depends  on  knowledge  of  three 
basic  parameters  for  each  scene  date.  These  are  the  deep  water  signal 
V$,  the  irradiance  attenuation  coefficient  K,  and  the  bottom  reflecti¬ 
vity  p.  In  the  water  depth  algorithm  presently  on  the  DIPS  it  is 
assumed  that  the  deep  water  signal  is  constant  throughout  the  scene. 
However,  spatially  varying  surface  and  atmospheric  effects  can  lead  to 
significant  errors  in  this  term.  The  extinction  coefficient  could  be 
highly  variable  both  spatially  and  temporally.  The  bottom  reflectivity 
could  also  display  large  spatial  variance,  but  temporal  changes  can  be 
expected  to  be  isolated  if  they  exist  at  all.  In  addition,  Landsat 
image  characteristics  including  striping  and  angular  distortions  will 
affect  predicted  depths.  Thus  there  exist  substantial  spatial  and 
temporal  complexities  in  the  determination  of  water  depths.  Under  these 
circumstances  of  parameter  uncertainty,  the  multi -temporal  technique 
becomes  an  attractive  approach.  The  multitemporal  procedures  as  des¬ 
cribed  below  were  developed  using  a  six  date  scene  set  covering  a 
portion  of  the  Bahamas  Calibration  Area.  More  specifically,  Landsat 
derived  water  depths  from  study  areas  3A,  3B,  3C,  and  3D  were  combined 
with  1980  ship  transect  depths  to  form  a  test  set  (See  Figure  1).  The 
specifics  of  the  test  set  development  are  discussed  in  section  4.0. 

3.1  BACKGROUND  ANALYSIS 

In  the  simplest  form  of  the  problem,  given  perfect  knowledge  of  V$, 
K,  and  p  and  with  noise  effects  spatially  uniform  and  comparable  for 
each  of  the  individual  scenes,  we  would  be  able  to  determine  the  average 
depth  or  "best"  estimate  at  each  pixel.  However,  under  more  realistic 
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conditions  we  have  imperfect  parameter  information  and  possibly  spatial¬ 
ly  varying  noise.  Under  such  conditions  we  desire  to  utilize  the  multi¬ 
temporal  data  in  a  methodology  which  will  produce  a  "best"  estimate  of 
water  depth.  There  are  essentially  two  types  of  uncertainty  which  in 
turn  suggest  very  different  approaches  for  obtaining  a  best  estimate. 
First,  with  V$  there  is  uncertainty  in  the  measurement  accuracy  and 
applicability  to  the  depth  determi nations  at  other  points  in  the  scene. 
If  the  variation  in  the  Vg  term  is  random,  then  the  errors  may  average 
out  to  some  extent  when  the  multi -temporal  results  are  combined.  If  on 
the  other  hand  the  variations  are  due  to  patterns  in  atmospheric  haze, 
then  it  seems  essential  that  such  haze  be  first  normalized  throughout 
the  scene.  It  is  recommended  that  such  a  haze  algorithm  should  be  added 
to  the  DIPS  processing  software. 

A  measured  difference  in  Vs  from  scene  date  to  scene  date  cannot  be 
used  to  improve  the  determination  of  K  or  p.  The  bottom  reflectance  is 
not  expected  to  change  temporally.  But  because  it  is  difficult  to 
estimate  algorithms  which  minimize  the  effects  of  changing  bottom 
reflectance  on  water  depth  calculations  are  being  developed.  The  zero 
depth  signal  VQ  contains  the  bottom  reflectance  coefficient  and  can  be 
either  estimated  from  the  data  or  calculated  using  solar  irradiance, 
sensor  responsivity ,  atmospheric  transmission  and  bottom  reflectance  as 
in  equation  (2).  The  values  of  V$  and  VQ  must  be  tied  to  a  subscene 
area.  In  fact  the  VQ  term  varies  from  pixel -to-pixel  in  the  scene  but 
cannot  be  directly  calculated  at  each  point  without  knowledge  of  p. 
Because  VQ  is  the  product  of  a  series  of  parameters,  variation  in  VQ 
from  scene  date  to  scene  date  is  not  directly  related  to  variations  in 
bottom  reflectance.  Knowledge  of  the  water  attenuation  coefficient,  K, 
could  on  the  other  hand,  be  useful  if  certain  assumptions  are  permissi¬ 
ble.  First,  if  the  K  value  does  not  change  from  date  to  date  and  there 
is  a  significant  and  known  change  in  the  tidal  state,  then  it  is 
possible  to  calculate  the  value  of  K  at  each  pixel  given  a  value  of  V$ 
and  VQ.  Second,  if  the  K  value  has  changed  temporally  by  a  known 
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quantity  then  it  is  possible  to  estimate  the  water  depth  independently 
of  the  bottom  reflectance. 

It  is  questionable  whether  either  of  these  approaches  is  applicable 
to  the  Bahamas  region  since  the  tidal  changes  are  small  and  on  the  order 
of  the  water  depth  errors.  In  addition  there  is  no  data  to  support  any 
uniformity  in  the  reported  tidal  state.  Variations  in  tidal  state  can 
be  expected  as  a  function  of  bottom  slope  and  depth  patterns.  There  is 
little  reported  data  on  spatial  variations  of  K  values.  It  is  likely 
that  spatial  variations  in  K  exceed  those  due  to  temporal  changes  with 
the  possible  exception  of  those  caused  by  passage  of  large  storms. 

3.2  MULT I -TEMPORAL  PROCESSING  METHODOLOGY 

With  this  background  let  us  now  explore  possible  methodologies  for 
obtaining  a  best  depth  estimator.  First  consider  a  case  where  one  has 
only  very  limited  water  depth  soundings  as  may  be  available  from  a  crude 
chart.  It  is  further  assumed  that  these  depths  are  suitable  for  pur¬ 
poses  of  checking  or  validating  the  results  obtained  by  processing 
remote  sensing  data  but  insufficient  by  numerous  parameter  estimation. 
Initial  values  of  water  depths  are  obtained  by  making  reasonable  assump¬ 
tions  for  K  and  deriving  Vs  and  VQ  from  each  of  two  to  four  individual 
data  sets.  Under  such  circumstances  an  approach  is  sought  which  will 
utilize  the  multi-date  information  to  obtain  results  which  are  superior 
to  those  from  a  single  date.  If  one  attempts  to  apply  an  iteration 
and/or  relaxation  process  over  the  parameters  of  interest  it  is  soon 
discovered  that  there  is  no  criterion  available  for  testing  convergence. 
A  plot  of  depths  calculated  on  one  date  versus  those  calculated  on  a 
second,  for  a  set  of  registered  pixels  and  a  range  of  depths,  can 
suggest  two  types  of  parameter  changes  (Figure  2).  A  regression  slope 
not  equal  to  one  suggests  a  change  in  K  while  an  offset  indicates  a 
change  in  VQ.  If  this  latter  tenm  is  calculated  rather  than  estimated 
from  the  data,  then  the  offset  may  be  due  to  changes  in  irradiance  or 
tidal  state.  Once  the  slope  deviation  has  been  removed  by  adjustment  of 
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one  date's  K  value  and  the  tidal  change  is  used  to  remove  the  offset, 
the  two  data  sets  are  in  agreement  outside  of  the  possibility  of  syste¬ 
matic  scattering  (spatial  variations)  about  the  regression  line  or  as 
residual  differences  related  to  spatial  location  or  water  depth.  If  the 
bottom  reflectance  parameter  is  perturbed  in  an  effort  to  obtain  better 
agreement,  one  realizes  that  all  such  perturbations  merely  move  the  pos¬ 
ition  of  the  point  along  a  line  parallel  to  the  slope  (note  that  such 
perturbations  are  taken  to  be  the  same  for  each  date).  In  this  case  the 
information  necessary  to  reconcile  the  bottom  reflectance  on  a  pixel -by¬ 
pixel  basis  is  not  present.  If  the  slope  in  the  original  plot  were 
greatly  different  than  one,  implying  a  large  change  in  K,  then  one 
could,  in  principle  increase  or  decrease  the  bottom  reflectance  value  of 
the  individual  pixels  to  bring  them  into  closer  agreement.  However, 
unless  the  atmospheric  spatial  variations  are  first  removed  from  the 
data  such  results  are  meaningless  and  such  variations  due  to  path  radi¬ 
ance  must  be  removed  from  the  data  before  beginning  the  depth 
processing. 

Presently  there  is  no  atmospheric  correction  capability  with  the 
DIPS  software.  Further  our  experience  with  the  Bahamas  data  set  demon¬ 
strates  only  very  slight  changes  in  K  value  from  date  to  date  and  insuf¬ 
ficient  ones  from  which  to  make  any  attempt  to  analyze  possible  spatial 
variations  in  bottom  reflectance.  Any  number  of  scene  dates  can  be  re¬ 
conciled  into  agreement  by  adjustment  of  K  and  the  offset,  and  the  resi¬ 
duals  can  be  used  to  indicate  any  systematic  differences.  Once  these 
are  renoved,  the  date-to-date  residuals  will  have  a  random  character  and 
no  further  parameter  adjustment  to  improve  agreement  is  possible.  In 
this  circumstance  the  average  depth  computed  for  each  pixel  becomes  the 
"best  estimator".  If  the  spatial  noise  properties  of  data  sets  are 
greatly  different,  then  weights  related  to  the  noise  amount  can  be  de¬ 
rived  for  each  data  set,  and  the  average  computed  as  a  weighted  average. 
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Under  a  separate  set  of  circumstances  the  multi -date  data  would  be 
accompanied  with  a  large  number  of  ship  survey  soundings  which  are 
suitable  to  both  calibrate  and  check  the  predicted  depths.  If  these 
depths  are  evenly  distributed  over  the  scene,  then  parameter  adjustments 
derived  from  these  water  depth  measurements  would  in  principle  be  valid 
over  the  entire  scene.  If  on  the  other  hand  these  depths  are  derived 
from  transect  data,  then,  of  course,  the  validity  of  extrapolating  water 
depth  remote  sensing  parameter  results  to  other  portions  of  the  scene  is 
in  doubt.  Given  the  ship  depths,  our  objective  is  to  not  only  adjust 
the  predicted  depths  from  date-to-date  but  to  make  further  parameter 
adjustments  to  minimize  the  differences  between  measured  and  predicted 
depths  with  a  least  squares  criteria.  In  this  approach  K  and  p  para¬ 
meters  can  be  estimated  initially  and  fed  into  a  system  of  equations, 
such  as  shown  below,  which  adjust  each  parameter  by  some  small  change  so 
as  to  minimize  the  difference  in  the  mean  squared  error  over  a  given 
series  of  N  available  scenes  and  M  pixel  locations. 

3Z  9Z 

Z-s  i  =  Z.--  +  (  )  AK.  +  (  )  Apj  j  -  1 . N  (4) 

J  J  9K  ij  J  3Pij  J  i  =  1 . M 

/N 

where  z. .  are  measured  water  depths, 

*  J 

z,.  are  Landsat  estimated  depths. 

•  J 

The  resulting  adjusted  parameters  are  then  fed  back  into  the  same 
set  of  equations  in  an  iterative  process  so  as  to  converge  to  a  best 
estimated  water  depth.  In  formulating  the  equations  one  must  be  careful 
to  make  the  number  of  independent  equations  substantially  greater  than 
the  number  of  parameter  unknowns  so  that  the  system  can  produce  a  stable 
solution.  For  example  suppose  there  are  M  ship  depth  locations  and  N 
scene  dates.  Then  one  can  write  up  to  N  x  M  equations  in  N  x  M  unknowns 
but  would  need  to  have  substantially  fewer  unknowns  to  obtain  a  stable 
solution.  If  one  allows  K  and  tidal  state  (T)  to  vary  temporally  but 
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not  spatially  then  K  and  T  produce  2M  unknowns.  If  the  bottom  reflec¬ 
tance  is  assumed  to  vary  spatially  but  not  temporally  then  we  have  an 
additional  N  unknowns.  The  variations  in  solar  downward  irradiance, 
which  can  also  be  thought  of  as  a  variation  in  the  Vs  term,  are  perhaps 
equally  variable  in  both  space  and  time.  However,  this  latter  parameter 
variation  leads  to  N  x  M  unknowns  making  it  unsuitable  to  this  iteration 
procedure.  Thus,  it  is  again  seen  that  a  procedure  is  needed  to  first 
remove  the  spatial  atmospheric  variations  prior  to  water  depth  proces¬ 
sing.  The  extendabi 1 ity  of  the  bottom  reflectance  parameters  to  other 
portions  of  the  scene  is  certainly  a  dubious  procedure.  However,  one 
could,  in  the  event  that  there  exist  large  K  changes  from  date  to  date, 
use  the  bottom  reflectances  as  determined  from  the  iterative  solution 
using  measured  depths  to  define  and  validate  a  bottom  reflectance  ad¬ 
justment  procedure  as  discussed  above  which  could  be  applied  to  all 
water  pixels  of  the  scene. 

Without  the  atmospheric  and  spatial  noise  normalization  procedure, 
ship  survey  soundings  can  be  used  to  adjust  the  K  irradiance  attenuation 
parameter  and  remove  depth  offsets  due  to  changes  in  tidal  state  or 
other  parameters  affecting  the  VQ  term.  In  this  case  the  parameter 
adjustments  can  be  expressed  as  shown  in  equation  (4)  above,  then  use  an 
iterative  process  of  first  determining  a  set  of  parameter  adjustments, 
then  substituting  the  adjusted  parameters  back  into  the  same  equation  to 
converge  to  a  least  squares  parameter  fit.  When  the  terms  in  the  itera¬ 
tion  equation  are  linear  in  water  depth,  the  described  iteration  process 
is  equivalent  to  linear  regression  analysis  and  the  parameters  defined 
by  regression  must  necessarily  be  the  same  least  squares  solution  as 
that  obtained  with  the  iteration  process.  Thus  linear  regression  can  be 
used  to  adjust  assumed  K  values  in  order  to  bring  them  into  line  with 
the  measured  depths.  Once  K  and  p  parameter  adjustments  have  been 
completed  for  two  or  more  scenes,  residual  differences  can  be  examined 
for  any  systematic  patterns  with  depth  and  those  can  possibly  be  removed 
from  the  data  with  some  further  parameter  adjustment.  At  this  point  the 
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predicted  water  depths  for  different  dates  are  essentially  equivalent 
except  for  random  differences.  The  best  depth  estimator  becomes,  as  in 
the  previous  case,  the  weighted  average  of  predicted  depths  from  each 
scene  date. 

Based  upon  the  analysis  of  the  available  Landsat  multi -date  water 
depth  .naps  in  the  Bahamas  calibration  study  area,  a  general  processing 
procedure  is  described  below.  The  analyses  of  these  data  are  described 
in  section  5.0. 

3.3  DESCRIPTION  OF  PROCESSING  PROCEDURES 

Preparation  procedure:  Using  existing  DIPS  software  and  procedures 


(1)  View  Landsat  band  4  and  select  desired  subscenes  for  multi - 

temporal  processing  (i.e.,  512  x  512  blocks). 

(2)  Locate  geometric  control  points. 

(3)  Register  by  warping  subscenes  from  different  Landsat  dates. 

(4)  Outline  deep  water  and  zero  water  depth  areas  and  calculate  the 

deep  water  signal,  V$,  and  zero  depth  signal,  VQ-VS, 
respectively. 

(5)  Use  available  supporting  data  to  make  best  guesses  of  K,  p,  and 

tidal  state  for  each  subscene  date.  At  this  point  assume 
that  K  and  p  are  constant  for  each  date.  Use  these  para¬ 
meters  to  estimate  water  depths  for  each  pixel  and  date 
within  the  subscene. 

3.4  MULT  I -TEMPORAL  PROCESSING  ELEMENTS 

Use  new  MTP  DIPS  software  to  resolve  date-to-date  depth  differences 
and  obtain  a  "best"  estimator.  Operator  selects  appropriate  command 
process  from  the  following  menu.  There  are  four  basic  commands  within 
the  menu:  LOAD,  POLYGON,  SCATTERPLOT,  and  APPLY.  Each  conmand  will 
have  several  operator  selected  menu  options  as  described  below. 


(1)  LOAD  -  Screen  subscene  for  selected  dates.  Operator  selects 
dates  and  the  program  loads  one  date  depth  file  into  each  of 
the  3  Comtal  image  planes.  Operator  has  option  to  load  two 
depth  files  and  their  difference  map  which  is  calculated  with 
this  routine  with  an  offset  value  of  128.  Once  image  planes 
are  loaded,  operator  can  proceed  with  selection  of  polygon  test 
areas. 

(2)  POLYGON  -  Define  study  areas  for  multi-date  analysis  using  a 
cursor  driven  polygon  selection  routine.  These  areas  may  have 
uniform  bottom  or  K  value.  Areas  selected  should  contain  loca¬ 
tions  of  available  measured  or  otherwise  known  depths.  Operator 
can  select  multiple  polygons  within  the  subscene  as  a  single 
study  area  set  using  this  command.  The  program  stores  data  for 
all  of  the  study  sets  selected  for  each  of  the  available  dates, 
including  estimated  water  depth  and  the  Landsat  counts  in  bands 
4  and  5.  All  of  these  data  are  stored  in  a  single  file  with 
appropriate  name  and  type  designators.  Operator  may  optionally 
exclude  one  or  more  dates  or  selected  portions  of  the  study 
area  from  further  processing.  Operator  can  use  this  command  to 
combine  one  or  more  study  sets  from  the  same  subscene.  Opera¬ 
tor  may  terminate  the  command  by  requesting  the  statistics  of 
the  study  set  (pixel  count,  mean  values,  range,  and  standard 
deviation  about  the  mean).  Operator  can  also  use  this  command 
to  estimate  the  uncertainty  in  the  depth  predictor  for  an 
individual  date  by  selecting  a  study  area  with  uniform  water 
depths. 

(3)  SCATTERPLOT  -  Adjust  K  and  p  parameters  between  dates  to  mini¬ 
mize  date-to-date  differences  in  a  least  squares  sense  and  sub¬ 
sequently  obtain  the  "best"  estimate  of  water  depth  for  each 
pixel  of  the  study  area.  This  command  has  several  operator 
controlled  options. 
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(a)  Input  measured  or  estimated  reference  control  water  depths 
as  polygon  areas,  transects,  or  points.  These  depths  have 
1 atitude/ longitude  coordinates  which  must  first  be  con¬ 
verted  into  line  and  point  coordinates.  Merge  measured 
depth  information  with  study  area  file. 

(b)  Scatter  plot  the  predicted  depths  for  two  dates  or  with  the 
measured  depths. 

(c)  Use  the  previously  calculated  study  set  statistics  to 
select  a  reference  date.  Input  the  reference  date. 

(d)  Set  maximum  K  and  p  parameter  values  and  delta  changes 
which  are  acceptable  in  the  least  squares  parameter 
analysis. 

(e)  Operator  selects  dates  from  the  study  set  whose  water 

depths  are  to  be  regressed  against  those  for  the  reference 
date.  The  slope  and  intercept  are  used  to  modify  K  and  p 
parameters  for  the  individual  dates.  If  adjusted  para¬ 

meters  exceed  operator  designated  limits,  the  operator  may 
eliminate  that  data  set  from  the  analysis  or  reset  the 
parameters.  Adjusted  parameters  are  then  used  to  predict 
water  depths  for  each  pixel  in  the  study  area.  Residuals 
from  the  reference  date  are  computed  for  each  of  the  date 
sets  utilized  in  the  analysis.  The  newly  predicted  depths 
(by  least  squares  parameter  adjustment)  are  then  averaged 
for  each  pixel  to  obtain  a  "best"  estimate.  Residuals  and 
"best"  estimated  water  depths  are  stored  in  the  study  set 
file.  The  adjusted  K  and  p  parameters  and  regression  sta¬ 
tistics  are  stored  in  a  parameter  file  for  the  study  set. 
The  operator  can  obtain  a  printout  of  the  study  area  file 
and/or  the  parameter  file. 

(f)  Same  analysis  as  in  (e)  with  the  reference  date  data 
replaced  by  the  measured  water  depths.  The  subsequent 
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analysis  would  be  performed  only  on  those  pixels  of  the 
study  set  for  which  there  exists  measured  water  depths. 


(g)  Operator  selects  dates  for  analysis  but  no  reference  date 
or  measured  data  are  utilized.  Average  depths  are  computed 
for  each  pixel  for  the  dates  selected.  These  average 
depths  then  become  the  reference  data  set  as  in  option  (e) 
and  an  identical  process  is  completed  to  obtain  a  "best" 
estimator. 

(h)  Operator  inputs  uncertainties  in  original  selection  of  K 
and  p  parameters  to  obtain  an  estimate  of  the  corresponding 
depth  errors  in  the  "best"  depth  estimate.  Errors  are 
computed  for  1.0,  2.0,  3.0,  5.0,  7.0,  10.0,  15.0,  20.0,  and 
30.0  meter  water  depths. 

(4)  APPLY  -  The  "best"  depth  predictor  as  derived  in  (4)  is  applied 
to  the  Landsat  data  to  obtain  estimated  depths  for  the  entire 
subscene.  Operator  initiates  APPLY  command  by  designating  the 
parameter  file  which  contains  subscene  name,  dates,  and  para¬ 
meter  values  needed  to  compute  the  "best"  multi -temporal  esti¬ 
mate.  Results  are  included  in  the  multi-date  file  structure 
for  the  subscene  and  can  be  subsequently  displayed  and  conpared 
with  previous  single  date  predictions  or  other  multi-date  esti¬ 
mates  using  the  LOAD  command.  Other  known  chart  depths  in  the 
subscene  can  be  checked  against  the  "best"  estimated  depth 
using  the  analysis  command  options  (a)  and  (b). 

The  balance  of  the  required  processing,  such  as  obtaining  hardcopy 
of  the  depth  map,  can  be  accomplished  using  existing  DIPS  software.  The 
MTP  software  delivered  and  installed  on  DIPS  consists  of  a  series  of 
modules  which  support  the  mul t i -temporal  command  operations  and 
interface  with  existing  routines. 

Software  has  been  developed  as  described  in  Appendices  A  and  B 
which  can  be  used  to  implement  the  above  procedures  on  the  DIPS.  Sec¬ 
tion  5.0  discusses  suggested  applications  of  this  software  on  the  DIPS. 
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MULT  I -TEMPORAL  DATA  SET  CONSTRUCTION 

Two  processing  steps  were  performed  in  order  to  prepare  the  six 
previously  processed  Landsat  scenes  and  the  1980  Ship  Transect  Data  for 
analysis.  First,  each  of  the  six  Landsat  data  sets  was  registered  to 
geodetic  coordinates  using  ERIM's  semi-rigid  Landsat  model  and  available 
ground  control  points.  Second,  the  1980  ship  transect  data,  referenced 
to  latitude,  longitude  geodetic  coordinates  by  SAI,  were  sampled  and 
merged  with  six  sets  of  raw  Landsat  data  and  six  depth  estimates.  This 

latter  step  created  a  data  base  of  about  400  pixels  for  four  test  areas 

around  the  Great  Bahama  Bank.  The  four  test  areas  used  for  this  study 

are  designated  as  3A,  3B,  3C,  and  3D  as  shown  in  Figure  1.  Test  areas 

were  selected  which  exhibit  a  range  of  water  depths  and/or  bottom  type. 
The  rationale  for  this  selection  is  more  fully  explained  in  section  5. 

4.1  DATA  SET  DESCRIPTION 

The  Landsat  data  sets  we  used  for  analysis  had  been  selected  and 
previously  processed.  Table  1  lists  the  scene  ID's  and  other  relevant 
information.  As  described  in  [1]  the  data  sets  were  processed  for  water 
depths  using  a  combined  ratio-single  band  algorithm  and  the  detector 
parameters  shown  in  Table  2.  For  all  data  sets  a  bottom  reflectance  of 
0.22  was  assumed.  This  value  is  used  along  with  other  parameters  to 

calculate  V  .  For  all  data  sets  water  attenuation  coefficients  of  K,  = 

_1  °  _1  ** 
0.0748m  and  Kg  =  0.326m  were  used,  corresponding  to  values  for 

Jerlov  Type  IB  water. 

Previous  water  depth  processing  results  exhibited  varying  water 
penetration  and  depth  uncertainty  owing  to  seasonal  changes  in  water 
clarity  and  to  cloud  and  haze  patterns.  Because  data  quality  is  impor¬ 
tant  for  subsequent  analyses,  a  qualitative  discussion  is  presented  in 
section  5.2. 
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TABLE 

1.  LIST  OF  1 

LANDSAT  SCENES  OVER 

THE  GRAND  BAHAMA 

BANK 

Satel I ite 

Solar 

Date 

Scene  ID 

Landsat  1  or  2 

Mode 

Elevation 

25  Feb  77 

5678-14102 

1 

Low  gain 

30° 

3  Feb  75 

1925-15015 

1 

Low  gain 

34° 

24  Dec  75 

5249-14435 

1 

High  gain 

28° 

29  Dec  74 

1889-15033 

1 

Low  gain 

30° 

11  Oct  77 

2993-14385 

2 

High  gain 

36° 

25  Jun  77 

2885-14444 

2 

High  gain 

54° 

TABLE  2.  VALUES  OF  VQi  AND  V$  FOR  THE  SIX  LANDSAT  SCENES 

VALUES  OF  VOi 


Scene  Date 

Solar 

V04(MSS4) 

Vn,(MSS5) 

Elevation 

25  Feb 

77 

30° 

22.4 

26.2 

3  Feb 

75 

34° 

25.1 

29.3 

24  Dec 

75 

28° 

63.2 

73.7 

29  Dec 

74 

30° 

22.4 

26.2 

11  Oct 

77 

36° 

88.6 

121.0 

25  Jun 

77 

54° 

122.0 

166.0 

VALUES 

OF  Vs  (MSS4) 

Detector  Number 

Scene  Date 

1 

2 

3  4 

5 

6 

25  Feb 

77 

15.2 

16.0 

15.2  15.2 

15.8 

15.2 

3  Feb 

75 

15.3 

15.1 

15.0  15.1 

15.6 

15.3 

24  Dec 

75 

45.5 

45.8 

46.0  45.8 

46.0 

45.5 

29  Dec 

74 

16.4 

16.1 

16.7  16.6 

16.4 

16.3 

11  Oct 

77 

33.4 

37.4 

38.5  37.4 

41.5 

40.6 

25  Jun 

77 

55.7 

58.0 

63.7  60.4 

64.9 

62.7 
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4.2  1980  SHIP  SURVEY  DATA 

In  July  and  August  1980  a  series  of  cruises  were  made  in  vessels 
operated  by  the  Johns  Hopkins  Applied  Physics  Laboratory.  Data  col¬ 
lected  from  these  vessels  included  echo  sounding  depth  transects,  sub¬ 
mersible  photometer  measurements  in  Landsat  MSS  and  TM  spectral  bands 
and  high  spectral  resolution  bottom  reflectance  data.  Because  of 
previous  difficulties  with  obtaining  reliable  ship  position,  a  special 
emphasis  was  made  in  this  survey  to  gather  accurate  coordinates  for 
survey  sampling  positions.  A  LORAC  positioning  receiver  was  used  in 
connection  with  a  series  of  geodetic  positions  as  located  with  a  satel¬ 
lite  positioning  system.  Details  of  the  reduction  of  the  navigation  and 
echo  sounding  fathometer  data  are  presented  in  ref  [1].  Bottom  reflec¬ 
tance  spectra  collected  with  ah  ISCO  spectral  radiometer  have  been 
previously  reported  [4]. 

All  depth  sounding  and  location  data  were  supplied  to  ERIM  on  mag¬ 
netic  tape.  Because  of  the  high  spatial  density  of  echo  sounding  loca¬ 
tions  relative  to  the  nominal  80  meter  pixel  size  of  Landsat,  the  mea¬ 
surements  for  areas  3A  through  3D  were  sampled  and  averaged  with  Landsat 
pixel  spacing.  The  latitude  and  longitude  of  each  derived  location  was 
then  assigned  to  a  particular  pixel  whose  center  coordinates  were  near¬ 
est  these  values.  With  this  procedure  we  were  able  to  obtain  a  repre¬ 
sentative  water  depth  value  for  each  pixel  which  was  intersected  by  the 
ship  transect. 

4.3  IMAGE-TO- IMAGE  REGISTRATION 

Before  multi -temporal  analysis  could  be  conducted,  each  of  the  six 
Landsat  scene  dates  had  to  be  co-registered  to  one  another.  We  first 
transformed  each  scene  to  geodetic  coordinates  and  resampled  each  scene 
by  nearest  neighbor  resampling.  The  registration  to  geodetic 
coordinates  was  required  to  merge  in-depth  sounding  information. 

The  registration  procedure  uses  a  semi-rigid  Landsat  imaging  model 
and  a  few  (5-10)  well  spaced  control  points  per  Landsat  scene  to  compute 
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two  twenty-two  term  mapping  polynomials.  Control  points  are  required 
because  the  satellite  ephermeris  and  attitude  information,  as  reported 
in  the  SIAT  file  of  the  Landsat  data  tapes,  is  not  sufficiently  precise 
to  assure  subpixel  registration  accuracy  to  a  geodetic  grid. 

Difficulties  were  encountered  in  correcting  the  scenes  because  of 
lack  of  well  spaced  control  points.  Four  good  ground  control  points 
were  located  in  the  Bahamas.  Because  the  western  boundary  of  the  scene 
covered  the  eastern  coast  of  Florida,  additional  ground  control  was 
sought  from  the  1:250,000  scale  Miami  and  West  Palm  Beach  sheets.  These 
control  points,  along  with  additional  points  obtained  from  chart  26320 
(scale  1:300,000)  were  later  rejected  as  being  too  imprecise.  Unfortun¬ 
ately,  this  left  us  with  only  four  points  on  the  eastern  edge  of  the 
scene  and  no  points  on  the  western  edge.  As  a  result,  the  model  east- 
west  errors  are  considerably  larger  (102.2m)  than  the  north-south  errors 
(29.7m).  Table  3  shows  the  results  of  the  modeling  efforts.  Notice 
that,  although  all  control  points  are  listed,  only  those  with  unit 
weight  are  used  in  the  model  application.  Similar  results  were  obtained 
with  other  frames.  The  conclusion  is  that  the  resulting  corrected  data 
set  matches  the  ship  transect  data  to  within  about  two  pixels.  This 
accuracy  should  be  adequate  for  most  bathymetry  analyses  except  in  cases 
where  there  is  an  abrupt  change  in  bottom  depth  or  reflectance. 

Nearest  neighbor  resampling  was  used  to  obtain  the  geometric  cor¬ 
rected  depth  files.  Use  of  cubic  convolution  or  restoration  is  not 
appropriate  for  these  data  since  non-linear  processing  has  been  applied 
to  the  Landsat  radiometric  data  values. 

Data  were  resampled  into  a  Universal  Transverse  Mercator  (UTM)  pro¬ 
jection  with  50m  pixels.  From  this  projection  it  is  possible  to  compute 
the  latitude  and  longitude  of  each  pixel,  using  well  documented  formu¬ 
las.  It  is  also  easy  to  compute  the  pixel  line  and  point  number  from  a 
given  latitude  and  longitude. 
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TABLE  3.  LANDSAT  IMAGING  MODEL  RESULTS 


Landsat  Ground  Control  Points  REV  10.0 
SCENE  ID  25  FEB  77 

RMS  ERRORS  EAST-WEST  102.2  NORTH-SOUTH  29.7  (METERS) 
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4.4  REGISTRATION  OF  SHIP  TRANSECT  AND  IMAGE  DATA 

Smoothed  ship  transect  depth  data  (see  section  4.2)  and  Landsat 
estimates  were  merged.  Data  were  stored  in  a  list  file  for  later  access 
by  the  statistical  analysis  routines  as  discussed  in  section  5. 


PROCEDURE  ANALYSIS  AND  EVALUATION 

As  discussed  in  the  previous  sections*  this  contract  work  involved 
both  the  development  of  a  multi -temporal  processing  procedures  and  its 
implementation  on  the  DIPS.  The  approach  taken  was  to  develop  elements 
of  the  procedure  based  upon  the  DIPS  software  capabilities,  theoretical 
considerations,  and  anticipated  quantity  of  multi-date  imagery.  Proced¬ 
ure  evaluation  was  made  with  a  single  multi -temporal  Landsat  data  set. 
Because  the  software  developed  for  the  DIPS  could  not  be  directly  imple¬ 
mented  on  the  ERIM  PDP/11  the  procedure  evaluation  analyses  were  carried 
out  on  the  University  of  Michigan  MTS  computer  system.  The  approach 
involved  using  the  MTS  statistical  analysis  package  on  the  previously 
registered  ship  and  multi -temporal  data  set  as  discussed  in  section  4.0. 
The  following  sections  describe  some  of  the  statistical  characteristics 
of  this  data  set,  typical  results  obtained  when  these  data  were  used  to 
implement  proposed  procedures,  an  interpretation  of  these  analyses,  and 
based  upon  our  experience,  a  recommended  set  of  initial  applications  of 
the  MTP  software. 

5.1  DATA  SET  QUALITY  AND  NOISE  CHARACTERISTICS 

For  purposes  of  this  description  the  Landsat/ship  data  set  as¬ 
sembled  consisted  of  four  separate  portions,  one  each  from  calibration 
subareas  3A,  3B,  3C,  and  3D.  Each  set  contained  the  multi -date  Landsat 
derived  water  depths  for  25  February  77,  3  February  77,  29  December  74, 
24  December  75,  25  June  77,  and  11  October  77.  In  addition  each  set 
contained  the  survey  ship  measured  soundings  and  TM  radiometer  measure¬ 
ments.  Each  of  these  multi-variate  data  sets  3A,  3B,  3C,  and  3D  con¬ 
tained  respectively  105,  231,  202,  and  238  pixel  locations.  Initially 
it  was  necessary  to  ascertain  the  relative  quality  and  noise  condition 
of  each  of  the  six  independent  water  depths.  Of  the  calibration  areas 
selected  for  this  multi -temporal  analysis,  3D  was  found  to  contain 
regions  with  little  variation  in  measured  water  depth.  For  this  reason 


3D  was  considered  a  good  candidate  for  noise  analysis  of  the  Landsat 
predicted  depths.  A  total  of  twenty  pixels  were  selected  from  an 
approximately  two  square  kilometer  area  within  3D.  For  each  of  these  a 
3x3  array  was  recovered  from  each  of  the  six  dates  with  the  center 
pixel  corresponding  to  that  pixel  selected  from  3D.  For  each  scene  date 
a  local  mean  and  standard  deviation  were  calculated  for  each  array  and 
used  to  estimate  a  standard  deviation  and  mean  for  the  entire  twenty 
arrays.  These  calculations  are  summarized  in  Table  4. 

TABLE  4.  STATISTICS  FOR  LANDSAT  EXTRACTED 
WATER  DEPTHS  SELECTED  FROM  AREA  3D 


Scene 

Date 

Mean 

Depth  (m) 

[x] 

Standard 
Deviation 
W  Ox] 

V* 

oV  /Vs 
(MSS4) 

°V/V 

(MSS4) 

Ratio 
Col.  6/ 
Col.  5 

25  Feb  77 

5.2 

0.82 

0.158 

0.034 

0.049 

1.44 

03  Feb  77 

6.6 

1.18 

0.179 

0.047 

0.069 

1.46 

29  Dec  74 

4.2 

1.11 

0.264 

0.042 

0.070 

1.67 

24  Dec  75 

6.3 

1.57 

0.249 

0.036 

0.083 

2.31 

25  Jun  77 

9.8 

2.50 

0.255 

0.033 

0.121 

3.67 

11  Oct  77 

8.4 

2.74 

0.326 

0.050 

0.167 

3.34 

The 

ship-measured  depths  for 

the  twenty 

pixels 

exhibited  a  mean 

of 

8.58  meters  and  a  standard  deviation  of  0.227  meters.  The  scene  dates 
in  the  table  have  been  placed  in  order  of  quality  from  the  best  to  the 
poorest  based  upon  (1)  visual  inspection  of  the  resulting  water  depth 
maps,  (2)  the  standard  deviation  of  the  predicted  depth,  and  (3)  the 
ratio  of  standard  deviations  in  MSS4  i.e.,  ratio  of  column  6  to  column  5 
as  shown  in  the  table  above.  Column  5  is  the  ratio  of  the  standard 
deviation  in  MSS4  over  deep  water  to  the  mean  deep  water  signal,  V$. 
Column  6  is  the  same  ratio  but  where  the  standard  deviation  and  mean 
MSS4  signal  are  averaged  over  the  twenty  pixel  arrays  in  3D.  The  ratio 
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of  these  two  quantities  (Col.  7)  provides  an  indicator  of  how  the  data 
vary  in  noise  properties  from  those  determined  in  a  deep  water  region. 
Analysis  of  the  TM2  submersible  radiometer  values  in  the  TM  green  band 
for  the  same  twenty  pixels  yield  a  standard  deviation  equivalent  to 
eight  percent  of  the  mean  value.  Some  of  that  change  is  due  to  changes 
in  water  depth  and  subsurface  downwelling  irradiance.  Comparisons  made 
in  the  table  above  suggest  the  following.  While  the  11  Oct  77  predicted 
mean  depth  falls  closest  to  that  measured  by  the  ship,  it  exhibits 
poorest  reliability  because  of  the  large  noise  components  which  are  far 
in  excess  of  those  due  to  changes  in  water  depth,  bottom  type,  or  noise 
associated  with  deep  water  signals.  The  sources  of  this  noise  are 
likely  a  combination  of  errors  from  image  to  image  registration,  pixel 
extraction,  and  atmospheric  conditions.  In  the  first  three  dates  of  the 
table,  on  the  other  hand,  variations  are  only  about  50  percent  greater 
than  those  associated  with  deep  water  variations  in  V$,  and  are  consid¬ 
ered  superior  to  the  last  three.  Since  there  undoubtedly  exist  some 
errors  due  to  each  of  the  previously  mentioned  sources,  the  reported 
standard  deviations  in  the  Landsat  predicted  depth  appear  reasonable 
even  though  they  suggest  one  meter  accuracy  at  eight  meters  depth  if  the 
offsets  are  corrected.  The  large  error  in  predicted  depth  is  due  to 
offset  which  suggest  some  difficulty  in  calculating  representative 
values  of  Vs  and  VQ  from  the  data.  Further  it  underscores  the  need  for 
multi -temporal  analysis  to  resolve  observed  differences  in  predicted 
depth  and  provide  a  best  estimate  water  depth. 

5.2  EXAMPLE  CALCULATIONS  USING  MULTI-TEMPORAL  PROCEDURES 

Example  calculations  using  suggested  multi-temporal  procedures  are 
presented  here  to  show  methods  of  operation  and  value  when  applied  to 
the  assembled  Landsat/measured  data  sets.  Basically,  calculations  were 
made  with  and  without  the  aid  of  ship  measured  water  depths.  The  large 
errors  in  mean  water  depth  shown  above  for  a  portion  of  area  3D  were 
found  with  each  of  the  other  areas  as  well.  This  finding  confirms  that 
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unless  ship  surface  truth  is  used  in  the  depth  analysis,  large  errors 
due  to  offsets  in  depth/signal  relationships  can  be  expected.  Statisti¬ 
cal  and  depth  analysis  results  are  summarized  for  areas  3A,  3B,  3C,  and 
3D  in  Tables  5  through  8  respectively.  Of  the  four  areas,  3B  showed  the 
best  distribution  of  water  depths  and,  therefore,  potentially  the 
greatest  opportunity  for  resolution  of  differences  in  predicted  and 
measured  water  depths.  The  following  discussion  is,  therefore  denoted, 
primarily  to  the  results  obtained  from  area  3B.  The  noise  analysis  con¬ 
ducted  above  indicated  data  from  three  of  the  six  scene  dates  to  be 
suitable  for  multi -temporal  processing.  As  a  result  most  calculations 
reported  were  made  using  these  three  dates.  For  purposes  of  this  analy¬ 
sis  the  231  pixels  extracted  from  3B  were  divided  into  two  groups,  1-75 
and  76-231.  The  first  served  as  a  calibration  set  and  the  latter  as  a 
validation  or  test  set.  In  general  two  types  of  multi -temporal  analyses 
were  performed  on  the  3B  data  set  --  (1)  A  least  squares  adjustment  of  K 
and  p  parameters  between  scene  dates  and  between  an  adjusted  average  and 
the  measured  water  depths.  (2)  A  straight  average  of  the  independent 
satellite  predicted  depths.  Initially  data  from  the  best  three  dates 
were  used  to  resolve  K  and  p  parameter  differences  using  date-to-date 
regression.  For  these  cases  all  231  pixels  were  used  to  produce  the 
regression  equations  (6.1)  and  (6.2)  in  Table  6. 

The  resulting  coefficients  and  constants  indicate  the  amount  of 

adjustment  in  K  and  VQ  necessary  to  bring  the  two  data  sets  into 

agreement  in  a  least  squares  manner.  The  coefficient  (a^)  dictates  the 

amount  of  adjustment  in  K  necessary  to  bring  the  two  data  sets  into 

agreement  (K*  =  a.j~*  *  K).  The  constant,  b^ ,  in  combination  with  the 

coefficient,  determines  the  adjustment  necessary  in  the  V_  term  (V  *  = 
2Kb  / a  )  oo 

VQ*e  v  v.  If  the  ratio  method  is  used  we  are  referring  to  a  K 

difference  and  a  ratio  of  VQ  for  MSS4  and  MSS5.  Using  equations  6.1  and 

6.2,  data  sets  2  and  3  can  be  transformed  to  estimate  V^.  The  remaining 

differences  can  be  attributed  to  random  noise  processes.  The  random 

differences  can  be  reduced  by  averaging  (V^),  and  (V^).  The 


TABLE  5.  SUMMARY  STATISTICS  AND  REGRESSION  RESULTS 
FOR  CALIBRATION  AREA  3A 


Pi  xel 

No.  51- 

■100 

All  Pixels 

VARIABLE 

MIN 

MAX 

MEAN 

STD  DEV 

MIN 

MAX 

MEAN 

STD  DEV 

Predicted  Depth  (m) 
February  25,  1977 

2.50 

25.50 

6.49 

3.90 

1.90 

25.50 

7.04 

5.17 

Predicted  Depth  (m) 
February  3,  1977 

3.80 

25.50 

7.08 

5.03 

3.70 

25.50 

7.26 

5.30 

Predicted  Depth  (m) 
December  29,  1974 

0.30 

25.50 

6.83 

3.37 

0.30 

25.50 

7.09 

3.61 

Snip  Meas.  Depth(m) 

5.34 

29.54 

9.15 

3.59 

5.34 

35.93 

9.57 

4.10 

Predicted  Depth  (m) 
from  Eq.  3  below 

6.20 

23.23 

9.15 

2.89 

5.75 

23.23 

9.56 

3.83 

Residual  Errors  (m) 
Eq.  3  below 

-6.52 

6.31 

0.10 

2.16  - 

13.00 

12.70 

0.00 

3.02 

Average  Adjusted  Depth  (m) 

from  Three  Dates 

4.18 

23.93 

6.58 

3.33 

4.10 

23.93 

6.90 

3.98 

Predicted  Depths  (m) 
from  Eq.  4  below 

7.42 

23.04 

9.32 

2.63 

7.36 

23.04 

9.57 

3.15 

Residual  Errors  (m) 
from  Eq.  4  below 

-7.75 

6.93 

-0.17 

2.35 

-7.75 

12.89 

O.OU 

2.65 

Correlation 

Standard 

Regression  Equation 

Coefficient 

Error 

_ 

5.1  V1(V2)  =  0.885 

v2  +  0. 

611 

0.909 

2. 

16 

5.2  Vj(V3)  =  1.066 

V3  -  0. 

510 

0.745 

3. 

46 

5.3  z  =  0.740  V.  + 

4.35 

0.804 

2. 

16 

5.4  7  =  0.822  (Adjusted  Average) 

+  3.74 

0.762 

2. 

35 

where: 

Vj  =  Predicted  Depth  (m)  25  FEB  77;  V2  =  Predicted  Depth  (m)  3  FEB  77; 
Vj  =  Predicted  Depth  (m)  29  DEC  74;  z  =  Actual /Measured  Depth  (m) 
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TABLE  6.  SUMMARY  STATISTICS  AND  REGRESSION  RESULTS 
FOR  CALIBRATION  AREA  3B 


Pixel 

No.  1- 

75 

Pixel 

No.  76- 

-231 

VARIABLE 

MIN 

MAX 

MEAN 

STD  DEV 

MIN 

MAX 

MEAN 

STD  DEV 

Predicted  Depth  (m) 
February  25,  1977 

1.80 

25.50 

8.05 

6.08 

2.70 

25.50 

10.99 

5.83 

Predicted  Depth  (m) 
February  3,  1977 

1.90 

25.50 

9.69 

6.68 

3.10 

25.50 

12.82 

6.92 

Predicted  Depth  (m) 
December  29,  1974 

3.10 

18.20 

8.11 

4.26 

3.20 

25.50 

10.77 

5.40 

Ship  Meas.  Depth(m) 

3.88 

39.36 

10.81 

6.93 

4.78 

43.33 

13.36 

7.73 

Predicted  Depth  (m) 
from  Eq.  3  below 

4.97 

27.10 

10.81 

5.68 

5.81 

27.10 

13.55 

5.44 

Residual  Errors  (m) 
from  Eq.  3  below 

-10.90 

20.76 

0.00 

3.99  - 

12.14 

16.23 

-0.19 

4.92 

Average  Adjusted  Depth  (m ) 

from  Three  Dates 

2.53 

21.44 

8.24 

5.03 

3.67 

23.88 

10.90 

5.23 

Predicted  Depths  (m) 
from  Eq.  4  below 

4.37 

25.71 

10.81 

5.68 

5.66 

28.47 

13.82 

5.90 

Residual  Errors  (m) 
from  Eq.  4  below 

-10.44 

17.07 

0.00 

2.38  - 

10.08 

14.91 

-0.46 

4.21 

Correlation 

Standard 

Regression  Equation 

Coefficient 

Error 

—  ■ 

6.1  V1 ( V2)  =  0.763 

v2  +  1. 

031 

0.879 

2 

.89 

6.2  V1(V3)  =  1.001 

v3  +  0. 

113 

0.859 

3 

.11 

6.3  z  =  0.933  Vj  + 

3.29 

0.820 

3 

.99 

6.4  z  =  1.129  (Adj 

usted  Average) 

+  1.51 

0.820 

2 

.38 

where: 

VT  =  Predicted  Depth  (m)  25  FEB  77;  V?  =  Predicted  Depth  (m)  3  FEB  77; 
V.  =  Predicted  Depth  (m)  29  DEC  74;  z  =  Actual /Measured  Depth  (m) 
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TABLE  7.  SUMMARY  STATISTICS  AND  REGRESSION  RESULTS 
FOR  CALIBRATION  AREA  3C 


Pixel 

NO.  1- 

50 

Pi  xel 

No.  51- 

202 

VARIABLE 

MIN 

MAX 

MEAN 

STD  DEV 

MIN 

MAX 

MEAN 

STD  DEV 

Predicted  Depth  (m) 
February  25,  1977 

1.90 

25.50 

6.54 

5.41 

2.30 

25.50 

6.75 

4.21 

Predicted  Depth  (m) 
February  3,  1977 

0.60 

25.50 

6.56 

4.99 

2.70 

25.50 

7.78 

4.60 

Predicted  Depth  (m) 
December  29,  1974 

1.20 

25.50 

5.37 

4.53 

2.20 

25.50 

6.48 

3.38 

Ship  Meas.  Depth(m) 

4.86 

35.05 

9.81 

4.97 

6.39 

24.90 

9.53 

2.55 

Predicted  Depth  (m) 
from  Eq.  3  below 

6.07 

25.08 

9.81 

4.36 

6.40 

25.08 

9.98 

3.39 

Residual  Errors  (m) 
from  Eq.  3  below 

-7.63 

9.97 

0.00 

2.42 

6.40 

3.49 

0.45 

1.91 

Average  Adjusted  Depth  (m) 
from  Three  Dates  1.76 

25.29 

6.07 

4.81 

2.84 

25.28 

6.90 

3.84 

Predicted  Depths  (m) 
from  Eq.  4  below 

5.84 

27.52 

9.81 

4.44 

6.83 

27.52 

10.58 

3.53 

Residual  Errors  (m) 
from  Eq.  4  below 

-5.01 

7.53 

0.00 

2.28 

6.62 

2.79 

-1.05 

1.91 

Regression  Equation 

Correlation 

Coefficient 

Standard 

Error 

7.1  V1(V2)  =  0.876 

v2  +  0. 

11503 

0.914 

1. 

84 

7.2  V  3 ( V  3 )  =  1.098 

v3  -  0. 

1125 

0.902 

1. 

96 

7.3  z  =  0.805  V,  + 

4.55 

0.876 

2. 

42 

7.4  z  =  0.921  (adjusted  average) 

+  4.22 

0.892 

2. 

28 

where ; 


V1  =  Predicted  Depth  (m)  25  FEB  77;  V2  =  Predicted  Depth  (m)  3  FEB  77; 
V3  -  Predicted  Depth  (m)  29  DEC  74;  z  =  Actual /Measured  Depth  (m) 


^ERJ _ 

TABLE  8.  SUMMARY  STATISTICS  AND  REGRESSION  RESULTS 
FOR  CALIBRATION  AREA  3D 

Pixel  No.  150-224  Pixel  No.  1-149,  225-238 

VARIABLE  MIN  MAX  MEAN  STD  DEV  MIN  MAX  MEAN  STD  DEV 

Predicted  Depth  (m)  1.90  7.50  3.76  1.20  1.90  6.60  4.63  1.08 

February  25,  1977 

Predicted  Depth  (m)  2.70  9.90  4.64  1.35  2.70  8.50  5.34  1.23 

February  3,  1977 

Predicted  Depth  (m)  1.60  8.50  3.82  1.23  1.20  6.50  3.98  0.94 

December  29,  1974 

Ship  Meas.  Depth(m)  2.35  9.58  5.48  2.60  2.50  9.72  7.56  1.69 

Predicted  Depth  (m)  3.85  8.76  5.48  1.05  3.85  7.97  6.24  0.95 

from  Eq.  3  below 

Residual  Errors  (m)  -3.50  4.94  0.00  2.39  -2.58  4.57  1.31  1.52 

from  Eq.  3  below 

Average  Adjusted  Depth  (m) 

from  Three  Dates  2.76  6.93  4.04  0.83  2.49  6.09  4.50  0.68 

Predicted  Depths  (m) 

from  Eq.  4  below  3.51  9.92  5.48  1.28  3.09  8.63  6.18  1.04 

Residual  Errors  (m) 

from  Eq.  4  below  -3.14  5.18  0.00  2.28  -2.83  5.12  1.38  1.57 

Correlation  Standard 

Regression  Equation  Coefficient  Error _ 

8.1  V1(V2)  =  0.571  V2  +  1.436  0.628  0.925 

8.2  Vj(V3)  =  0.650  V3  +  1.800  0.566  0.980 

8.3  z  =  0.876  Vj  +  2.19  0.404  2.39 

8.4  z  =  1.537  (Adjusted  Average)  -  0.736  0.491  2.28 

where: 

Vj  =  Predicted  Depth  (m)  25  FEB  77;  V2  =  Predicted  Depth  (m)  3  FEB  77; 

V^  =  Predicted  Depth  (m)  29  DEC  74;  z  =  Actua 1 /Measured  Depth  (m) 
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result,  which  we  will  refer  to  as  (adjusted  average),  can  then  be 
related  to  the  measured  depths  in  an  effort  to  make  final  adjustments  to 
K  and  p  parameters  as  shown  in  Figure  3.  For  this  case  the  regression 
equation  is  given  as  equation  6.4  in  Table  6. 

Predicted  depths  are  shown  for  both  the  regression  pixels  (1-75) 
and  the  balance  of  3B  (76-231).  The  residual  errors  between  this  model 
and  the  ship  measurements  are  shown  in  Figure  4.  The  residual  patterns 
appear  random  except  for  groups  of  points  along  parallel  lines  oriented 
at  a  sixty  degree  slope.  These  residual  patterns  are  associated  with 
the  use  of  quantized  signal  levels  used  to  predict  discrete  depths 
rather  than  continuous  levels.  This  effect  will  be  most  pronounced  in 
deeper  waters  where  there  are  just  a  few  Landsat  raw  data  count  changes 
over  a  large  range  of  depths.  In  these  cases  a  single  depth  is  pre¬ 
dicted  for  pixels  having  a  range  of  measured  depths.  The  residual  is 
a  simple  linear  function  of  the  measured  depth.  Outside  of  these  pat¬ 
terns,  the  residuals  appear  to  be  random.  When  this  analysis  process 
was  repeated  using  all  six  dates,  the  standard  error  of  the  estimate  in¬ 
creased  from  2.38  to  5.27  meters.  This  increase  is  expected,  given  that 
the  latter  three  scene  dates  are  of  relatively  poorer  quality  as  discus¬ 
sed  in  section  5.1.  A  further  comparative  analysis  was  made  by  using 
only  the  first  and  best  scene  date  (25  Feb  77).  The  resulting  regression 
equation  is  given  as  equation  6.3  in  Table  6.  Plots  of  this  regression 
analysis  and  residual  errors  are  shown  in  Figures  5  and  6  respectively. 
The  standard  error  of  3.99  meters  is  approximately  3  times  that  obtained 
for  the  three  date  case  above.  This  comparison  suggests  that  the  pri¬ 
mary  effect  of  using  the  multiple  dates  was  simple  reduction  of  random 
noise. 

In  the  second  type  of  analysis  performed  on  these  data,  Landsat 
predicted  depths  were  simply  averaged  on  a  pixel  by  pixel  basis  with  no 
parameter  adjustments  from  those  originally  assumed.  Averages  for  three 
and  six  date  cases  are  plotted  against  measured  ship  depths  as  well  as 
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Figure  4.  Scatter  Plot  of  Residual  Errors  of  the 

Regression  Estimate  (see  Figure  3)  versus 
the  Measured  Depth. 
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depth  differences  in  Figures  7  through  10.  Correlation  coefficients  of 
0.840  and  0.855  were  obtained  between  the  three-date  and  six-date 
averages  and  the  measured  ship  depths  respectively.  The  six  date 
average  produced  a  slightly  tighter  grouping  of  plotted  points  in  the  5 
to  15  meter  depth  range  than  the  three  date  case  indicating  this  average 
was  less  sensitive  to  actual  depth  changes.  Thus  while  the  straight 
averaging  process  will  tend  to  reduce  absolute  errors  in  mean  depth  as 
described  above,  the  average  as  a  best  predicted  depth  will  show  greater 
absolute  errors  as  water  depths  deviate  from  the  mean.  If  one  can 
effectively  reduce  the  constant  differences  between  predicted  and 
measured  depths  then  the  adjusted  averaging  process  appears  to  be 
superior  to  straight  averaging  of  multidate  extracted  depths. 

5.3  MULTI-TEMPORAL  ANALYSIS  PROCEDURES  ON  THE  DIPS 

The  results  described  above  should  not  be  construed  as  an  evalua¬ 
tion  of  the  MTP  software  capabilities;  rather  they  are  results  obtained 
with  MTP  type  operations  which  were  considered  appropriate  to  the  avail¬ 
able  data  set.  As  previously  stated,  evaluation  must  be  made  on  the 
basis  of  analysis  of  several  sets  of  multi-temporal  remote  sensing  data. 
The  following  descriptions  are  intended  as  representative  analysis  pro¬ 
cedures  which  could  be  carried  out  on  the  DIPS  with  the  aid  of  the  ERIM 
developed  software. 

In  each  of  the  following  examples  it  is  assumed  that  necessary  pre¬ 
paration  procedures  (image  to  image  warping,  etc.)  have  been  carried  out 
as  described  in  section  3.4  and  the  DIPS  operator  manuals.  For  each  date 
the  best  available  parameter  estimates  have  been  used  with  the  DIPS 
DEPTH  routine  to  convert  the  Landsat  signal  levels  at  each  pixel  in  a 
512  x  512  subarea  to  an  estimated  water  depth.  At  this  point  the  opera¬ 
tor  can  call  up  the  MTP  menu  and  have  the  following  selection  of 
operations: 


u.P  DEPTHdl) 


Scatter  Plot  of  the  Residual  of  the  Unadjusted  Average 
Predicted  Depth  versus  the  Measured  Depth  (see  Figure  7) 
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Scatter  Plot  of  the  Unadjusted  Average  Landsat 
Predicted  Water  Depth  for  All  Six  Available  Dates 
Versus  the  Measured  Depth. 
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Multi-temporal  Program  Depth  Analysis  Sub-Menu 

1.  Load  image  files  to  display. 

2.  Select  image  plane  for  viewing. 

3.  Conduct  simple  operations  between  image  planes. 

4.  Select  smoothing  operation. 

5.  Select  polygon  areas. 

6.  Make  scatterplots  and  calculate  regression  statistics  for 
pixels  within  polygon(s). 

7.  Extract  statistics  from  polygon  area(s). 

8.  Apply  regression  coefficients  to  adjust  depth  algorithm 
parameters. 

9.  Estimate  parameter  error  propagation. 

10.  Exit  menu. 


The  following  four  examples  are  considered  representative  of  the 
types  of  water  depth  problems  for  which  the  MTP  software  could  help  to 
enhance  the  depth  estimates.  In  each  case  the  objective  is  to  extract 
the  water  depth  information  from  remote  sensor  data  given  available 
surface  truth  measurements  and  water  calibration  depths. 

(1)  Assume  we  are  examining  a  small  area  with  uniform  bottom  type 
but  with  unknown  reflectance.  Suppose  further  there  exists  a 
large  K  difference  between  the  available  scene  dates  as 
ascertained  with  the  SCATTERPLOT  routine.  Use  regression  to 
estimate  the  K  difference  pairwise  for  the  scene  dates.  Return 
to  the  DEPTH  routine  and  treat  the  pairs  of  dates  as  pairs  of 
wavelengths  in  the  ratio  algorithm  which  eliminates  the  need 
for  a  specific  bottom  reflectance  coefficient. 

(2)  Assume  that  the  K  value  is  constant  but  unknown  for  the  avail¬ 
able  scene  dates.  Suppose  that  for  at  least  two  of  the  dates 
there  exists  a  known  tidal  state  change.  Use  the  SCATTERPLOT 
and  REGRESSION  routine  to  confirm  the  offset  in  predicted  water 


depth  due  to  tidal  state.  Use  the  regression  slope  information 
and  the  APPLY  routine  to  remove  any  differences  in  K  type  para¬ 
meters  from  the  data.  Use  the  APPLY  routine  to  calculate  the 
predicted  depth  difference  between  scenes  (pairwise)  and  the 
STATISTICS  EXTRACTION  routine  to  compute  the  average  differ¬ 
ence.  Adjust  the  K  parameter  by  the  ratio  of  this  average  to 
the  known  tidal  state  change  and  recalculate  using  the  DEPTH 
routine.  Actually  we  do  not  have  to  assume  a  constant  K  from 
scene  to  scene  since  the  differences  can  be  assessed  from  the 
REGRESSION  coefficients  and  adjusted  individually  by  the  tidal 
state  factor. 

(3)  Assume  we  are  examining  a  large  area  with  variable  bottom  re¬ 
flectance.  Use  the  SCATTERPLOT,  REGRESSION,  and  APPLY  routines 
to  adjust  out  scene-wide  K  and  r  type  parameter  differences. 
Group  the  available  scene  dates  into  two  groups  according  to 
known  high  or  low  tidal  state.  Adjust  each  scene  date  using 
the  APPLY  routine  to  add  (or  subtract)  a  constant  from  each 
pixel  depth  so  as  to  transform  the  data  to  a  state  of  normal 
high  or  normal  low  tidal  state.  Even  though  the  residuals 
appear  to  be  random  errors  about  the  tidal  difference  they  may 
contain  systematic  spatial  components  due  to  bottom  reflectance 
variations.  A  depth  difference  map(s)  can  be  computed  using 
the  APPLY  routine  and  subsequently  loaded  into  the  available 
Comtal  image  planes  using  the  LOAD  option.  If  the  difference 
maps  display  patterns  which  correlate  with  bottom  features  as 
determined  from  aerial  photos,  ship  surveys,  or  knowledge  of 
coastal  processes,  then  such  difference  maps  can  be  used  with 
the  MTP  software  to  essentially  adjust  the  bottom  reflectance 
on  a  pixel  by  pixel  basis.  The  depth  difference  maps  should 
have  similar  features  and  are  essentially  equivalent  since  the 
scenes  have  been  normalized  for  K  and  tidal  state  differences. 
Random  features  in  these  difference  maps  suggest  that  the 
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influence  of  spatially  varying  noise  due  to  atmospheric  condi¬ 
tions  is  evident.  In  the  former  case,  bottom  reflectance  ad¬ 
justments  can  be  accomplished  as  follows.  First  the  depth 
difference  maps  (high  minus  low  tide)  should  be  smoothed  and 
averaged  (for  usable  scene  dates)  to  reduce  random  noise  com¬ 
ponents.  Pixels  in  the  resulting  difference  map  which  have 
value  greater  than  the  average  (normal  high  tide  minus  normal 
low)  suggest  a  bottom  reflectance  which  is  greater  than  that 
assumed  in  the  original  calculations.  Alternatively  those  with 
lower  value  suggest  a  lower  value  in  bottom  reflectance.  Depth 
deviations  due  to  bottom  reflectance  are  calculated  by  simply 
subtracting  the  mean  difference  (using  APPLY).  These  deviations 
can  then  be  subtracted  from  the  individual  scenes  to  remove  the 
unwanted  effects  due  to  bottom  reflectance  variations. 

(4)  For  this  case  surface  truth  measurements  are  available  to  cali¬ 
brate  water  depths  extracted  from  remote  sensing  data.  First 
the  SCATTERPLOT  and  REGRESSION  routines  are  used  to  analyze  the 
multitemporal  data  set  as  in  the  previous  examples.  Having 
done  so,  the  APPLY  routine  is  used  to  adjust  date-to-date  dif¬ 
ferences  due  to  algorithm  parameter  variations.  At  this  point 
the  multidate  sets  are  essentially  equivalent  and  any  remaining 
date  to  date  differences  are  likely  due  to  random  noise  compon¬ 
ents.  The  best  depth  estimate,  in  this  case,  is  a  pixel-by¬ 
pixel  average  over  the  available  scene  dates.  This  latter 
estimate  is  the  best  one  can  do  without  supporting  surface 
truth  calibration  data.  The  extent  to  which  such  truth  data 
can  be  used  to  improve  the  predicted  depths  depends  on  its 
quality  and  applicability.  The  level  of  representativeness 
dictates  the  spatial  area(s)  of  the  subarea  where  calibration 
depths  can  be  used  to  make  wther  adjustments  to  the  remote 
sensor  extracted  water  depths.  Three  types  of  conditions  seem 
important.  (1)  If  the  measured  depths  are  evenly  distributed 
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across  the  subarea  and  if  they  cover  a  sufficiently  wide  range 
of  depths  then,  resulting  calibrated  predictions  can  be  applied 
to  the  entire  subarea.  (2)  If  on  the  other  hand  the  measured 
depths  are  taken  from  a  single  small  area  then  it  is  doubtful 
that  they  could  be  applied  elsewhere  in  the  subarea.  (3)  If 
multiple  calibration  areas  are  used  and  each  representative  of 
a  separate  bottom  type,  it  may  be  possible  to  iterate  over 
values  of  bottom  reflectance  to  obtain  a  suitable  fit  for  each 
calibration  area.  Depth  algorithms  calibrated  in  this  way 
could  be  applied  to  other  locations  in  the  subarea  with  similar 
bottom  types. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  project  discussed  in  the  report  documents  ERIM's  effort  to 
develop  a  multi -temporal  processing  algorithm  for  DMA's  Digital  Image 
Processing  System  (DIPS).  The  project  utilized  a  Landsat  Multidate  data 
set  which  includes  the  Bahamas  Photobathymetric  Calibration  Area.  The 
processing  algorithm  developed,  however,  is  not  dependent  on  the  speci¬ 
fic  use  of  Landsat  data  but  rather  can  be  applied  in  principle  to  any 
mul titemporal  data  set.  Development  and  evaluation  of  this  algorithm  as 
discussed  in  section  5.0  and  elsewhere  has  led  to  the  following  conclu¬ 
sions  and  recommendations. 

6.1  CONCLUSIONS 

(1)  The  multi-temporal  algorithm  developed  under  this  contract  and 
its  accompanying  software  could  not  be  fully  implemented  on  the  ERIM 
PCP/11  computer  because  of  critical  differences  in  the  DIPS  hardware  and 
software.  This  situation  precluded  any  multi-temporal  image  processing. 
In  view  of  this  situation  a  data  set  was  assembled  consisting  of  coinci¬ 
dent  ship  survey  data  and  multi-date  Landsat  signal  values. 

(2)  The  available  Bahamas  Landsat  data  was  found  to  vary  widely 
from  data  to  date  in  terms  of  its  utility  for  water  depth  extraction  and 
for  multi-temporal  processing.  The  observed  variation  is  considered  to 
be  due  principally  to  atmospheric  and  system  noise.  Of  the  six  avail¬ 
able  data  sets,  three  were  found  to  be  of  comparably  good  quality  and 
three  of  relatively  poor  quality.  When  we  attempted  to  utilize  any  of 
these  latter  scene  dates,  the  predicted  depths  were  less  reliable  when 
compared  to  measured  ship  survey  soundings.  In  view  of  this  experience 
it  is  concluded  that  preliminary  quality  review  and  noise  analysis  must 
accompany  the  selection  of  comparable  multi -date  imagery.  Further  this 
experience  suggests  that  a  typical  multi-date  Landsat  data  set  will 
consist  of  two  or  three  scenes. 
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(3)  During  the  present  study,  three  multi-temporal  depth  procedures 
were  implemented.  In  the  first  procedure,  ship  transect  data  was  used 
along  with  data  from  the  three  best  dates  to  obtain  an  averaged  depth 
file  by  first  adjusting  parameters  of  the  depth  algorithm  for  each  date 
to  minimize  the  mean  square  errors  between  calculated  and  ship  data, 
then  averaging  the  three  revised  depth  files.  In  the  second  procedure, 
the  depth  files  from  the  three  best  dates  were  adjusted,  then  the  files 
averaged.  Results  were  compared  to  ship  data.  In  the  third  procedure, 
depth  files  from  the  best  three  dates  were  simply  averaged,  with  no 
parameter  adjustments,  and  results  compared  to  the  ship  data. 

Results  were  evaluated  by  assessing  the  bias  (mean  error)  and 
standard  error  (mean  square  error)  between  the  ship  data  and  the  result¬ 
ant  average  depth  files.  Results  are  shown  in  Table  9.  For  the  first 
procedure,  the  bias  is  identically  zero  as  a  result  of  the  least  squares 
regression  normalization.  The  standard  error,  for  the  points  examined, 
is  2.38  m.  Because  this  error  depends  on  depth,  the  standard  error  for 
a  set  of  points  different  from  those  used  for  this  analysis  will  gener¬ 
ally  not  be  the  same  as  the  standard  error  we  obtained.  For  procedure 
2,  there  was  a  bias  of  1.51  m  and  the  same  standard  error,  2.38  m.  Bias 
occurs  because  ship  data  were  not  used  in  the  normalization  procedure. 
Procedure  3  produced  a  bias  of  1.94  m  and  a  standard  error  of  4.13  m. 
These  numbers  are  poorer  than  for  procedure  2  because  no  parameter  nor¬ 
malization  was  performed.  The  differences  between  the  biases  and 
standard  errors  of  procedure  2  and  procedure  3  are  an  indication  of  the 
improvement  brought  about  by  parameter  normalization.  The  increase  in 
bias  of  procedure  3  when  going  from  three  scenes  to  six  may  be  a 
reflection  on  the  poorer  data  quality  of  the  additional  three  scenes. 

Another  evaluation  of  procedure  1  was  to  compare  standard  errors  of 
one,  two,  and  three  data  average  depth  files.  Results  are  shown  in 
Table  10.  The  fact  that  the  three  data  standard  error  is  lower  than  the 
one  date  standard  error  shows  the  improvement  to  be  obtained  using  the 
multi-temporal  procedure.  The  results  for  two  dates  appears  anomalous, 
but  few  definitive  conclusions  can  be  drawn  based  on  this  data  set 
alone. 
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COMPARISON  OF  BIAS  AND  STANDARD  ERROR  FOR 
THREE  MULTI-TEMPORAL  DEPTH  TECHNIQUES 


Procedure 

Bias  (m) 

Standard  Error  (m 

(1)  Parameter 

normal i zati on  with  ship  data  -  3  scene 

0 

2.38 

(2)  Parameter 

normal i zati on  without  ship  data  - 
3  scene 

1.51 

2.38 

(3)  Straight  average 

no  normalization  -  3  scene 

1.94 

4.13 

-  6  scene 

2.71 

4.13 

TABLE  10 

COMPARISON  OF  PROCEDURE  1  USING 
ONE,  TWO,  AND  THREE  DATE  DEPTH  DATA 

Number  of  Dates  Standard  Error  (m) 

1  2.89 

2  3.1 

3  2.38 

(4)  Because  of  the  complex  character  of  the  Bahamas  multi -date 
Landsat  data  set  it  is  not  possible  to  predict  general  performance  of 
the  algorithm  for  other  such  data  sets  and  for  other  types  of  multi - 
temporal  data  from  which  bathymetric/hydrographic  information  could  be 
extracted. 

(5)  The  multi-date  algorithm  has  been  designed  with  flexibility  of 
precise  procedure  to  allow  the  DIPS  operator  to  investigate  various 
processing  procedures  to  enhance  not  only  the  accuracy  of  water  depth 
predictions  but  also  the  image  detection  of  submerged  hazards. 


(6)  The  application  of  Kalman  filtering  theory  was  briefly  investi¬ 
gated  as  a  basis  for  a  multi -temporal  algorithm.  The  Kalman  theory 
presents  a  very  generalized  least  squares  formulation  adaptable  to 
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imperfect  parameter  information  and  various  Gaussian  noise  variables. 
This  approach  was,  however,  considered  infeasible  because  of  a  practical 
requirement  for  large  number  of  scene  dates.  Because  of  the  expected 
limited  number  of  dates  available  for  any  one  scene,  it  is  concluded 
that  the  multi -temporal  algorithm  must  be  so  formulated  to  rely  more 
heavily  on  the  spatial  variations  within  any  one  scene  and  less  on  the 
actual  date-to-date  variations  for  any  scene  location  (pixel). 

6.2  RECOMMENDATIONS 

(1)  The  algorithm  developed  for  multi -temporal  remote  sensing  data 
should  be  evaluated  against  other  Landsat  multi -date  sets  and 
also  those  that  may  be  derived  from  aircraft  and  other  sources 
of  high  resolution  spatial  information  systems. 

(2)  Research  and  development  should  be  initiated  to  construct  an 
algorithm  which  normalizes  satellite  and  aircraft  radiometric 
data  on  a  pixel -by-pixel  basis  so  as  to  extend  the  applicabi¬ 
lity  of  water  depth  predictions  beyond  the  immediate  area  of 
surface  truth. 

(3)  Since  the  software  delivered  to  DMA/HTC  has  not  been  completely 
checked  out  we  recommend  DMA  staff,  familiar  with  DIPS,  to 
initiate  a  test  using  the  Bahamas  data  set.  Documentation 
files  provided  with  the  delivered  software  are  sufficient  to 
allow  installation  and  operation. 

(4)  Since  DMA  has  a  requirement  to  upgrade  the  DIPS  as  improved  and 
special  purpose  algorithms  are  developed,  it  would  be  advan¬ 
tageous  to  have  a  DIPS  simulator  on  the  ERIM  PDP/1 1  to  provide 
a  means  for  complete  checkout  of  future  software  and  to  allow 
development  of  special  options  to  operator  processing 
procedures. 
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APPENDIX  A 


SOFTWARE  DESCRIPTION  AND 
INSTALLATION  INSTRUCTIONS 

This  Appendix  contains  listings  of  documentation  files  as  provided 
on  the  ERIM  generated  magnetic  tape.  These  listings  include  INSTALL.DOC 
which  provides  detailed  instructions  to  DMA  DIPS  operators  on  the  proper 
installation  of  the  ERIM  MTP  software.  Also  included  are  documentation 
files  describing  the  overall  MTP  software  (OVERVIEW.DOC),  a  checkout 
procedure  (CHECKOUT.DOC),  and  a  sequence  for  running  the  various  MTP 
menu  options  (RUNNING.DOC). 
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r* aomeist.ooc  ,HH  prim  nov.  id*? 

thjs  "PPAOmfi  ST.nno"  file  IS  THE  first  DOCUMENTATION  FILE  YOU  SHOULD 
READ  EOS  the  E»I«  mTP  mulTI-TEmPORAL  PROCESSING  DEPTH  ALGORITHM  SOFTWARE , 
er-p  TMSTji.L*  Trov  SEE  INSTALL. one  AND  INSTALL,  oho.  OTHERWISE,  USE  README1ST.DOC 
TC  nlPf-T  f O'l  TO  ALL  THE  OTHER  DOCUMENTATION,  AS  FOLLOWS* 

U  install. noe  provides  documentation  and  a  command  file  to 
install. CMD  INSTALL  the  ERIM  SOFTWARE  AUTOMATICALLY  and 
painlessly  (well,  almost,  we  H0PE5, 

?•)  OVERVIEW, oor  PROVIDES  AN  OVERVlFw  OF  THE  ERIM  MTP  SYSTEM. 

S'  CHECKOUT, OOC  --  DISCUSSES  USE  OF  SAMPLE  DATA  FOR  A  RUNTHROU C.H  TO 

SEE  IF  YOU  OFT  THE  SAME  RESULTS  WE  DO.  THE  FOLLOWING 
RUNNING, DOC  SHOULD  HE  P P A 0  IN  CONJUNCTION  WITH  THIS. 

a  )  RUNNING, DOC  --  DESCRIBES  GENERALLY  WHA*  YOU  WANT  TO  DO  WITH  THT 
VARIOUS  PARTS  OF  THE  MTP  SYSTEM  FROM  AN  OVERALL 
ALGORITHMIC  VIEWPOINT,  and  DESCRIBES  THE 
SIGNIFICANT  OPTIONS  AVAILABLE  in  The  PROGRAMS. 

IT  WILL  BE  USEFUL  TO  READ  CHSCKOUT.DOC  WITH  This. 


TN  APfUTTON  TO  T*F  0vER*lL  DOCUMENTaT  t  ON  FILES  REFERENCED  ABOVE,  THIS 
TS$?MTv»Ttnv  TiPf  S  “  C  l'  L  o  INCLUDE  THE  FOLLOWING  FILES. 
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INSTALL. DOC  RHM  CRIM  NOV,  198? 


DOCUMENTS  PROCEDURES  FOR  INSTALLING  ERIN  MULTIpTCMPCRAL  PROCESSING 
DEPTH  A  \  A  i_  SOFTWARE  ON  TH£  POP  11/35  AT  DMA, 

7HE  P»IM  MTP  DEPTH  ANALYSIS  SOFTWARE  CAN  8E  INSTALLED  ON  PMA’S  ll/«5 

R Y  DOING  T«£  FOLLOWING  stepsi 
n  NOV>^£LLO  ?P,2/fRIM 

I  OG  ON  (Vi,?)  f T UIC  USED  POR  THE  ERJM  SOFTWARE) 

?)  ■trS> prp  «,«r*/ne. 

CLEAN  nsc  any  p;i  ES  LEFT  FROM  THE  PREVIOUS  USE 
J1  "C^XO'INT  HT0JDMA 

mount  the  FCJM  DISSEMINATION  tape  (NOTING  VOLUME  LABELfDMA»  800  BPIl  AND 
BLOCVSTIE  HSS512.,  THE  DEFAULT  ON  THE  11/45,  FOR  EASIEST  COMPATIBILITY), 

<11 

COPY  IN  A  COMPLETE  SET  OF  N£W  SOFTWARE,  COMMAND  FILES,  DOCUMENTATION, 

AND  T?ST  DATA  FROM  THE  TAPE, 

5)  "C P >  r  JN'ST  ALL  ,CHO 

USE  Tma  INS7ML  COmmano  FILE  TO  AUTOMATICALLY  COMPILE  AND  TAS'.RUI'.O 
ALL  NEW  ANO  MODIFTFO  PROGRAMS, 
fcl  MCR  >a  YE 

MCR»MF.LLO  S,? 

LOG  ON  THE  PRIVILEGED  MASTER  UIC  TO  COPY  SOME  STARTUP  COMMAND  FILESl 
7,  MP.i»>oTo  /NV  =  tJ0,Pl  STT0.CMD,STTt,STTP,ETTe,ETTl,ETT?,  IMAGE 

COPY  IN  A  NEW  VERSION  OF  THE  DIPS  COMMAND  FILES  to  INSTALL  OR  REMOVE 
ALL  Tm?  ORIGINAL  OTPS,  ERJM  MODIFIED,  AND  NEW  ER*M  TASKS. 

El  mch>*Im.»gE 

finally,  INITIALIZE  dips  WITM  the  NEW  ERIM  TASKS  IN  THE  USUAL  WAY, 


THE  s"i.l."»T*i'5  DIFC;JKF7nN  1k’  ?TMTUA°  T?  Twr  CDMmf»wTS  and 
SCLF-DT  -f.TaTir'N  t  sj  thF  TnGTa.L.C-D  Fit*,  Q.V. 

thf  ,cmo  ml*  install*  all  brim  mult i-t*m»op*l  processing  programs. 

t«!<  S',«-TW4«C  is  4SS'J“*n  tO  RE  put  On  (30, P5,  WHICH  SHOULD  3*  CLEANED  OFF 
f;-'SY.  7hi-.n  log  'n  (  j-i ,  pi  AND  nc  *“CR  >  *  I n?t  all  . 

-A  Pc,,  ----a-v,  cji.tj  vilL  hav*  TC  bf  COPIED  TO  I3,?J  SUB  S  EG'JF  NTi_ v  5  . 

IT  is  sss.'-rn  that  ri.R)  h AS  hc,OL»,  SVGT.INC,  PmCDm.’nc,  anO  t^f  sy37.0PJ 
PcSlDtNT  r  0  "  0  N  ,  ,„ILE  n?R,1(3Dl  HAS  TC  SL I  3 , OLD  (NONE  SUPPLIFD  ON  ThE  fPI*  TAPE), 

SINCE  EVERYTHING  HAD  TO  »fc*  P'JN  ON  THE  ONE  UIC  (30,?)  AT  ERTM,  THE  INSTALL 
"'"“(VI  Ml.'-  E'.I'*-  ,/uJ  AND  «.TK:>  FILER  COULD  not  RE  TMTFD,  it  *S  POSSIBLE 
*"  -•  *  '  '=•  '■"■  ’"I':""  EDITING  ERRORS  CDNV'PTINC-  FOR  I  N  3 7  ALL  »  *  T  0  V  ON  ME  )'u'n/Rr, 
r-.T-;  .  i'.Y  IN  utt nc  T  NO  RJ'.eS  ON  CT“'.’  JIC'S.  I  APPOL0CT7"  ECP  i\“ 

I'  ■•’VVf.rr.-j'.fv  Tm;;  -ay  CAUSE, 

v=  *.DRJ  and  a.LjT  files  are  Included  ON  the  f RTM  TAPF.  B"T  IT  »3  OEM,  TO 
I'iS.'P"  that  Yfv.i  wav=  a  complete  comp  TI3Lfc‘  SYSTEM,  TO  REGENERATE  Them  fro*  Fortran 
S'  -ROE  O'"  iv  the  *.e«p,  ♦  ,*«B  COM“ANO  FILES  AS  15  DONE  IN  T-'E  INSTALL. CM  ML*. 

E  V  r  I  - L  V  -icy  E  V  I  H  T  A  3 Y  G  7  N  ?  7  A  L  L E  0  ATS; 

1)  ema't  --  ♦.catyf-rplo*  AND  H  FORES  SION 

? )  E  apply  --  C'YP’.'Tfc  nijltI-temporai.L''  ADJUSTED  N'w  DEPTH  I«AGE 

(ALSO  USEFUL  TO  SET  AN  AVERAGE  OF  UP  TO  5  DEMw  IMAGES) 
SI  *ERROR  --  CALCULATE  ERROR  PROPAGATION 
m  0  '  T  F  7  f  D  dips  R  0  y  T  t  •  "  1  a  •?  F  t 
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V;  0 E  "■  T m  --  *RT"-7'-'M.  DIM  S  7  •, ? L C  •  T 1 1‘ E  D'  Mu  4 _ s " o  ; t i'«  ,  *:th 
,TD  (v.,T:N  OEP'i'-1  pj 7 NT)  and  i -DEEP  'P"p"h  Xw*-3“ 
GENERATION')  “ODIFRh  TO  PROVIDE  OUTFITS  F*R 
*■11“  R  R  C  G  E  A  “  S 

SOwF  DIRS  ROUTINES  ARE  RFI.ISFP  msrf  f  <M  CONVENIENCE  VI.  NEW  EM"  IN'TFPA*  ACES  I 

6)  F  0  n  *d  T  H  N*4  Tv  T  F  R  r  c  E  TC  CALL  EYIFMNG  F  v  ?  C  *■'**'.  S  P'.’tNCS 

7'  Fif-i.''  --  7  1  NT  C  •*.••'-  TO  CALL  E  7  I  S  T  T  V  G  '7 '•  L  Y '■  C  N  D*  *  IMTIDN 

iv  f  T  '■'  a  L  v  or*;'7  v.‘L  “'tM  Tasrj  ;-F  r."i  ‘fd  nkT  JT  ■'  i  p  7  I  'n  ; 

A)  o  |  •  ,*  t  A  LOAD  TNG  :••#■*•.-,  '  7  L  '  D  v  T  ;  l 
S’,  oc  ITT  »  --  C  7  C  L  !  *  n  I-A'-L  PLANES  F*!  0I*'bLAy 
1M  SA-CTA  --  I  RAGE  CALCULATIONS  (E.S,  ADDING  TWO  7  -  A  G  E  M 
It)  HISTTA  --  HISTOGRAM  IMAGES, 


SURROijT IN*  S  Kpp  AND  I M  0  c  E  P  7 v  THE  ORIGINAL  DIPS  DEPTH  TASK  MUST  BE 
M  CD  I  M  r  D ,  CO-MAN?  FILES  KOP,r!iP,  I'DfcEe,F»P  AND  DEPTH, TK8  (ONLY;)  A’E 
INCLUDED  (K[1P.F«P  AND  I“DEEP.*a»  INVOKE  DEPTH.TK5  FrR  DEBUGGING  K0P  AND  IXQEEP) , 

•  LSD  FDR  CDMPLFTFNFSS  Twr  ORIGINAL  DIPS  DEPTH, FTN  (AS  " OE P THOL D  ,  FTN " )  AND 

MATCHING  OEPTm.CBJ  AND  OFPTH.lGT  AR?  INCLUDED  ON  THE  E R I M  TAPE,  HOWEVER  IT 

WOULD  RE  SAFEST  TO  USE  T«E  ORIGINAL  DIPS  DEPTH  FROM  (S,?J  IN  CASE  jt  HAS  ANY  REVISION 

FINAILV,  you  mAY  WANT  TO  PUT  K CP, OBJ,  IhDEF®,OBJ,  AND  AREA, OBJ  IN  THE  MASTER 
LIBRARY M,P] hc.OLB  CONSISTENT  WITH  PRIOR  DIPS  USE  (AND  SUBSEQUENTLY  USE 
(5,?1  0EPTH,F/1P  ANO  DEPTH.TKB  NORMALLY), 

I*  YOU  WISH  TO  PUT  all  EMM  punoRAw$  UNDER  (J,i)  ALONG  WITH  THE  ORIGINAL 
DIPS  PODGRAMS,  THERE  PR0R4BLY  WILL  BE  NO  PROBLEMS,  JuST  CHECK  FOR  EXPLICIT 
REFERENCES  TO  Thf  i ! I C  1 30 , ?1  ,  I  INTENTIONALLY  LEFT  THEM  IN  FOR  TrF  STTJ,l,3 
AND  tuA  r.F  COMMAND  F  ILFS  AND  MOOIFIEO  1 3  Z  ,  e  J  OE  p  t  H  .  TK  B  FILES,  AND  may  have 
OVERLCOkEO  30mE  ELSEWHERE. 


f  [5*,?)  TNST&l'-.CmD  cBJM  ‘TV.  !°S? 

)  T  H 1 5  >!LE  INSTALLS  Ml  FPIM  “UL7 I -TEMPORAL  PROCESSING  PROGRAMS. 

I  Th:s  tTFTiil.)'  ]S  ASSur'i-0  TCI  «6  PUT  ON  r$',,2),  W“IC“  SHOULD  Sc  ClFANEO  OFF 
I  i!»jT,  Then  LOG  ON  A  NO  00  mcR>PINST  ALL , 

f  (A  PEL  CO““ANO  PILES  will  have  to  BE  COPIED  TO  [5,2]  SuBSEOUENTLY)  . 

; 

1  TT  ’  F  *  SS;jvF  p  that  [T,21  “AS  “p.013,  S v  ST , INC ,  P^COM.INC,  ANO  THE  $YST , OB J 
;  ■■fO;orvT  pO““0\,  whue  UPS,  12?:  “*3  TC3LlS.0LR  (NOE  SU°PLIc.,  »v  THr  EGIM  TA“E)  , 

1 

r  fT’.op  *vp®y *“1 NG  “AO  TO  BE  PIN  on  thj  ONE  L'lC  [JO,?]  AT  ERIM,  THIS  INSTALL 
j  file  plus  *,f«p  ano  «.T<B  r  I  l  E  S  COULD  NOT  BE  TESTED,  IT  IS  POSSIBLE 

I  T“Epp  asp  sqmc  MINOo  poitInG  ERRORS  CONVERTING  FOP  INSTALLATION  ON  THE  D“A  11/45, 

I  *  A  P  T  T  C"L  *  P'.  Y  IN  PEPt’ENCING  PILES  CN  OTHER  UIC'S,  I  APPOLOGIZE  FOP  ANY 
I  TNCONVENIfnCES  ThTS  “AY  CAUSE, 

I  <.,'■!  t  ant  ».i.st  r ;  l  p  S  ape  InDLUOPO  on  the  ERJM  TAPE,  But  it  T$  BEST,  TO  INSURE 
t  ‘hat  v ”ij  H.W*  A  C C “ r L P T r-  CT“® ATIR'.E  SYSTEM,  TO  REGENERATE  T“*M  FROM  FORTRAN 

i  r>o ,'»rp  coi-p  ano  t*<e  .,p«p  ano  *,t*r  command  pilps  as  is  done  in  this  pile, 

* 

I  *  me  following  WAS  presumably  ALREADY  OPEN  DONE  OR  YOU  WOULDN'T  BE  USING  Twis  FILES 
C“cLLO  V,P/F“IM 
J  “  ' 1 N  T  M  T  T  t  0  “  A 

t  't-'i  f.,  »] 

;  *  I  »J  S  7  A  l,  I  ,  CMO  , 

f 

t  CO'IPTL  p  and  TASKBUILO 
I 

S“UJ.FcP 

*  E  .V 1 T  H  ,  P  4  P 

*  ■  •••r.Y.p'J® 

P  •;  2  t  ■>  C  4  ,  A  R  E  »  /  •'  S  P  «  A  R  p  A 

«•  ':jtt,p«= 

*-  A  »  P  V  ,  P  «  P 
*EP“R'’otPap 

,  SUNPTUT INES  KOP  AND  TMOEEP  IN  THE  ORIGINAL  DIPS  DEPTH  TASS  MUST  A<.5D  BE 
,  “OUPTPO,  C 0 “ “ A n o  PILES  KOP.PaP,  IM0SEP,F«p  *N0  DEPTH, TKB  [Onlyj]  APE 
,  :  .,-li  0  *  f  “  0°  ,  p  ap  ANO  I M  n  t  e  ,  P  a  P  INVOKE  DEPTH. T<q  pop  DEBUGGING  <DP  ANn  ImoeER), 

.  a  ■■  r  n  (  “R  c-!'!5lft<tnfss  T“f  ORIGINAL  0 1°S  OEPTh.PTN  [AS  "DEPTHCLO.FTN"!  AND 

!  '*  A  T  o  - 1 N  G  “».3Tw,-pJ  A  0  DEPTH,  lST  ASfc  INCLUDED  ON  THE  E  R  I  “  TAP*.  wOw'FVER  IT 

I  o 1 1 L  0  °E  f  a  f  F  s  T  TO  USE  T“E  ORIGINAL  DIPS  DEPTH  PROM  [3,2)  t 

ra»  DFPTH.ot-PTw/.SP*  fS,?J  DEPTH 

*4®  *DP , V D° /- FP* K 0® 

c,:r>  I“r  =  EP,  Im0FPP/hSP*TM0EEP 

«nsPT“,’sa  ••  FINALLY  TASSBUILO  THE  ORIGINAL  OTPS  OEPTH  WITH  NEW  ERIM  sop,  IMOEEP 

r 

I  NOVI,  I  OS SING  ON  fS.P],  DO  I 

;p;p  iS,P]  /\ivt rsH(?i 5TTP.C-D.3TTI , ST T2,ETTT,ETTi ,ETT2,  IMAGE 
)  and  1,'SF  a S  PREVIOUSLY, 

I 

i  FINALLY,  *th  M4r  WANT  TO  OUT  FOB, OBJ,  IMOEEP, OBJ,  ANO  AREA. OBJ  IN  The  “ASTER 
,  i (  ?*»,»*  ?s,?i  “o.olp  CONSIST*  T  with  PRIOR  DIPS  USE  (AND  SUBSED'JEntly  uJF 
I  r  5,  2]  r,p?*H,P  .’!P  ,  oe?T“,Tk®  NORMALLY), 


ALL  PROGRAMS  WRITTEN  and  MODIFIED  BY  ERIM 

1*  OTPS  MAIN  MENU  ROUTINE  MODIFIED 

i*  OTPS  SMOOTHING  MENU  REUSED 

1»  OTPS  POLYGON  DEFINITION  ROUTINES  REUSED 

1*  T“! ?  rCDL“ ROOFED  AREA  IS  NEEDED  9Y  ESCATT.TKB 

j«  N'E-  L p  1  “  SCAT-fcRPLCT  AND  REGRESSION  ROUTINE 

1,  new  t P I m  ROUTINE  TO  APPLY  ADJUSTMENT  to  [AVFRAGEO)  depth  IMAGES 
j  *  NEW  PRIM  ROUTINE  TO  CALCULATE  ERROR  PROPAGATION 
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Terim 


OVERVIEW. DOC  <hh  tRt**  NOV,  IRflP 

SEE  Of admaist.oOC  9f-F0R*  USING  THIS  DOCUMENTATION  PILE, 

TH  T '  ->Vr  =  VTEH,r»oc  DOCUMENTATION  FILE  IS  INTENDED  TO  PROVIDE  A  PRRAO 
avjc.  0 e  T-r  psTm  M T P  MIJI.TI-TEMPORA..  PROCESSING  nrpTf,  ANALYSIS  PROGRAMS, 

"'"Tails  C^nSUi  t  pFiT'FMST.noc  AS  a  DIRECTORY  to  THF  remaining  DOCUMENTATION, 
and  RUNNING, TOC  P OR  Ov"h  all  INSTRUCTIONS  ON  HOW  to  PROCEED  USING  the  Er:m  SOFTWARE, 


ENTIRELY  new  ERIM  TASKS  INSTALLED  ARE! 

« 1  FSCATT  --  SCATTPRPLOT  ANC  REGRESSION 

P)  HARPI.Y  --  C P m c i j T t  MULTI-TEMPORALLY  ADJUSTED  NEW  DEPTH  IMAGE 

(ALIO  USEFUL  TO  GET  AN  AVERAGE  OF  'J P  TO  5  DEPTH  IMAGES} 
1)  gtpRCR  --  CALCULATE  erro*  propagation 
MODIFIED  DIPS  ROUTTnf.S  AR£i 

<i)  MUJ  --  main  menu  PROGRAM,  MOOIEIEO  TO  ADD  ERIM  memj 
S)  DFPTH  ..  ORIGINAL  D I °S  SINGLE-TIME  DEPTH  ALGORITHM,  WITH 

*•"»»  ( f V 0 w N  depth  bGI'T}  and  IMDEER  (depth  image 
GENERATION)  MOOIEIEO  TO  PROVIDE  OUTPUT'S  FO» 

E R I “  PROGRAMS 

SOME  OTPS  ROUTINES  ARE  REUSED  hfRE  FOP  CONVENIENCE  VIA  NEW  ERJM  INTERFACES! 

M  ES“CTh  ..  new  Interface  TO  CALL  EXISTING  SMOOTHING  ROUTINES 
7)  f  POLY  --  NF/  INTERFACE  TO  CALL  ev  *  ST  TNG  P  J  L  v  D  )  N  R  r  :  N  I  T  I D  N' 

SNT  FINALLY  SPMF  ORIGINAL  *“  I  ®g  T  A  G  <  T  COULD  DE  REUSED  WITHOUT  adapt:-".': 

*)  "'I'TTA  -  -  LO'Or-G  I"- GLS  D'vTD  CDmtjl 
R)  OCI-TA  - -  CYCLING  image  plane?  FCR  DISPLAY 
)?)  RAYTTA  --  IMAGE  CALCULATIONS  CE.G,  ADDING  TWO  IMAGES) 

in  histta  --  HlGTCGPtM  IMAGES, 


VERIM 


th«  FIRST  NEW  9  0:JT  I  ME  JS  "muJ,"  Th?  INJTIALLV  RUM  ROUTIN'  WHICH  0IS»LAY3 
The  “l'v  MEN'S.  it  rs  IDENTICAL  TO  Thf  ORIGINAL  OIPS  "MUJ"  WJTM  The  ADDITION  OF 
ANDTHFR  Si;bwc\ii  FOP  TME  FR  I  M  '"’TP"  Ml'L  T  I  -  T  EMPQR  A  L  PROCESSING  DEPTH  ANALYSIS 
fVlwAFF  IP,  JUST  after  the  ORIGINAL  DIPS  D'PTh  ANALYSIS  IN  The  MAIN  M F N U )  • 

RE'DRE  YOU  APE  ready  to  DO  the  e R i m  MULTI-TEMPORAL  DEPTH  ANALYSIS,  IT  13 
NECESSARY  TO  DO  ALL  T'UGF  WARFI'iG  AND  REGISTRATION,  TRUE  DEPTH  CALCULATIONS, 
DEFINITION  nr  POLYGONS,  A  *  0  ANY  OTH( 3  DESIRED  ANALYSIS  THROUGH  THF  ORIGINAL 
OI^S  SINGLF-TIHF  oeRTh  ANALYSIS  'OR  EACH  OF  THF  IMAGE  OATES  TO  P-E  USEO  AS 
INPUT  for  The  n*!Ji  TT-TfmRORAL  DEPTH  ANALVSIS,  THE  new  FRIM  SOFTWARE  IS  NOT 
INVOLVED  HgRt,  IJ'EOT  for  Two  MOO  I E I C A T I  ON 5  TO  SUBROUTINES  UNDER  The  DIPS  DEPTH 
° R ^ r R A h  TO  ALL''"  PUTTING  OUT  INFORMATION  NEEDED  FDR  T«E  ERTW  PROGRAMS: 

II  fOPTipv'ALlYS,  IN  HDD,  THE  K\dhn  DEPTH  ROUTINE,  P.JT  CUT  T«E  TRUE 

(,FP»W  F, DINTS  IN  A  PSE'JDO-I-'A.r.F  -Hi,’  CAN  PE  INPUT  TO  THE  £R  JM  FSCATT 
oCATTFROLDT  and  REGRESSION  ROUTINE,  IE  COMPARISONS  TO  TRl)6  DEPTH 
W I  L 1.  BE.  WANTED  LATER.  AMO 

PI  (MANDATORY  FOR  MTP  PROCESSING! ,  IN  T  MOEE  p ,  THf  ROUTINE  WHICH  NORMALLY 
r.FNERATFS  A  SINGLE-TIME  depth  DISPLAY  ON  THE  COHTAL.  *LS0  RUT  OUT  AN 
I " A o F  file  ODMTAINIMG  The  TRUE  SINGLE-TIME  DEPTHS  NEEDED  AS  INPUT  FOR 
T‘»p  “  U  L  T  I  -  T  F  M  P  0  R  A  I.  CALCULATIONS,  TO  EAR  GREATER  PRECISION  Than  T”F 
Hr  (,L  *T$?L»Y  fNDR-ALlY  ,1  meters,  PUT  CAM  SE  SPECIFIED  CY  USER)  , 

TwTS  OUTPUT  FILE  CAN  PE  SHIPPED  IF  NO  MTP  PROCESSING  IS  INTENDED, 

THE  T-<D  most  important  NEW  ROUTINES  FOR  THE  ERIM  MULT  I -TEMPOR A L  DEPTH 
PROCESSING  A’E 

n  Fsru*'.  WHICH  PRINTS  SCATTERPLOTS  and  CALCULATES  REGRESSIONS  between 
DIFFE-i-NT  SINGLE-TIMF  depth  images,  OR  AVERAO.F,  OR  REFERENCE  IMAGES, 
Dw  true  0  F  p  T  w  ;  NORMALLY  FOR  POLYGONS  OF  SELECTED  BOTTOM  REFLFCTANCE 
types,  OS  OPTIONALLY  bftwfEN  WHOLE  IMAGES  AS  A  FINAL  CHECK  OF  tm& 
op  ?UL  TS, 

?!  P.  apply,  which  APPLIES  The  REGRESSION  COEFFICIENTS,  SELECTED  AS  A 
REB'itr  of  uStNG  fscatt,  TO  SO«F  COMBINATION  of  t«E  INDIVIDUAL 
5  T  N  D  i_  F  » T  T  w.  E  CFBTH  IMAGES  TO  GET  AN  IMPROVED  DFPTh  IM4GF.  any  ARBITARY 
wftouti-d  AV’FSA0,E  OF  SINGLE-TIME  IMAGES  Can  BE  USED  AS  Tw£  Input, 
hardly  AI.SD  provides  a  Convenient  way  TO  OBTAIN  a  buff  AVERAC-ED  I«agE 
(4  straight  average  of  several  SINGLE-TIME  CEcTH  IMAGES!  TO  L’SF  AS 
ONE  OF  T»H  POSSIBLE  REFERENCES  f;p  ESCaTT, 

FERVOR  IS  ANOTHER  NJW  ROUTINE  WHICH  MERELY  CALCULATES  THE  RESULTANT 
RR-Rf. GATED  fRWPRK  WHICH  RESULT  FROM  APPLYING  Tw»  REGRESSION  TRANSFORMATIONS 
TG  USER-SPECIFIED  FPPOPS  IN  THE  ASSU"E0  PARAMETERS  ( K  AND  RmQ)  , 


BO«F  OI.d-f  EXISTING  DIRS  PROGRAMS  ARE  ALSO  SUPPLIED  IN  Thf  eriM  M T 9 
m'nu  for  C ?NV’'  NJFNCE  JUST  RFC  »  U  5  E  they  ape  LIKFLY  TJ  SE  FREUUENTL Y  USED  IN 
r3N JUNCTION  WITH  the  MULTI-TEMPORAL  STAGE  OF  THE  PROCESSING,  TWC  ARE  CALL'D  By 
NF. W  fhtm  interface  MENU  PROGRAMS,  THOUGH  BASICALLY  UNMODIFIED  DIPS  ROUTINES: 

1)  ESMOTH  --  SMOOTHING  «dut I mc  s  TO  SMOOTH  ANY  (51EX51?!  DEPTH  IMAGE 

P  f  F  DP  r  USF  IN  ESCATT  AND  FA?PLv  ( » A RT I  Cm.  ADL Y  RECCMm'NOEO 
if  USING  F  5  r  4  - t  TO  COmraRE  »  OFcrw  image  TO  TRUE  DEPTH 
POINTS! 

— !  p  =  n|  y  ...  T  -  -yfctns  ADDITIONAL  POLYGONS  OF  KNOWN  OC’TCM  REFLECTANCE 
TYPES, 

AND  FOUR  EXISTING  CIPS  ROUTINES  COULD  8  E  C»LLE0  A;  IS  wITHOijT  A  N"w  JnTERFACF: 

X)  •' T u t T a  --  loading  image?  into  thf  twrfp  c 0 m t a L  Image  Planes  for  OISRLAY 

i'  ’ r  ’  ’  T  A  —  S'L'CT  TNG  W-JOh  IMAGE  •>LAN£  '0  DTS  =  lam 

f;  ■<  i  y  t  T  a  -  -  pfw-  -•  **  m  c  \  g  —I  'CELL  A*  EC  jS  0  f  c  f  a  *  I C  n  S  'add:*  2,  SCA-.C*  g,  ’*,  t 

u>.  : m.'.g-s  (hut  also  f*"»ly  to  got  a  yeigh-  v  a .■  r  a  ; e  of 


•G  c*  VASES. 
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VERIM 


r,-EC’<'"'JT , o?r  khw  erim  nov,  i98? 

SE*  bfahmfj ST.DOC  aEc0?F  USING  THIS  DOCUMENTATION  PILE, 

■’■-'is  checkout. hoc  documentation  file  describes  song  sa-ple  data  ano 

1  ,ITP'JT  FRo-i  Thf  F  R  T  M  OSQORA-S  WHICH  CAN  3E  COMPARED  WITH  RESULT?  C8TAINEC  AT 
-)M4  rj  v  f  9  t  r  V  yoijP  INSTALLATION  ANO  RUNNING  PRCCEOURFS,  INCLUDE®  ON  ThJS  TAPE 

iff  r,cuE  input  images,  polygons,  ano  sample  outputs  from  the  er;m  prccra-s 

■  - ■ C H  vO"  CAN  COMPARE  WITH  YOUR  RESULT*, 

f  files  omA'.T  “A  Through  dm  a  R . I mG  ABE  SInGLE»"ImE  0E®Th  I-aGES, 

"fl  T 1  EAC-1  other,  focm  T  h  £  o  j  m  j  m  i  AREA, 

FILE  D“»AV1 ?S4, I“G  IS  AN  UNWEIGHTED  STRAIGHT  AVERAGE  CF  T-ESE  F0uR 
SINOlp-TThf  IMAGES.  IT  *A$  OBTAINED  8Y  RUNNING  ThF  EABPLY  PSOGRAu  f  V  I  A  ThF  F9IM 
MTp  MENU  S  F  |_  F  r  T  I  o  N  FOR  A  R 3 1_  V  }  ,  IN  P'JNMNG  EABPLY,  a  INPjT  IMAGES  .'ERF  REQUESTED. 

*  •Tr.M’S  *FRF.  mffa  jl»f  j  (-ENCE  EQUALLY  I  ,  0 , 1  ,  0 , 1 . 7 , 1  ,  Z  ,  NORMALIZED  TO  ,?5  Each  TO 
•=.  «  TO  I  ,  m  C  V  f  P  A '.I  *  AmP  /NITv  TRANSFORMATION  SELECTED  B*  DEFA..T  f  M L  T  I  P  L  I  E  R 
AOOI'IVF  COnc’ANT  ?,?)  ,  0°V:CU?Ly  T“E  INPUT  FILE  NA“CS  a'FRF  E  N  T  r  b  F  d  AS 
'-•:!,!UG,  OMAP,I“G,  ETC,,  ANO  T  «E  OUTPUT  FILE  N  A  “  F  1$  D  N  A  A  V  1  2  3  «  .  I M  G  , 

T«F  f  1 1  M!-PCiR.°OL  IS  A  DEFINITION  of  POLYGONS  Detained  FOR  THESE 
images  at  Om4  p,.,oiNG  The  Ofmo  on  OCTOBER  1  ,|R5t. 

ThF  CO,  5  F I L  f  S  E SC  A T T , 0  A T M  .  ,  ,  fl  A»E  'H*  PRINTOUT?  RESULTING  From 

R.-nT'G  FSCA'T  F“!»  THF  A  1  3  A  G  E  IMAGE  U  *  A  A  V  1  J  3 U  ,  I v  G  AS  A  REFF  =  ENCE  ON  T-E  Y. 

'.*!$»  vfrs’,'5  Each  of  tm£  F2  JR  Individual  In°JT  IwaGES  PM  a  i  ,  I  mg  , ,  ,  Dm  a  a  .  I  mo,  Jn 
CN  'mf  x.axTS,  nof-ally  T-IS  appears  cn  the  LINE  PRINTER,  PUT  IT  LAS 
5  *  *  *  0 u  f  r  rn  A  OISF  F  HF  ay  REASSIGNING  UNIT  3  fCR  SEE  ESCATT.TK0  ASG*LBI3I. 

’0  Avon  INCLUDING  SAT,'°»TEO  VALUES  FRC“  DEEP  WATER  or  AN0MAIIFS,  TmE 
•N=  *  'ita  SO.V.E  *•«  I  TMT’FO  to  P-R50  ON  The  j.tnjs  (85«  -CLIO  have  wCRKEO)  , 

*•■•  »)•,  »:•,>.  -G  arc  RR-N’FO  In  the  captions  'JNPE5  T»f  SCA-Tf  RPlO’S. 
no*:  *-*?*  IN5  *  Image?  -*E=E  OFFTH.PBrCEGSEO  a T  frjr,  wI'h  -.0  r'ROSLE“D  at  the 
L  0  > c  c  LNO  OF  ’HE  jaT*  scalp,  hq*cVE9(  pms  image*  GEvERA’ED  «V  ThF  £R I m» m do  I®  IED 
IMo-.po  BOUTIN*  UNDER  TMC  OTPS  “E B  T  H  PROGRAM  4T  n"l,  IT  wO'JlO  RE  NECESSARY  TO 

-  fLAVO  ooos)  -NO  !  UNDERFLOW  CODE  5  In  normal  SING  F.TImF  0ErT“  tmag»S 

r '  •  ,  c -  V  T  e 1 .  r  -earns  t!"v  r  „  c  fP7"„wh.jtfTE0  SO3  ROUTINE  ‘.ora  T-F  "tet  0>  F  ?  '  H 

1  .  0":  •  C  ;lD  AL'A’S  have  to  3E  SKIPPED  AS  A  SC-OAT-.  C:"F). 

AREAS  paThfr  Thau  t-e  fstjoe  Image  FIlE  REa  WERE  SELCCTE', 

'NO  n  ( E  -V_  gIv'N  as  T*'F  POLYGON  file,  the  fc'mst  g  y  :■  S  T 1 0  N  wmEThff  al_ 

=  ".m-o  ?  ,c-'E  *  a  T  f  o  a  a  ?  answered  NO,  Tm£s  as  the  DESIRE'1  POLVGONS,  Af.  ?  0  *  T  0  M 

Tycc  a  r  f  a  OO-'FS  f  p  r  *  n  T  f  o  gv  ESCATT)  WERE  INDIVIDUALLY  fN'EREO  EXCEPT  "7"  FIR  CEE3 

‘  *  *  f  *  ,  Thwst  h  .  %•  T.  =  •.  ?Sf\  i.  IS  Tf->  I N  Thf  CAPTION?  AT  Thc  tf’y'y  PF  Twr  «-aT’FP?lOtS, 

i  ■•*,  ante  codes  instead  0s  CCJnT  *  *m»olS  -'E  PICKED  f."’  *  -I  clO*. 


t  H  f  -  E  J  ,  s  o  0  u  v  A  T  ID  N  pfcjIRfm  Foa  IN  0  A  R  P  i_  v  TO  vri1'  EiC1-  Oc  Those 

r  — ,  i '  r.  s • • o t c . t  •  m  t  , . c c r h  :“£?;?  ovai  ,  i  m  g  ,  ,,dmal,img  tn  t“E  * -  a  t i 3 ,  to  t«f 

A  A  G  F  ’  i-a'.e  o„A  W1  ??u,  [MJ  l  ?  E  o  AS  A  I'FFCJfMt  oN  t„E  y-AXIS,  IS  GIVEN  AT  ’h« 
a-T»-u  HE  -„E  crjtTtRP'.rT  CAPTIONS  hY  T-E  REGRESSION  '.INF  Y  »  A«x*«.  (AlThO.Gh 
,  o-  t.  Twr  ff  ?A“OLF  SCiTTFROlOTS  R'.’V  IN  farlIER  VERSION  0?  Thf  RS'GR  *" , 

’_C  FI*'.'  FSOjYT  ALSO  GIVES  *UE  I N  V  f  a  5  E  REGRESSION  LlNe  *  •  A  *  *Y*P  * ,  IN  CASE  THE 

.  •  c  c  ?  ,  •  -  r.  I  '  t  G  c  TS  0  !  T  CN  yhF  y  -  A  X  I  *  I  ,  T  w  I  S  RCG3ESSI0N  FIT  IS  OF  COURSE  OEsENOENT 
*  .  *-*  o-oioe,  O.alI'y,  ANO  E?PFF«-\Y«Ttv*NS.SS  CF  Th?  bolyGO’  S  ..  IT  I »,  us  *0  Thf 
:  ••*}  T  TO  "SKI  s  ")r  yi  s  :PHlIEG  P  ;l*3“N5  CD/E»ING  A  GOODLY  PANIC  ~r 

i"U“  T  y  0  a  s  :  *  0  rrp'MS,  ano  T-aT  hf  VERIFIES  THE  R  E  A  SON  A  ?L  E 1  E  ?  S  Cr  T«f  '  R  A  N  SF  0R“  A  T  1 0*0 
Th'.<,e  R  c  *  U  L  *  S  *  H  0 »  n  in  The  SCATTERPlOTS  ARE  TYPICAL  Fa]Rly  RE*SCNA»lE  PESLL’S, 

finaily,  GIvfn  a  satisfactory  transformation  as  DSTitnEo  as^ve  frc.m  the 

ESC*tT  0  'NS,  THF  EABPLY  PROGRAM  i  g  bi,m  AS  A  *  o  y  f  ,  '.SING  EITHER  ONE  SINGlE-TI“E 
-  r  a  -  w  I  WAGE  ,  -9  AN  ( 0  ~ ’ T  0  G  A  L  L  *  W  F  I G .  T  F  D I  AVERAGE  OF  SEVERAL  AS  T  u  F  INF  *  ,  A  N  0 

T.r  y  r  A  «  x  ♦**  T'.|'.SFM*ITI"*'  aEL'CYEO  FwOw  T-F  *'  f  5  T  SCA'TERP'.ots,  *0  ?E*.f-A*E  '“E 
:«!•»■■  C'.*3J*  0  s  3  T  “  :  M  A  G  F  t  SINCE  A  SATISFACTCRy  cwCTCE  OF  -RAVSFORMATICN  is 
•0  S-M?  FIT^’  DFmfnOF'jt  ON  T«F  OFFFATOR’S  SFlFCTION,  NO  PAR'IC'JLAR  5a“PlE  IS 
I  NC  L  U  of**  hf  k 
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punning, doc  '3MW  C!?rM  nov. 


SEE  PBAOMFj*.T,DOC  BEFORE  USING  T-ls  DOCUME  NT  AT  I ON  FILE, 


■  a  »  4 


"  u  ?  c  3  . 
l'i.  >!'• 


■inn.ooc  opcu“Entati"n 

$  m r'  :J 1  D  p  q  D  C  f  F  D  TO  U5e 
IS  4  N  INTEGRAL 
■E5  FACILITY. 


IS  INTENDED  TP  DESCRIBE  T M  OVERALL  DETAIL 

Tub  £rtw  -l ^ tp "  “ULTI-TE“  =  C9Al  PROCESSING 

Par*  o f  0  u  a  •  s  0 1 3  S  *iTco  “Fatm  PROCESS  I  v  3 
0'!i*S  p"  =  '.  1  / a 5  c o“=l-TER  pt'iLITy,  this  reDjTRFj  as  integrate"  uSE  Cc  t-f 
0=TGIM_  nru-.  pr"CFS3TNG  TwRorOH  COMPLETE  S I  S'  GL =  *  T I  ME  DEBTH  IMAGES,  FOLLOWED  ST 
icr  If  THc  VBA  F5;»  sniTAAPE,  when  PININ';  thf  3EC0N0  EPIM  part,  IT  mat  8E 
wpi.^piji.  3CR  TJF  :\AP?t,  TO  CONSULT  Tuf  SAMPLES  DESCRIBED  IN  T-E  "  CmEC  *  OUT ,  0"C " 

Of  STATION,  ASO  FOR  “"RE  OETAJl  TH4N  CAN  BE  PROVIDED  -£RE  THE  INDIVIDUAL 
R  P  “  R  A  “  S  S-CUL"  0  S  CONSULTED, 

tub  f;o5T  cnr.FS  QE  T-E  p="CESSING  REGJlBE  USE  OF  the  ORIGINAL  01  =  5 
S"Fyw4»-  jn  T-'E  USUAL  way. 


TT  TS  ASS'JMEO  THAT  there  ARE  GODO  pROCESSABlE  IMAGES  COLLECTED  FOR  TnO 
" R  OAT£S  OVer  The  jamb  AREA  FOR  Thf  ERIM  MUL T I • T E MP OR  A L  PROCESSING  TO  SE 

CF  ;SE,  it  13  manOATPRv  that  Thf  3F  images  BE  SPA’IALLv  =EGISTEREO  to  EACH  "ThER. 
.'-?,.E  I "  “at  \ "  R  m  »  l  l  *  PE  PREFERRED  TP  0"  T“IS  *T  REGISTERING  AND  »A  =  3INr,  each  TO 
T  h  ■'  SJMf  STAA'.AR"  “AP  COORDINATE  SYSTEM,  THJS  TS  NOT  NECESSARY  Fpe  T-r  ERI“ 

»  L  G  0  »  I  ’  h  m  O  r-  R  s  E  . 


IFTF?  PEG!  STPation  AND  SELECTION  pE  m^ychIng  area  5 1  ?*S  1  ?  Si  I1*  IMAGES, 

Tuc  v1  C  R  -  4  L  nTWS  s  :  NRI  E «r  !  MB  H  E  3  T  “  °R"CESSING  S“0'JlO  BE  COmp-.ETEO  THROUGH  OBTAINING 

T--  -'rpth  -M4-.ps,  t a* o  or  t^ee  o:rf£ufnhc3  should  ® e  no’eo  “E’e  if  t«»s  s"aoe  is 

TP  aer  FCI.LOwF)  ay  TUB  frjm  “:JL  T I  -  Tp“RpR  AL  PROCESSING, 

1!  IF  IT  m:-,ht  BE  PE  S I R  £  0  TO  CO“PA3E  any  Cp  the  DE°Th  I“ASES,  aEEO»E  OR 
aFTFR  APOUtONAL  M  ilti. temporal  RRPCESSING,  TO  TRj£  DEPTHS  *n  th£ 

F.SI-  ESCaT'  SCATTER® lot  program,  than  DURING  Thf  TRUE  DEPTH  PROCESSING 
T“F  c  f  e  R  a  t  p  p  should  R  E "Uc  S  T  h  R I t i n  G  THE  TRUE  OEBTHS  IN  * - E  NEW  EBI“ 
r  =  Sr  I'OI -T“A'.F  FILE  forma-. 

»j  ,Y4\  SELECTING  o' l  TP-"'  A  -  A  S ,  Thf  "PfRATC3  SHOULD  <  E  E  p  in  “IN"  that 
h f  4 T I.  L  BE  WANTING  0 ‘  3 R E S E n t  a '  I  ■; E  « 0 T T c m  types  and  DEpT“3  0 v E R  AS  „  I "  E 
A  range  De  "  E  R  T  m  S  A  a  w  e  DA-J  s  E  L  T  A  D  L w  IDENTIFY,  IN  ORDER  TD  GET  GOOD 
S'a«lE,  StGMETOANT  R E SUl  T S  FRO*'  T-E  SjESEOUENT  SC.ATTERPlOTS  AND 
RFGRrSSTON  CALCULATIONS, 

31  AND  E I N  A  L  L  T ,  4";N  CALCULATING  THE  OEPTh  J-AGE,  THF  OPERATOR  MUST  NOT 

CN.Y  display  IT  ON  thf  C  0 M  T  a  l  ,  BUT  ALSO  RSGUEST  T“E  NE*  “ORE  °®EO!SE 
:>rrrw  j  y  <  t  ni-tO'jr  ff  t  P  “jr-jc  A',  ??  A  Z  *  ?  T  -*  ^  Z  ~  5l'*'V‘Tr!^A‘;CEi 

: ' i  n_ t : \ ?-  fi  mpv.s  c::;nt  i  ^ a r t ^ •  j ^  *‘r-T-  :r 

RG.P  ‘'E'fRSI,  IF  5  i  e  a  T  f  f  d«p*hs  as;  REOjI’ED  Thf  CREbaT"*  G-'.'.D 
o  1 0  y  f£3HAPS  ,  R  "FT-PS/CO'JNT  {!▼  mIGht  BE  PR”feRai>lF  t0  .3?  Thf  ji"E 
F"R  all  depth  jmaGES  to  avoid  E'jTJRE  CCNFjSION  InT£b«T)<IlG  T-E“5. 

4  f  t  f  c  TUB  SINr-LE-TI“E  DEPTH  images  APE  OBTAINED,  TWE  "pfrathr  IS  t£A"T 

TD  fnTfp  The  s  R I “  “...'I-Te-pdRal  PROCESSING  “E‘u,  T“Ese  aPE  T»0  » E »  ?'tBS  T0 

Tmf'm.-.TT.»E“P~Ral  ‘  -  0  *  •  s  S I N  P  rr  D°Ta!N  AN  IMPROVED  "FPTh  I“A'.E,  FIRS"  IT  *"  5.N 
TD  do  RPL""S  and  pf  CP'SSION  C  t  l  C  .  L  A  * ! "  n  S  OF  T-E  r'lvi:.*..  31'  *  .  E  •  T  I  “ 

r  f  r.  Y  .•  TM  a-;bs  ve  j*  .  B  S"“«  T,  FI.  EFT:  "  STaNOAPO  C»  »SF£RfNCE,  RESTRICTIN',  *-fc  Q»TA 
♦  0  R  "  l  *  G  D  N  a  P  E  A  S  S  R  L  S  C  *  E  D  r0  RFopcyB\T  a  Vi  R  I  f  T  T  OE  BOTTOM  FcF'.E"1nCE  Tvsfs 
IV1  n;3'Hj,  Thm  DbeR  iTo®  Than  Tya“INES  T“E  SC  A  TTE  R Bl.  3T  S  ,  DECIDES  nr*  TC  REJECT 
ANY  f>  •'  •_  y  0  N  S  'r'iH.d  D  i  t  A  val'JE  a  that  Sffm  TP  a\  unfairly  BIASING  T-E  PES..TS, 

Ft  Fj'.c  r-'«  <"fp  AG  neoESSahy,  once  S4TISFIF-J  HfRB,  The  SECO'. "  S*ED  Is  A 
T  I  "F  mf;han!*g  0  N  ■.  I  N  G  T  ►  F  E  A  P  P  L  v  RR"GRa“  On  CAC-  I ~  I  V  I  0  .  i  l  !<f.T 

y  *  -.G-.v-*  r  -t  *  2  T  4  i y  43?  «~2  avfragf  Or  ••**!  .<’Tr  tub  SrLEO*EO  fEGR-TSIOn 

•r.-nO  Ol"4TT"n  “  a  y  c  They  y*2i  0LOSrLv  “I’l"  REcrRFNCE, 

THE  READ"'  “AY  -ERE  WANT  TO  CONSULT  T-B  "  W  'j  N  N  I  N  G  ,  0  0  C  ”  0  0  C  J  “  c  N  T  A  T  T  0  N  EIlE 
> '  "  fampi.E  p.'S  *T  D!S"USSES, 

4  IV*'I4_  S*:=  p-  "OFF  RJN'INO,  T»i  '-JT'FGS.'"  4’  *  f  r f  Y  i  ;  •  -  -. 

PR“GR4“  ‘"s'*  I".  T"  "  •POSE  A  S/I'AjlE  F-I'FR'NpT  I“;gf,  T-jb  tan  et  vfj* 

INI'  !  *■  ;  4  •_  N'.G.r-T’‘,r  I  v  AG  r  ,  J,  '<OIr'“TED  'p  •„  6”,  -  'f*>  4  V  E  0  A  '  r  "r  f  '  - 

T'-’VID.A.  I  “  *  g  £  f  ,  o '  PJSGIOL*  'R  .  E  '  :  p  '  H  ROIL’S  En:0"ED  IN  An  t  =  1“ 

...  ..c  er.  r  r\  N(  PRIOR  DIB?  <?’  KNOWN  OEprH  POINT  R  o  U ? ;  N  E  cN"FR  The  OIpS  0£®Th. 


thkn  The  e*»tM  “ULTI-TEMPO^AL  PROCESSING  PROGRAM  ISCATT  IS  RUN,  VJA  THE 
SCATTPRPI.OT  CHOICE  ON  THE  EWJM  OTP  MENU.  EACH  INDIVIDUAL  SINGLE-TIME  CEPTh  image 
SHOULD  oe  RUN  VERS"S  TH£  REFERENCE,  USING  ALU  POLYGONS  PELT  TQ  BE  RELIABLE  AND 
Relevant  OVER  A  GOODLY  SELECTION  op  bottom  TYPES  and  depths  (EXCEPT  "deep  water" 

ANO  UNRELIABLE  OR  NON-ijnIFDSM  DEPTH  POLYGONS},  THE  PRINTED  SCATTER°LOTS  3H0UL0  BE 
PYAMINFO  for  DEVIANT  POLYGON  areas  (IT  WILL  8£  DESIRABLE  TO  REQUEST  POLYGON  AREA 
BOTTOM  'YPE  symbols  INSTEAD  OF  COUNTS  ON  THE  SCATTER»LOT  TO  SEE  THIS},  AN'O  POR 
H  T  E  0  DATA  VALUES,  OVFRFLOw  (PSS),  UNDERFLOW  fn,  AND  LAND  (0}  DEPTH  COOES,  IF 
T n,!  SCATTERPLOT  IS  TOO  CLUTTERED  TO  ALLOW  DISTINGUISHING  ALL  THE  INDIVIDUAL 
pr,v^;y.'  *  t  e  *  <s  f  THE  OPERATOR  CAN  MATE  RUNS  USING  A  F  E  VC  POLYGONS  AT  A  TJve,  THFN 
T--  OPPRAIDR  SHOULD  RERUN  SELECTING  OUT  UNDESIRABLE  POLYGONS,  AND  RESTRICTING  Th£ 

Oat  t  L  0 ' c  -T  A'ij  upper  ROUNDS  AS  NECESSARY ,  PARTICULARLY  DELETING  0,  J,  ANO  353  COOES, 
TP  TOC  «»NV  POLYGONS  MAVF  TO  BF  REJECTED,  ANO  THE  RESULTING  REGRESSION  FITS 
UNREALISTIC  OR  unstable,  the  OPERATOR  CAN  ATTEMPT  TO  USE  THE  POLYGON 
MFN)  F N T p Y  UNDER  THE  FRIM  MULTI-TEMPORAL  DEPTH  PROCESSING  MENU  TO  SELECT 
ADDITIONAL  BETTER  POLYGONS, 

ONE  A I  YftRMATTVF  WAV  TO  00  THIS  RECPESSTCN  STEP  INVOLVES  USING  THE  TRUE 
DP-'TUS  rvcnr;)  tN  A  f  P  ?  E  UOO)  - 1“  A  ‘c  AS  THE  REFERENCE,  IN  THIS  CASE,  Thf  SCA7TFRPL0TS 
G->~  'LO  HP  R" N  OVER  THE  WHOLE  0)3X513  IMAGE  ARE*  AND  THP  LOWER  BOUND  SFT  TO  1  OR 
“I'.H'R  ON  Th:  true  DEPTH  IMAGE  (SINCE  »  DENOTES  LACK  OP  ANY  TRUE  DEPTH  DATA]  ,  TWE 
NUMBER  op  P  t  y  f  I,  S  INCLUDED  IN  THE  SCAfTERPLOT  WILL  THEN  BE  LIMITED  TO  THOSE  ENTERrO 
AS  T=UF  f'EPTH  POINTS,  IT  WOULD  HE  ADVISABLE  TO  00  SOME  SMOOTHING  OP  THE  ACTUAL 
CALCULATE  DPPTH  IMAGF  TO  REDUCE  THE  EFFECT  CP  NOISE  IN  THE  CALCULATED  I M AGE  AT 
•h:  fPRR  S U m A ; ' l_ V  VE c Y  SPARSF)  TRUE  DF^TH  POINTS.  THE  SMOOTH  ROUTINES  ARE  INCLUDED 
-ns  fR-m  •  ■  p  n  U  TO  FACILITATE  THIS. 

This  REGRESSION  FIT  IS  dp  COURSE  OfRENDENT  ON  THE  CHOICE,  QUALITY,  AND 
REPRESENTATIVENESS  of  the  polygons  --  IT  IS  up  TO  THE  OPERATOR  TO  make  SURE 
MR  S'JRRi  TES  »ni YGONS  COVERING  A  GOODLY  Range  OF  known  BOTTOM  types  AND  DEPTHS, 

AN"  -HAT  HE  VERIFIES  THE  RE  A  S0NARLcNF  SS  OF  THE  TRANSFORMATION).  THE  RESULTS 
in  the  sample  SCAlTERPLDTS  ARE  TYPICAL  F  A I RLy  REASONABLE  RESUL'S, 

AS  TMf  S'CONO  KEY  STEP,  DIVFN  A  SATISFACTORY  T R A ‘!SF OR « A T I C N  AS  OBTAINED 
"ROM  T nr  f_SC  ATT  RUNS,  The  EAFP'.y  RRCGRa«  IS  RUN,  USING  EITHER  ONE  SINGLE-TIME 
r1T'!  IMAGE,  DR  AN  (OPTIONALLY  WEIGHTED)  AVERAGE  OF  SEVERAL  AS  THF  INPi",  AND 
The  y-i*y*b  transformation  $f.I  ECTEO  FROM  THE  BEST  SC*TTERPlOTS ,  TC  G-nfpaTE  THE 
•m-ctved  OUTPUT  PEaTH  Image. 

FINALLY,  tH£  ESCATT  SCATTsp»lCiT/RFGPESSION  PROGRAM  C*N  BE  R'JN  AGAIN  ON 
T  H  -  RsCL'T-ND  Mgi.  T-.TFMPDRALLY  O-wlV-n  orpm  !«*(•.£  VERSUS  A  REFERENT,  t\"  ?M» 

R’iSU'.T?  CT'RARPO  with  THOSE  FOR  INDIVIDUAL  SINGLE-TIME  DEPTH  I“tr.t;s,  AS  A  WAY 

of  measuring  how  reasonable  the  ofrth  calculations  seem  over  the  entire  I“age, 

TO  DO  THIS,  the  whole  I-AGE  AREt  SHOULD  HE  USED  INSTEAD  OE  THE  SELECTED  POLYGONS, 

B  IT  THV  lD-yR  and  URPFR  POUNDS  can  BE  RESTRICTED  ON  EACH  AXIS  TO  SELECT  meaningful 
S h a l L 0 ‘NO  MODERATE  qfpThS  whf.re  REASONABLE  RESULTS  SHOULD  9E  EY°ECTEd. 


APPENDIX  B 


SOFTWARE  LISTINGS 

This  Appendix  contains  FORTRAN  IV  listinqs  of  the  two  principal 
proqrams  contained  in  the  MTP  software  package.  These  listinqs  include 
EAPPLY.FTN  and  ESCATT.FTN.  A  list  of  all  of  the  MTP  individual  proqrams 
is  contained  in  Appendix  A  (README1ST.DOC).  Many  of  these  proqrams  are 
modified  proqrams  from  the  DIPS  software  package.  Listinqs  of  these 
latter  proqrams  are  available  from  the  delivered  MTP  tape. 


67 


2Tri. 


Z 

ft- 

o 

I 

►Hi 

15 

ft- 

M 

4 

UJ  z 

x 

3  O 

or 

M 

o 

X  ft- 

u. 

ft-  <c 

CO 

►4  X 

z 

3  ft 

or 

4 

O 

UJ 

or 

X  u. 

X 

ft- 

uj  tn 

ft- 

X  z 

o 

cr 

ft-  *4 

4 

ft 

c 

kJ 

CO  ft- 

CD  Z 

z 

UJ 

4 

ft* 

o  cr 

♦ 

_J 

41  4 

T  u: 

i n 

*> 

UJ  z 

•O  JL 

M 

>H  u  o 

if 

4  _J 

c:  h 

_l 

4  If, 

>- 

►  (9 

a:  ft)  • 

4 

z 

►4  UJ  CO 

Z  X 

UJ  •-» 

ft  u> 

4  M 

15  >- 

ft-  15  IS 

rr 

«  -J 

3  UJ  4 

X  tai 

X  ft 

ft  Or  T 

ft- 

M  OL 

Z  ft-* 

ft  cr 

■4 

•4  4 

0.1  UJ 

X 

ft)  X 

o  aJ 

ft-  >- 

IL  H  )- 

UJ 

ft  (D 

O  O  ft 

X  M 

U« 

_J  l*.» 

*4  3- 

« 

r  uj 

u.  u  f 

or 

>- 

13 

15  ft 

15  • 

_J 

ft-  *4 

4  UJ  ^ 

a  x 

a 

3  X 

cr  »-  uj 

O' 

a 

ft  •-* 

ID  ft-  T 

CL  «0 

4 

Z 

>  4  H 

•o 

r> 

UJ 

• 

M  X 

4  CJ  f- 

JO 

% 

t 

»- 

CO  t 

-4  a 

¥> 

i 

ru  ft 

«  UJ 

at 

y- 

ft 

«4Ul 

4  >-  — 1 

o  ro 

_J 

o 

z 

IT.  O 

ft  12 

o  iO 

4 

k. 

o 

* 

■  Z 

r  rr 

z 

M 

AJ  ft- 

r>  »4 

UJ  r« 

4 

ki 

•- 

UJ  Ul'  ft) 

ft- 

X 

a 

i r  o. 

€5  Z  w 

•  • 

X 

•4 

M 

ft- 

4  *-* 

►«  r- 

»- 

z 

ft 

UJ  3 

X  X  J 

ft-  CJ 

0. 

U 

ft  O 

*4  or  4 

—ft  o 

k» 

3T 

CO 

o 

U*  Z 

3 

c. 

tn 

UJ 

X  o 

ft-  »>  ►* 

X  o 

4 

o 

ta.' 

3  U-  12 

UJ 

ft. 

»- 

ft  CO 

6.  C  ft* 

r  «o 

►- 

• 

O  i-4 

►-  ft 

M  H 

X 

o 

1 

> 

3  <0  O 

ft  > 

Uj 

• 

UJ  UJ 

O  ft 

UJ  U' 

X 

-ft 

z  or 

UJ  Z 

ft 

►4 

_J 

o 

3  U* 

>- 

ft 

w 

UJ 

a 

-J  ft- 

UJ 

ft- 

r  tu 

ff*  74 

■5 

&  CO 

•0 

O  ft- 

ft- 

•s 

Q  4 

z 

4  4 

C  UJ 

cr 

*4 

1 

U  -1 

Z  CD 

3T 

ft  “  * 

4 

U 

1 

>- 

• 

CJ 

CO 

o 

O' 

_r 

4  U: 

J 

► 

Cl 

4 

4 

CO 

(D 

.1 

o 

ft- 

u 

UJ  > 

<■» 

a 

UJ 

ft- 

ft  _l 

ft 

0 

ft- 

a 

D 

Ul'  4 

»- 

4 

4 

4 

Z 

X  z 

■v 

Ul 

o 

Uj 

4 

3  4 

cd  o 
3  u 


►<  J  *  * 

_#  *  * 
rr  ♦  * 


o 

fcu  •- 

to  II 

3  Ui 

k. 

«r>  u. 

UJ  tu 


uu  uu  uuuu  u 


IJUt'L'UU  U«J  Ll  U  U  u  C  LI  U 


Z  ft-  k.1 
+-  rr  -*  • 
UJ  c  u- 
um  I  q; 
ft  ?  -« 


r  u  x  t 

o  *>  c 

U  U  J 

•  z  ♦  < 

U  <  U  H 
=  U  t  J 

l) 

OO  JtJ 

u'  z  _• 

«_»  <  4  l»t 

4  X 

u  »-  o  »- 
W'J/ 
ft  ft  H  X 

ft  Z  I-  h 

•—  U'  *-* 

•n  _j  :♦ 
M  u. 

z  i3  o  r 
u>  4  u* 

X  I  ►  li 
UJ  _J  l3 

>-  uj  r» 
<  <»  o-  '• 
ft-  o  u  .•» 
•Ml  >  C 


•  •  •  ft 

X  T  X  1 

O  C  O  i_' 

u  u  n  l 

•  •  «  « 

U  U  U  U  4." 


a 

7  T 
ft-  ‘.iJ 

ft.  a 
u  J 
c  «  w 


t  a  ft-  uj 

D  U 

7  L  7 

O  3 

r  o  ^ 

1  W  ft  v. 

►h  CD  *->  —• 

MX*-  C* 

ftt  —  r>  x" 

x  a  •-* 

z 

4  *  *4  7-  ^ 

—  -*  2  —  <- 

w  *  z 

►-  v  -i 


r->  *  O 

ft-  7  r~  U 

Ai  3  O  3 

—  ft  t  ft-  a* 

mir  2  z  r 

ft-*  3  C  « 


ft  ^  3  '  *’  7 

73  z  <o  •-  .3 

oi  ji  X  r>  •  ii 

i*  c  *-  r  *-•  "* 

M  T  (J  L*  7  ^  ►  ■# 

7  ft-  /  ^  * 

*  J  / 

ft  U  7  C  *-  L 1  ► 

L*  u‘  7*.  >  c# 

ft-  ft  a  ft  *<  ' 

U-<  i*  •  U  1  Uj  “  »  C» 

:  j  c  e  * 

«i  4  X  u  -  •  U  1  4.' 

fr  ft  ft  ^  «*  »  ► 

«  U  7  7  ^  '  >- 

ft  ft  *H  *-  (»  U'  If 


o  f.  —  ft  ft*-  ^  l»' 


c-  *S  r  f  f  C  f. 

fv  r  r  c 


4)  •  «*«• 

in 

«*  y 

pa  0  or 

ftj  H  kl 

^  o  to 

&  X 

O'  X  3 

0  O  2- 

te*  Cr 


PO 

2 

te 

*0  U. 

tei 

tei  o 

-i 

V) 

z 

teJ 

m  «r> 

3 

cr  0 

4 

15 

U  tel 

X 

PA  •-* 

3 

4 

H  « 

i\ i 

a. 

O. 

te  4[ 

O 

0 

«•  • 

<» 

te.  X 

tel 

tel 

IS  * 

tei 

• 

o 

tei  m 

*- 

X 

o* 

•A 

C5 

Ui 

cr 

O 

X 

«* 

« 

M 

< 

•4 

X 

Mi 

cj  te 

15 

z 

te- 

m 

X 

»r 

4 

X 

3 

M 

>o 

• 

*-« 

Z 

tei  te 

tel 

tei 

IT 

tei 

o 

«k 

:»  z 

*-• 

ir 

3 

15 

3 

i- 

Hi 

tf 

M  te« 

• 

r*- 

«r 

4 

3 

c. 

te 

C 

3 

te 

Ai 

X 

*- 

o 

CJ 

H  te. 

B1 

•4 

»-• 

H 

X 

15 

U- 

M 

*0  0 

• 

z 

»- 

te 

o 

X 

x  a 

IT 

M 

3 

cr 

CJ 

»- 

A* 

x  4  uj 

te 

A 

K. 

3 

te-* 

3 

» 

•  15  S 

Z 

te- 

A 

Or 

m 

A 

X 

u 

*  a  -* 

• 

* 

•-• 

Z 

ui 

• 

Z 

M 

hi 

2  Z  tt 

— * 

IT- 

•-* 

te- 

te- 

te 

M 

O 

o  «  Ui  « 

ar 

cr 

Z 

13 

z 

« 

H  >- 

o 

tr 

u. 

tel 

• 

• 

te* 

* 

It* 

te  Ui  4 

• 

X 

rv 

w 

te^ 

r*  « 

»-« 

u> 

O 

o 

•■N 

X 

tei 

*-4 

r» 

*•» 

Ui 

«r  « 

tei  ►* 

• 

z 

te 

r 

« 

CD 

M 

X 

4  4 

3  t-  « 

UJ 

cr 

te 

•w 

•-• 

3 

te  « 

o 

4 

o  * 

tei  4 

0 

r 

A  3 

CD 

cr  te 

• 

•- 

UJ  2 

M  4 

«  u  < 

» 

tn 

te- 

te- 

a 

X 

tu 

X 

X  4 

pa 

►* 

< 

XU 

0 

or 

3 

UI 

o 

o  ♦ 

15 

o 

4  % 

— •  4 

•-  3  * 

r 

z 

IT 

ui 

r? 

X 

<J 

te-  * 

4 

15 

Z  *. 

4 

*• 

Z  te 

X 

•< 

«r 

p- 

-r 

i- 

X 

Ite 

UJ 

cr  4 

rv 

te*  *-•  UJ  O 

X 

3 

Kt 

z 

u. 

o 

•  O  te 

►te 

►-< 

X 

4  4 

K1 

42 

CJ 

o: 

a. 

AJ 

UI  UJ 

» 

O' 

13 

u: 

M  V 

4 

cr  4 

M 

CJ  3  4 

D 

• 

c> 

o 

te 

or 

z 

te*  N 

z 

»  4 

te- 

3  T 

4 

Z 

iii 

M 

c 

M 

te4 

#— 

• 

« 

u; 

3 

te 

X  4 

te  X  M 

3 

15 

cr 

X 

• 

a  ^  in 

te-  « 

r- 

A 

M 

M 

A  4 

te. 

V 

X 

•x 

Mi 

a? 

Ai 

3 

h  te 

te  te 

z 

■ 

*•  • 

cr 

»— »  ^ 

u»  *  te 

cr 

4 

X 

rv  te* 

N— > 

ID 

w 

• 

•  te 

Ui 

Ui 

UJ 

e> 

p^  te/ 

K 

te>  4 

OCD 

a. 

#—  * 

4 

M 

Kl  ^ 

X 

te- 

te- 

1-  te 

• 

5 

2 

cr  p~ 

cr 

te  4 

r 

CJ  •  A 

c  • 

S  4 

*"■ 

«* 

cr 

4 

3 

cr  « 

► 

4 

Ki  • 

■ 

4 

4  4 

7 

*  4  te 

U‘  < 

*  « 

z 

•4 

_» 

z 

tei 

CJ 

X 

U  i 

• 

•  * 

*»■ 

z 

te- 

• 

cr 

u»  ? 

Ui 

X  4 

X 

»:  ~ » 

te 

te  * 

• 

<« 

a 

r^> 

cr 

Or 

te  • 

A 

--  *”i 

PD 

A  A 

o 

ft  * 

v 

X  C’ 

4  » 

A •  4 

3 

V 

>■ 

A* 

o 

* 

n 

*-  • 

M  « 

* 

4 

w 

u.' 

>  O 

V 

C  4 

te 

tei 

^  O'  r 

Q  K  i 

U  * 

M 

X 

te  « 

»- 

«-« 

u 

w. 

u 

v- 

^  ♦ 

V. 

*0 

■  U  1 

* — 

o 

~ 

te  U, 

te 

U  4 

CJ 

tei 

I 

X  tei 

•a  A 

A  4 

►* 

» 

a 

X  M 

te 

0 

3 

X  ♦ 

» 

♦A 

IA 

X 

UJ 

z 

UJ 

tei 

c* 

3 

o  — 

ft 

a 

te.  ♦ 

u 

t 

cr 

ID 

to  m 

IT 

x-' 

•"»  4 

X 

3 

•-* 

_! 

7 

X  a 

1-  7  O 

-j  *“ 

h  4 

I 

0 

te- 

X 

X  te 

fD 

i- 

P- 

i- 

z 

4 

o 

z 

m  <r» 

M 

4  3  U 

ID  ff 

<1  ♦ 

IT 

►-« 

►— 

te 

te- 

a 

o 

O  te 

4 

O 

<1 

r 

z 

15 

z 

M 

i 

** 

te 

3  U.l 

X  •  % 

*  ui 

y  * 

rv 

w 

“■ 

3  IT 

4> 

Ai  K> 

aJ 

CJ 

CJ  * 

UJ 

X 

3 

X 

IA  U. 

o 

»- 

A. 

a 

7 

3  te 

or  *  * 

«■  CJ 

Pr  • 

3 

o 

3 

ir> 

CJ 

« 

A 

A 

c  • 

A 

►-• 

CJ 

z 

r 

►- 

O 

<  > 

O  X  M 

X  LJ 

O  * 

Cj 

z 

«  n 

« 

u 

b 

U  te 

► 

O 

CJ 

n 

O 

o 

> 

X 

cj 

O  te 

u  a  a  ft  «j 

u  4 

o 

•-* 

CJ  c 

Ki 

4 

M 

M  • 

u 

4 

u 

c 

CJ 

Ai 

AI 

ki 

« 

■o 

o 

Kt 

CJ 

K» 

K< 

• 

3 

A  tn 

a 

•*r 

• 

• 

« 

• 

« 

te 

« 

te 

te  • 

-J 

ti_ 

>. 

X 

X 

> 

i 

T 

5 

> 

y  ♦ 

IT 

0 

O'  4 

a 

♦ 

or 

O 

o 

CJ 

O 

n 

C  > 

rJ 

CJ 

c  * 

•r 

K* 

PD  4 

t 

• 

o 

IJ 

•3 

o 

CJ 

(  • 

U' 

u 

CJ  • 

<r 

cr  4 

> 

z 

♦ 

w* 

K>  * 

1" 

« 

« 

« 

« 

te 

« 

te 

« 

te  • 

X 

te 

cr 

c  o 

O 

D  4 

1- 

t.' 

CJ 

A' 

A  CJ 

Ai  U 

u  c.» 

o 

C_' 

CJ 

o 

U 

Cl  CJ 

CJ 

CJ  IJ 

4' 

p*' 

K* 

CJ  u 

l-1 

CJ  CJ 

• 

4 

> 

tx 

.J 

* 

V 

r. 

A 

n 

e» 

ir  ^ 

r- 

«' 

C' 

—  A 

•n 

5? 

cn 

*Tt 

ti 

0 

f 

IT 

V 

«r 

IT 

K  1'- 

IT 

\' 

If"1 

-r  *1* 

J 

-r 

u 

o 

4 

« 

r: 

f 

te- 

c.- 

i. 

e. 

t*  r? 

r 

f . 

C. 

e  «• 

r 

CJ  c. 

r. 

i*. 

te< 

6 

e 

te 

c. 

r. 

te 

r: 

ci  te 

te 

r» 

I, 

1^ 

te  te 

te 

te  te 

te 

70 


00*>T  J  •  rUMITfH  |«  UNIT  FOR  READING  INPUT 

0  TYOF  OflB,  T ,  J,  INFC,  NAME 

?<?«0  F03M*Tf/'-:tMAC-E*,  JJ,  ',  UNIT',  13,',  IREC*',Xa,',  NAME*  ' , A 1 5 

0(>FNf!.lNlT*J,Ni“E»NAME(i)  ,  REAOONUV,  AC  CESS*  'DIRECT', 


F'JPTP  AM 


'Trim 


«HI 

—  tr 


•4? 

G  * 

«n 

U> 

—  *A 

tei 

4- 

UJ 

x  6 

2T 

M 

49 

—  r 

M 

z 

•  • 

_l 

4— 

a 

X  X 

4- 

<*  4 

ft 

UJ 

7 

r»  *-* 

(lJ 

♦ 

ir 

3> 

a  * 

> 

« 

■< 

w 

O 

* 

i: 

UJ 

A-  »* 

♦ 

O 

r>  3)  ^ 

a 

* 

«* 

o  o  ~  ^ 

o 

« 

4- 

*• 

li.  U  ?  7 

<T 

« 

2* 

•  ' 

MM 

_J  _l 

« 

a 

cr  cc  tr*  in 

f 

« 

z 

»- 

M  •  «. 

Z 

•  C? 

♦ 

►-« 

r> 

n< 

V  X 

_ I 

r  —  z 

* 

A 

• 

x  ■»  rv  ru 

49 

e>  %  iaj 

« 

_j 

4- 

• 

m 

7 

«o  • 

* 

n 

u 

ir 

GK  59 

►  * 

o-  U.'  « 

« 

-«* 

O 

a* 

A* 

*>«:*-.*-* 

* 

z  — 

* 

o 

ft  ©*  in  in 

— « 

UJ  O 

* 

UJ 

UJ 

• 

■ 

«.  «■ 

• 

a  c 

* 

n 

«o 

•M 

u- 

WWG  - 

1C 

»- 

« 

o 

o 

4— 

p«- 

®-  O.  ^ 

»■ 

4-  * 

« 

_j 

_J 

~J 

•A 

*-  >  K  H 

#-* 

« 

o 

o 

3» 

* 

»-  4-  ^ 

9T 

4H 

« 

*) 

• 

*— > 

r*  M 

♦ 

* 

4. 

UJ 

O 

r*  «-•*  « 

4- 

_>  * 

« 

—• 

— -• 

•n 

cr 

z 

A-M  a 

3) 

•  • 

« 

w-> 

♦  ♦  K 

C 

P  C 

♦ 

mm 

#-» 

• 

4-4-3) 

a 

uj  or 

* 

a 

• 

a 

a  or  O 

•  c 

« 

•A 

IT 

2T 

4-  4-  ti. 

<c 

Ci  o 

* 

*•» 

IT 

a* 

•>  ft  3> 

G?  UJ 

« 

4- 

AJ 

4- 

4- 

cr  a  c 

—  a 

« 

w— 

•  -J  M 

*-• 

*-> 

m 

« 

#-» 

4- 

r1 

■  #  • 

—1 

K  ' 

« 

4- 

Af* 

*? 

■ 

4Jli>*-*- 

« 

G  x 

« 

7 

M 

Or 

* 


»  •• 

n-  ft 


0. 

I 


— »  _f  m-  4- 

•  •  ■<  -a 


X.  41 

_J  I 
^  n 


a-  r. 

a  i 


*.  r 

c  ** 

R  a  r 

*  O  n  £'  cr  < 


z 

o 

«_) 


15* 

G» 

er  ps 
<r  <s> 

u  c  c  <r  ' 


r. 

e. 

O' 


* 

o 


U»  M  z  u1  £ 

2  i>’  ^  /  •- 

►  +-*  u* 

*-  a  »<'»  -!• 
z  a  e  z  c 

c*  -»  _*  o  -J 

4_  C’  *“>  4.1  C'  L. 


1  7 
U  UJ 


r  e  ^  • 

F-  fi  R  P. 


73 


IV-PLUS  J  7 1 3?U  3  12-DEC-82  PAGE 

0PL7.FTN  /TPJPIOC^S/WP 


?7371fc  77777? 

77377?  77.71772 

77  71  27  7  .'7777 

77777a  7777 7? 

77.27  1  ?  P-VUCI 


1 3?  *  1 3  12-0CC-S2  PASS 


Ju  1  iv 


X  '  I.  <  . :  •  £  '  A  Xrffis  ‘  X  A  fin '  AhHi  *  In.  IS  8  -  "i »  -a 

fX3a'evMd't‘QIilc‘fij(i'»HdT»'f.»3«A'\v  )«»  rr  *  1 1  3  a 


FCUTA^Lt  C0'">f  --  I  V  JT  T  ‘LI7ATI3M 


'• A ”E ' i,  n  • 

iMlMtrjl  *  8L4NK 


FORMAT  (SX,  'ENTER  COMTAL  IMOEX') 


I 

o 

7 

~z  a 

-  »  C' 

3  a 

v  —  u  I 

a  i- 
Ui  uj  a>  i 

a  ui  uii 


J  3  3  1  >. 

o  c  t_  c.  :j  i 

CJ  O  U  U  LI 


U  U  U  C  U  I 


#-* 

« 

II 

feJ 

*— * 

7 

Cl 

«• 

_J 

ft 

X. 

IrJ 

c 

a 

CJ 

<t 

x 

* 

»-« 

'■ 

«k 

•  1 

_i 

•-< 

ft. 

7  - 

• — 4 

t- 

r>  x 

» 

X 

•-«  o 

c 

O 

CJ 

u 

CJ 

H 

X.  * 

O'  UJ 

CJ  M 

5. 

Ik*  nr 

o 

c 

*-  lA- 

*•* 

o 

Z  X 

*.  — 1 

*- 

UJ 

M  O' 

UJ 

»  c 

<r 

Cl 

•A  lC 

-  O' 

x  Cj 

«-•  » 

«  1 

tv  uj 

O'  K1 

•V 

^  U-* 

<r  »- 

►-  Z 

cr 

c- 

< 

UJ 

«« 

X  V 

z 

tu  ►- 

2 

ft  z  c; 

OL  ►» 

o 

t  «*  Uj 

II 

>-  ft 

o 

L_  O  34. 

h  3 

u. 

—  CJ 

»-S 

ft 

ft 

o_  r  » 

•-» 

a.;  a. 

<  •  a 

X 

O 

Cl  u> 

o 

Cf 

CJ 

C 

O' 

r*- 

c.  c 

r 

—  CJ 

* 

ft  CJ 

c.. 

;  >  j 

1  5 

2 

i .  rj  t 

c  . 

c  C. 

ecu 

c . 

ft 

C’ 

i  c.>  cj  cj 

CJ 

O  CJ 

c.» 

<  M  <T 

z  k-  ».) 

►-  ^  V) 

r  > 

u  u.. 


cj  r.  c 

Cl  Cl*  l 


*  ♦  CJ 

*  ♦  e  — 

*  *  r  tv 

i.  o  c  ci  u  ^ 


»r  a 
r*  k  **«  r- 

e  h  t  -  cs. 

G  &  «v  M 


TYDF  *■"  i ,  y  V  f  I) 

FORMAT (/?*  ,  »ENTE»  SCALE  SELECTION  fc?  * ,  i  , » - ax : s •  / 

rx ,  'rFTin-  ally  follcaf?  f*  l t*'ER  »^r  upper  lI"'TS’/ 
?V,*?f  0  A  ■'  A  TO  PLOT  E.3.  "!,0.25^"'/ 

?x,'cwnrF  of  scales  s«ci  ,r>  te  i  t t r a t e o  ext'/ 

?X,'  1  »  .5  Pfc  T£P  RtSOLL'IOM,  PA.'  DATA  SPREAD  p5?'/ 


2ERJM. 


N 


M 


is 


&  &  G1 

sain 
<v  y* 


Ui 

<  y  < 

qc 

o 

U*  UI  UI 

Ui 

•< 

a-  a  a 

a 

CL 

ana. 

a 

rj«  « 

3 

<  <  < 

K  K  H 

ft 

tn 

P5 

<  y  y 

ft.1 

c 

GOO 

X 

« 

-4> 

o 

V  X  X 

_? 

•i  i 

C' 

y  y  < 

y-* 

*- 

•—  m  ir 

0>. 

X  X  X 

Ui 

% 

O 

• 

nrr 

O 

> 

*- 

LO 

v- 

13  ft’  ft 

■£> 

A* 

»  y  ^ 

y  »~ 

O 

r 

40 

7  2  7 

X  o 

_> 

•-* 

-4 

K  H 

o 

9 

o  o  o 

_J 

<L 

m 

ft 

*r  is  is 

f- 

Cl 

ft  ft 

•r>  *  • 

it/ 

o 

3  3  3 

o 

O 

»- 

• 

C 

r- 

<U  rs  «*' 
•  y*  »-« 

2 

• 

3  -J  •-• 

►- 

K  *-  *- 

A 

o  ::  o  *~ 

■< 

e 

■ftt 

1“  I  J 

y* 

co  lj  if,  x 

C3 

•< 

»- 

*- 

#  *-4  H 

►— 

Ui  Ui  UJ  M 

ft. 

O  •- 

y£ 

yi 

_1  _J 

3 

cr  ft  a  3 

«c- 

A/  -ft 

O 

►y 

O 

►y  c  a 

M  U« 

a_ 

ft 

*c  o 

*->• 

^  JD 

•-  3  K 

»- 

tr  a-  a  2 

o 

ft 

X 

ft 

X  •  • 

^  y  k 

3 

O' 

lu  y  u» 

»- 

o  ft 

o 

*-• 

o 

Max 

o  <c 

o 

in 

K  1-  H  - 

o 

1“  o 

—1 

■3  0  0 

jWr- 

•• 

Ui  Hi  ai  y-> 

13 

►- 

o 

1- 

o_  •  • 

a  *-•  c. 

CO 

L  L  L  M 

o  •> 

►y 

-i 

H 

3S  6 

3  W  ^  H 

H 

<\i 

r>  r> 

lO  »~ 

X 

• 

X 

•  •  • 

•  X.  o 

X 

•* 

a  n>  ^  j. 

<1  2 

►*• 

►y 

ft 

*-• 

ft  ►-  K 

i-  i  ♦'  e 

-y 

*  *  »  M 

*-» 

-1 

l-  to 

’  *_, 

• 

H 

•  •  • 

•  c  «-> 

> 

«  ft  ft  o 

e.  -ai 

6*  _i 

OC  w 

& 

ft 

a  «:  «-» 

w 

is 


3 

G 


O 


l&i 

*9 


O 

O 

w» 


7 

o 

13 


O 

a 


7 


tO 

3C 

U 

k,  r 
—  «* 
If  x- 

fv  Or 
<S  f- 
>  >» 

•0 

.J 

a 

i 

>  2 
v->  ►- 
U 

ft  * 

-rt  »- 

a  *- 
»  « 
tf  11 

Cj  »•> 

U  lit 


O  tf  J  z 
-  C 


U-  <* 

K  Ui 
•  Ql 


ft  >.  • 

:»*-*►- 
0  3  3 

_j  r.  • 


ft  X 

ft.  w  w 

t  j  JI  X 


*■  M  >  »  IT 

ft  r  «.«»»- 


*  a  u  •-»  zs 


*-  b-  »■* 


a  *«i 
ui  v 
L>  ft 
Cl  £3 


•  ft  ►- 
i'  JP 
*—  <j  ft 
ft  CJ  »./ 

h.  u>  a 


i  Q  f 
»-  2  f. 
,J4  X 
0 

3a  u 
wk;  a 

It  J£  >- 

•“•  c.  »- 


>«  3  •-  M  9 

A,  —  x  —  AJ 

H  L  3  ■£*  *--  *-• 

<  g  o  <  a. 

I  m  J  U  V  ft 

a  o  •-  a  i'  u 

G  ()  U-  >  O  U 

U.  «k  M  I  -  L.  < 


h  c.  a. 

x 


►-  o  z 

»•■  14.  I  *— 

Y-  ft 

►-*  O  »-*  t*> 
.,i  n  ft 
l_  .  uj  «i 

u  y  or 
a.  *+  - 


u.  o 

^  * 


•  t/5 
OC  ft  I 
O  <l  3 

•  v- 

■o  ft  1 

O  l 


ft  ft 
X 


s  I 
c»  -o  . 


y« 

•>*  M 

«4> 


& 

*V 


X 

o 


»o 

X 


ft 

z 


c 

U' 


ft) 

cr 


3 

O 

X 


•o  • 

*-  tn 

X  7 
I-  O 

C 


w  tf  c.  a  ro  ai/  ucr-eior,  *:f\, 

K  N  o  v  (<  c{,  *.  «r  c,  o>  t* ■  t*  r*  (j  fr 

6  t  M>  f,  b  K  M.  f.  C>  fv  t  r«  <- 

«>  f>  e*  g  c>  o  &  ft.  ei»  c*  *♦ «.  rv  a 


^  Cj 

cr  a 
C.  C\ 

&  & 


ft 

K 


O 

o 

40 

ft 

4 

o 

ft 

3 

3C 

k-  a 

o 

• 

o 

■4 

IV 

»- 

_r 

r 

Ul 

Ui 

Cl 

Ui 

-» 

1 

z 

-> 

u.«  ft 

Y 

Ui 

M 

iis  r- 

*  * 

*— 

»*- 

U- 

h  »_ 

> 

-y. 

... 

_ , 

o.  y 

O- 

X 

3 

Z 

a 

*-y 

lO 

o 

/• 

_ _  , . 

c. 

ft>  tf » 

IS 

M 

3  Q 

ft 

Ui 

C 

>■ 

M 

ft 

a 

a 

_l 

«-« 

z 

3  3  »-y 

M 

* 

o 

ft 

o 

y 

ft 

V- 

G  -* 

a 

V' 

cr 

-  3 

u.. 

7 

»-■ 

> 

V.‘  1- 

o 

ui 

O 

ft 

o 

3 

U  *  CJ 

u 

cr 

iT 

I- 

O 

ft  y 

to 

y 

U 

o 

ft 

*-< 

k~* 

y-s 

13 

3  u 

IV  ft 

>-  •  - 

ft 

i 

UJ 

— 

. ' 

v  .'  - 

ft 

M 

Or 

ft  xi 

o 

ft 

X 

ir 

O 

y- 

*4  ¥~ 

c» 

>- 

*  c- 

O 

w 

> 

X. 

ft  i- 

UJ  ►- 

tn 

3 

C 

to 

•  o 

>  O  ft 

ft 

V 

ft 

H 

a  c* 

-«  13 

> 

C 

i* »  C 

J 

• 

c 

X 

3 

a 

« 

a  c 

I* 

» 

»  c.  r» 

O 

r. 

• 

_ ( 

X 

=1 

y-  -  - 

L-  ft 

a 

f' 

> 

c* 

ft 

OJ 

A 

—• 

•  n  * 

<.  • 

« 

>. 

rc 

ft 

o 

- ' 

o- 

4-  it 

r  ft 

r 

r* 

»— 

cr 

y* 

1- 

3  .  /  ft 

»-  ft 

c 

«— 

►- 

C 

k- 

»-■ 

»». 

•  o  • 

• 

ft 

<3 

c  , 

1- 

ft 

> 

U' 

y  a  .-r 

>  »-• 

« 

*— 

> 

fly 

ft  u 

7 

U 

fr 

t-  ■ 

»-  ;  *-• 

• 

i — 

0 

i 

1 

> 

O 

> 

O 

o 

ft  '  -  ft 

uL  It. 

ft 

a 

G 

►  s 

o 

»- 

ft 

yi 

^  *-* 

1 

;< 

ti 

c: 

.-J 

«r 

o 

a 

40 

A 

-t  u 

U 

r.- 

u.  ►- 

3 

o 

X 

u  ^ 

ft 

*  • 

• 

• 

♦ 

k 

r  ■  3 

3 

ur 

♦ 

*  ♦ 

• 

« 

* 

cr 

O 

♦ 

♦  « 

*» 

* 

o 

o 

Cj  U 

1J 

CJ  O' 

CJ  t-v 

cj  cj 

C-> 

ft 

x 

tc- 

cr 

(.  -< 

.*r 

O' 

o 

r  f-. 

g; 

o 

t. 

c 

w 

ft 

t 

c 

ft 

&  ft. 

81 


84 


JT  TS  POSSIBLE  That  TME  USES  DELETED  n.L 
pO T  T  0M  SELECT  TYPES,  IF  so  THE  LEAST 
V'  CA'-'  00  I  ?  WAS  V  HIH/-ES. 
fvO  TEP«I*»Tfc  thc  piwrT 


«\J 

«& 

I 

o 

U. 

a 

t 

a 


*r 

<v 


«  «  «  «  * 
«  «  «  *  « 

u,1  o  y  lu  y 
«  to  (T) 
Dtr  D  33 
!«• 

>  Ci  >  X  x 

to  *  to  to  ir 
r>  *  r*  3  3 
to  *  V?  to  to 
Q;  *  ^  o-  or 
U-‘  — •  U  :  U:  UJ 

>  >  >  > 

X  *->X  >1  M 
to 

O  '■OUlS 

r>  *  r>  r»  33 
tr<  ••  tc  id  n. 

tu  C*  Uj  Uj 
O  U  Q  Cl  O 
A 

*  to  c  *  « 

*  kt  *  *  « 

«  X  «  «  « 

u 

*  X  *  *  * 


•O 

a. 

H  f\i 

—i  _nn 
u.  w  x# 

5-  CV 

x  ru  x  * 
I  _/h% 

XU.  H 
V  «*«Y 
•V  X  » 
X*  -J  m- 
O'  »  U  _J 
U  >  U. 

H  J  X  > 

?  k  m  x 

u.  >- 

X  x» 

*  x»  —  « 

*  o  « 

*  *  « 

K  C  ^  (5  K 

o 


K»  w 
O  h  h  H 
<  ft  tai  « 
uj  r  iu  k  £ 
ft  a  u  »-•  at 
>  C-t.'  a  o 

X  U  «x  ft  fc. 


u 

UJ  UJ 

UJ 

to 

to  t/J 

to 

3 

3  3 

3 

X 

X  X 

X 

CO 

a?  co 

to 

3 

3  3 

3 

to 

ft)  CO 

to 

a 

tt  cr 

or 

UJ 

Li  j  UJ 

UJ 

> 

>  X 

> 

X 

X  X 

X 

u  o 

(9 

3 

ry  t> 

3 

C£ 

tt' 

UJ 

UJ  UJ 

UJ 

X 

o 

O  Q 

o 

tt 

t>- 

■n 

ft  ft 

ft 

r-\ 

5. 

« 

ft  ft 

ft 

k- 

2 

ft 

ft  « 

ft 

Qr 

ft 

UJ 

ft 

ft  « 

ft 

> 

— < 

— • 

— •  — • 

2 

ft 

ri 

^  r' 

•— * 

y 

ft 

•“>  *>  •"»  ^ 

*■) 

•>  3  3 

3 

•L  * 

«-  •> 

•. 

» 

«» 

w 

ft  I 

ft  ft 

X 

ft 

5 

ft. 

ft 

X 

3 

»-<  w 

v-> 

k— 

X' 

'-' 

ft 

»- 

■—•XX 

X  X 

X 

X 

X 

X 

a 

ft 

lO  M  M 

M  H 

M 

M 

M 

H 

X 

V-J 

r. 

— 

>« 

ST-  It  »• 

It  ■ 

ft 

tt 

tt 

tt 

«. 

ft 

X* 

ft 

w 

O 

^  O 

*— * 

w* 

X  3  — > 

3  — > 

3  -> 

•“? 

•-> 

X 

— > 

»— »  VJ 

VJ  w 

*-< 

»— 

w 

x » 

X* 

OZ2 

X  X 

2 

2 

X 

X 

X 

w 

W  *-*  t-t 

1^  v-A 

x* 

X 

X 

X 

>- 

-J  .J 

a.  a 

__i 

_j 

L. 

<4. 

2  X 

2  X 

i; 

X 

«l 

X 

-> 

2 

*•-. 

r'  x*  ** 

»—  ft 

•  - 

M 

"ft 

w 

•- 

X 

t.  ft  y 

ft  1 

2 

2 

ft 

ft 

X 

3 

2 

• 

Q. 

• 

ft.'  X  r^r\r>4 
ft  Cl  r>  r-  <>S  O 

•  3  ~t  3  ->  rv 


tt 

Xn 

ft 

2 

V 

£> 

to 

fO 

tt 

2 

tt 

0 

cr 

u 

t. 

«. 

U 

0 

IT 

3: 

»* 

K» 

—• 

_l 

•ft 

2 

cr 

w 

ft 

a. 

XI 

ui 

• 

•• 

O 

X 

UJ 

X 

A 

ft 

X 

2 

0 

l. 

X 

O 

O 

X 

ft 

X 

o 

u  x 

ft 

-J  Z  2  X  X  tO  uJ 

t>  *>  »-*  H  *H  c-  3D 

x  _j  _  a.  a  2 

x  ?  v  2  x  on 

M  M  4  H  «t  H 

5  I  I  2 

l.  ••••DC 

•-<  X  X  X  X  13  o 
-jtr.it!> 


jr-s  r'  r>  <-» 

#->  #-  r> 

3  3  3  3 

/  2  x  x  iU 

m  k  ►.  *-  3 

j  j  p  a.  ir 
2  x  /  X  M 

v  >  s.  x  2 

r  *  *  •  UJ 

I-  I-  l-  u 

t-0  _t  o 


2  *1  3  1 


O 

O 


o 

o 


L  L  r  5 

>  X  X  X 


UlLU.li 
»>«  ».<  M  M 


3  3 

ft  s  ? 


X  *»  >  > 

M  M  H  XI 


a 

X 

X 


>  X  3  3  X 


*■  c  *s 
'  X  3 


UJ 


Cr¬ 
uz  2 

>  M 

Z  _J 

X  X 
‘  *4 


‘  x  *. 


ft  •.  t 

T 


*■  & 
>  3 

3  Z  r-» 

•  -*-!►- 
^  X  _J  «s 

^  Z  K' 
K  >  H  C 


U  fi 
n  k» 
2  O 


tt  w 

1  U*  2 

>  X* 

Z  -J 

X  X 

•  -4 

—  X 
tt  «  t 

ft"  »•• 

>*■  3 

3  Z  ' 
X  -J  » 
X  x- 


IO  x 
lO  <J 
UJ  X 


U  u. 

X*  M 


Z  CL 
O  v 
U  ►  - 


rv  K» 


-J  * 

7  x 

1  cr 

*"  — 1  e. 

X-  •»  w  pr 

-J  LV  :.'  rr  <y 

X  «  2.  ul  c 

x  L  •<  K' 

O'  u  x  *_  O'  «— 

to  _J  Z  x.  UJ 

UJ  c  O  K  UJ  ►- 

if  <3  U  2  (i  •  - 

t~*  ova 

-J  4_J  X  ft 

u> 


_»  r.  _j 

~  K'  « 

€L  U  'v  O 
—  *  "v  t  , 

K.  x  v-  K 


u*  #-  ft  a: 

fl  ►  •  iv  o 
>  O  O  v 
►-  t  it  »- 


c  CJ  c-i 

«-J  —  C-  O  o  C- 


o  c,'  o  f-  r_  c  o 

U  ’J  o  t->  t  ;  l  ’  U  i 


1  A  O  O  O 

c  COO 


HITD 

O  X  O  o 


o  c  o 

tJ  C  f.  i  t 


c  t  .’  t_»  u 


A  t- 
X  ^ 
A  U 


t*> 


A« 

«  S<  f: 


:J  IT 

«•  t'  *<  Ct) 


ft  r 

tf  p 

r. 


r.  —  a  f- 
C  o  (  (' 


85 


TC1  PLOT  PtpLfCTANCE  TTP  L 


2erjm 


A 

% 

•0 

kl 

CL 

>- 

fr¬ 

ill 

U 

Z 

<< 

A 

4- 

«G 

O 

1 

kl 

l.> 

_i 

Hi 

k 

O 

WJ 

• 

Or 

A 

■ 

»* 

1M 

r-> 

— • 

«r 

o 

(T 

t~ 

4~ 

in 

? 

O 

• 

*• 

3 

o 

ki 

•* 

O 

3 

A 

U 

*■-» 

cr 

** 

M 

A 

y- 

►> 

• 

• 

•-» 

4- 

U 

•r 

4- 

o 

z 

2. 

tt 

-J 

t 

y~ 

a 

a 

Wj 

V. 

% 

z  c 

tn 

«. 

«  ■ 

St 

X 

o 

r 

u 

fri  A 

_«  • 

• 

it.  c 

--  *-• 

y- 

«* 

— 1  -  J 

»-  f- 

y-  -> 

if  J 

in  a 

«*  r> 

o.  • 

•  & 

•  •• 

hu  > 

UJ  7 

A  a 

(i  U  » 

IJ  *- 

»-  M 

K- 

►  o 

U  U 

it 

>  S 

4-  k  2 

^  ►- 

•"*  z 

to 

3 

-i 

a. 

•* 

* 

i 

r-  i*>  * 

* 

>  x 

* 

»>'  t- 

—  ~  t-> 

u 

i-> 

?  • 

«  4- 
U  t- 
^  *1 

CT  l/l 

*1  •». 

cr 

o  «. 

C  C 

©  .' 

r  <•*. 

...  V-. 

v  *  •■  . 

W 

V*  It* 

r-w  «. 

a  e 

C. 

X 

Ai 

_J  IT 

C 

B 

*  —> 


C» 

•Q  «  J 

►-  Z  O 
Vy  »-•  O 


a  •  ^  7  • 

^  o  '■>  Ai  *-«  a  a 

—  U.'  *■’  w  k  UJ  o 

*A*  w  I  *  •■>  «  u. 

4  r>  I  M  h  ^ 

*"*  •  •"*  %ft  I  O'  P  r\i  UJ«£Ci' 

«-» W  _J  a  <j  —  W  T  tn 

•  ^  7  u_  _>  0  3  l-  U'  «  . 

*  -J  -J  _J  «0  4  •-  __r  7  + 

K‘ w  <  i»j  «*  m  u  ^4  aj  »* 

V  l/)  »•  l_>  3  >  »  *“ >  *-  CA  l.'  •_  UJ>-»~ 

X  HtfJZ  j 

H  X  O  *"  M  it  ^-<OS(V-liXwM  H  K1  |l 
4  li.  ►*  M  {T  M  /y  <T  u.  K  Jj. 

*  >-  ►“•»-•  Z  >  ^  7  »•  *-  X  »  «  J*  v 

o  <  4-  *-<  O  O  14/  o 

u  •  nu.  ho  u  le¬ 


ft-  f\  fn  st  »r.  f, 

r~  c  f  ;»  .  .  r  r.  r: 

r,  A  rv  i\;  a.  a  a  a* 

f'  t  (v  M5j  K  f,  A 


«i  ff  f  «  rv  r,  d  IT 
.•  -•  ••  ~ 
fl  (V  fr  fi  f  fv  a  A 
tv  6>  t  R  14  H  t.i  (> 


A  A  f..  A 

A  r;  i  r.. 


Jim 


on  tiJ 
W  _» 
2  U<  O 
*l>  J  T 


O  O 

x  u.  <o 

^  k- 
C  M 

fr  2  t 

O  O  M 
t  13  J 
>- 

»“  Jft 

«•-*  c_- 

I  ft  O 


O  —  i 

c  • 

A-  o 

U  I 

•  I 

o  K  I 

»-4 

t-  a  - 


(t  O  UJ 

<  >  —  o 

k  k  > 


> 

O'  ~S  N 
O  %. 

U  o  0-* 
*-  —  f\J 


ft  ft, 

Ul  *fcj  o 

a  o  ?  k 

>-  u  r 

►*»>•<  u  ^ 


ki 

c»  o  o  r  r  <  c  l  toe  cj  «' 
O  CJ  1J  U  L'  u  u  U  O  IJ  <  , 

•  *•«««« 

UlJUUL  tJUl.lL’  O  4-»  C»  U 


»  4C  • 

V 

La, 

4 

LaJ 

►"  P~ 

• 

2. 

4 

H 

o 

►*  3  >-•  ~ 

4 

*-• 

4 

llu 

>  -f 

o  o 

4 

4 

u 

_r  X 

*  t  k>  -  r 

♦ 

U- 

4 

V) 

r>  _j 

*-  O  W  *_  o 

# 

o 

4 

ri  *.■ 

DUK  DU 

H 

4 

4 

UJ 

fv 

O  M  o 

* 

13 

4 

—I 

C j  4  4 

i  n  . 

* 

2 

4 

o 

f\  —  — 

D  M  D  O  W 

4- 

* 

*-« 

4 

X 

•4 

L)  m  y 

M 

4 

r> 

4 

:* 

»•  r> 

2  *• 

« 

to 

4 

c 

*  *"•  m 

r>  x 

4 

UJ 

4 

13 

►- 

K1  O  -  M  O' 

*~*  o 

* 

u 

4 

2 

a  ^  O' 

* 

o 

4 

M 

o 

•  '  -  c 

t 

*— • 

—  •-  r>  ~  • 

>  - 

4 

a 

4 

o 

r  wo  cr  jo 

o  p.' 

♦ 

4 

f>-  v>  j  cr 

u 

* 

4 

III 

O'  o  cr 

4 

4 

p 

CJ 

Ul  U  UJ 

•  •*' 

♦ 

4 

* 

• 

a*  ►-  K 

r.  <r 

* 

4 

tv 

►- 

O  ►'  T  ft  H 

<7 

4 

4 

"X 

r> 

>  tv  t'  V  ft 

«■  cr 

« 

♦ 

U 

1-  i  u  ^ 

•-  * 

4 

4 

O 

_  I 

*— 

cr 

♦ 

4 

U.' 

tJ  ^1 

u« 

<T  *-• 

♦ 

4 

I 

ft 

e> 

nr 

4 

4 

»- 

>4  4 

UJ 

* 

4 

U' 

U 

_j 

uJ  r- 

4 

1 

2  ft  7 

_  i 

ft  ►- 

4 

4 

_< 

*-  7 

-«* 

>-  tt 

• 

4 

u. 

LUJ  C.J 

J 

4 

4 

ll> 

r  c  c>r> 

* 

4 

O 

L’  U  ft  U 

• 

4 

*-i 

«  «  *  *  « 

4 

4 

c- 

:'  y  :  j  : 

4 

4 

Li¬ 

t  j  c  t  or- 

* 

4 

ft 

*  ’  t » tj  o  u 

u 

4 

4 

o  r 

4  4 

4  4 

4 

CJ  O  1J  C‘  C>  o 

u  c  u 

L»  1 1 

O  U 

U  l 

«..  » 

til'  U  111! 


fv  s»  \r  «  ■  *-  r»  r  f 

f'  n  M'  n  n  I  ^  •« 

<V  <v  tv  <v  f  1  *  •  r. 

f*  Ci  ft  ft.  r-  r,  p  r 


TA* 7  _CC  ’  TO  IMirjr 


-press  PQIVGP\  4I?£AS  TO  fHECK  HHETHt  9  THIS  PIXEL  N  13  IN  THIS  *OL VGON 


'Trim 


uj 

m 

3 


Ll>  uj 
c n  o 
3  3 


X 

>- 

X 

X 

c 

_j 

</) 

O 

to 

Cl 

3 

3 

Zli 

3 

-j 

«o 

(TJ 

to 

C/3 

UJ 

or 

O: 

Q 

a 

X 

bJ 

UJ 

UJ 

u> 

> 

> 

> 

> 

a 

X 

x 

X 

X 

u 

UJ 

n 

7 

a 

C5 

tr/ 

u 

t—4 

r> 

r> 

3 

3 

o 

-J 

rti 

30 

rr 

O' 

z 

UJ 

UJ 

u. 

(t » 

UJ 

to 

o 

O 

o 

o 

M 

c 

X 

* 

« 

* 

* 

X- 

X 

« 

« 

« 

♦ 

« 

« 

« 

♦ 

a 

z 

X* 

o 

« 

« 

* 

• 

V 

— » 

— • 

— • 

— • 

Vi 

X 

UJ 

X 

•* 

r 

-J 

-1 

o 

r 

ft- 

X* 

z 

*-* 

X 

X 

fO  ^ 

x-s 

z 

UJ 

to 

<o 

•r 

X 

M 

o 

<r 

a 

K'  rv 

a 

n 

_» 

«r  ^ 

*— 

X 

UJ 

o 

CD 

x 

x 

o 

3t 

X 

X 

X 

O  Or 

tr 

z 

M 

o 

o 

>  UJ 

UJ 

a 

CL 

e 

o 

> 

> 

_l 

o 

n  z 

z 

Of 

m- 

tr  » 

- 

X 

r 

<- 

./  « 

#— « 

#~s 

M 

o 

X 

rv 

ftl 

n. 

a. 

a 

o 

z 

%_y 

*-* 

®  - 

X 

3 

X 

X 

•  •-« 

u 

— J 

n 

X 

►— < 

X 

•-• 

o 

u 

X* 

—1 

3 

15  >- 

X 

a. 

X 

G 

O 

•  X# 

ft 

M 

« 

* 

*-* 

3 

3 

rv  • 

3. 

• 

— • 

<d 

— . 

_J 

• 

• 

_J  ~ 

PA 

Us 

Jj 

X 

X 

k  C 

cr 

~  ft 

z 

O' 

e» 

3 

>-  * 

V  • 

z? 

*-« 

u> 

• 

• 

x 

r? 

Xt 

•-f 

T 

ft 

A' 

rv 

X* 

> 

>  ^ 

\- 

X 

(J 

7 

z 

*J  X 

x 

7 

x« 

or 

• 

• 

Rv  f 

C’ 

o  »- 

O 

<s 

X 

ft 

or 

tx 

x  - 

o  •- 

O 

K 

o 

r 

J 

z 

K> 

cr 

r-< 

UJ 

• 

■ 

». 

»~N 

u* 

*- 

z 

Z 

1 

/-a 

Z-' 

Z 

_J 

a 

G 

i 

» 

a. 

X 

— • 

x> 

• 

X 

X 

♦ 

rv 

V' 

Z 

UJ 

o 

UJ 

»-* 

Z 

X 

Z 

x» 

M 

ft. 

U 

X 

b 

u 

u 

*— i 

X 

> 

T 

•" 

" 

T 

_• 

X 

X 

X 

r-~ 

X 

a. 

O' 

X 

ft 

r 

M 

x» 

X 

G 

«. 

z 

M 

M 

Xt 

3 

a 

0. 

m 

c 

►— 

• 

• 

-i 

_J 

M 

C 

z 

UJ 

B 

M 

X 

>- 

c 

r 

UJ 

3 

3 

II 

O 

e» 

a 

a 

r 

• 

o 

z 

x> 

• 

• 

• « 

7 

u 

X* 

> 

— 

a-- 

a‘* 

_i 

X 

*- 

> 

c 

<x 

z 

Z 

_j 

_  i 

-J 

t r> 

/ 

— 

rv 

X 

*— 

-t 

c? 

V 

_i 

o 

7 

z 

7 

<J 

u- 

u 

u 

a 

*-+ 

u 

t_>  X 

—1 

_ 1 

r 

X 

►- 

a  x 
Z  Z 


rv  —»  . 

Z-  7  . 


7  t  <  ♦  K 

7 

♦  M  V  >  M 

>r  >  > 

>11 


u 

t/; 


rv 

z 


r 

Cl  I 


3 

D 


C/~>  IT 

M  W 

n  » 

y  > 

>*X»*> 

3  7  ?  :• 
«r  o  tA  Z 


3  G 

1  7 

U 

u 

X 

>  7 

>  (  „ 
n. 

*/ 


o 

>  ♦ 

a 

ID  — 


rv  ^ 
pa  to 

•  > 

h  Z 
15  D  ^ 

•  C 

x  Li  ♦ 

rv 

z  ■  a 


o  tr,  n  ^ 
c-  >  _i 
^  <  c 

Z  >  X 

-k.  c  u  c. 

h  u  >  r 


o  u  r 

►  -  c. 

<••  >-  ^ 


I 

>  z 

X  b 


o 

v  o 


OC.CG 
lj  l:  v.  o 


UJ 

u 

a 


a 

u 

Uh 

i 


(-)(->  IJ  o 


o  u 
U  to 

U  U1 


rv 

(y 


rv 

& 


V  ir  i r 

rv  rv  rv 

r  r  6- 


\r  z' 

C’  1/ 

a,  rv 


h-  eC  cr  r.  -- 

v  tr  ir  ^  j 


*V  rv 

f.  ey 


\  At'  A 


IT  X 
7  -r  j 
A  A  A 


89 


wE  WANT  TO  »E*-IECT  SVH30LS 


m*svi..n?)  «  symbol ;polvto*«) 


7 

U 

to 

« 

ft 

ft 

O 

2 

Ui 

ft 

ft 

e> 

2 

♦ 

♦ 

ft 

V 

U 

M 

ft 

ft 

_i 

*> 

_J 

♦ 

ft 

« 

a 

ft 

ft 

ft- 

a 

» 

UJ 

2 

« 

ft 

ft 

«r 

2 

<2 

♦ 

T 

ft- 

LJ 

ft 

ft 

ft 

a 

O 

ft- 

I 

n 

ft 

ft 

_J 

X 

2 

■7 

r. 

*» 

♦ 

ft 

ft 

«• 

r 

ft  ' 

O 

k-i 

V- 

Li 

lU 

ft 

ft 

N. 

L« 

X 

O 

K 

cr 

12 

« 

ft 

ft 

♦ 

*-> 

fL 

X 

ft- 

ft- 

*— s 

O 

■< 

ft 

ft 

nj 

2 

a 

H 

2 

% 

% 

*.> 

X 

« 

ft 

ft 

X 

■V 

_j 

X 

ft-  M 

-V 

• 

M 

ft 

ft 

* 

« 

rv 

-*< 

O  ft 

2 

co 

ft 

ft 

ft 

ftC 

X 

ft-  ^ 

M 

O 

2 

OJ 

ft 

ft 

ft- 

2 

« 

4 

♦ 

2  4- 

O 

O 

— • 

ft 

ft 

ft 

UJ 

X 

JMX 

ft 

>- 

L5 

tr» 

ft 

ft 

CO 

00 

X 

X 

Ui 

•  a 

co 

_J 

>- 

ft 

ft 

ft 

I3» 

ft 

X 

ft 

u 

CO 

<v 

M 

♦  •-» 

UJ 

Q. 

C? 

_j 

« 

ft 

ft 

1 

I 

a 

X  T- 

X 

n. 

<« 

ft 

ft 

y- 

X 

X 

ft- 

- 

« 

-V 

■* 

q: 

2  X 

0l 

M 

(X 

rr 

ft 

ft 

X 

X 

cr 

u. 

—*  O. 

UJ 

*  UJ 

ft 

« 

*- 

ft 

UJ 

Tj 

*— * 

> 

X  X 

CO 

2 

> 

ft 

O 

ft 

2 

«0 

cr 

n 

r' 

f? 

»  *r 

♦ 

ft- 

r' 

• 

0 

ft 

ft 

.  « 

* 

»-* 

♦ 

X  X 

ft 

2 

ft 

a 

ft 

_J 

'v 

V. 

X 

X 

o 

a 

CL  - 

fO 

M 

a 

Q- 

ft 

ft 

a 

ft 

M 

n 

X  X 

*-• 

f-3 

r> 

ft 

I*  • 

*— 

r  • 

• 

»- 

X  LL 

ft 

U. 

ft- 

*_ 

« 

■* 

♦ 

t. 

ft- 

*■ 

UJ 

V* 

_i 

X  M 

ft 

_ 1 

ft 

9- 

ft 

•  - 

Z 

«T> 

ft 

*  ^ 

■ 

ft 

ft 

ft 

-* 

ft 

CO 

ft— 

► 

»-* 

a 

<s 

? 

O 

X  T 

V 

ft 

c 

a 

ft 

ft 

cc. 

O 

-J 

fu 

ft 

2 

Q.  - 

.-> 

=y 

2 

2 

ft 

ft 

CO 

ft 

UJ 

a 

ft- 

fO 

Ui 

X  X 

Li 

M 

Ui 

Ui 

9- 

ft 

ft 

u 

z 

2 

r-^ 

z 

o 

^  G- 

Q- 

ft 

ft 

ft 

Ui 

ft 

cr* 

>v 

*h  a 

7 

(. 

s 

o 

•-» 

• 

2  2 

% 

•* 

« 

2 

ft 

X 

ft 

Ui 

0 

2 

< 

t-* 

X- 

*- 

ft 

— • 

ftft-« 

CO 

— ■ 

— ■ 

ft 

IT* 

ft 

*- 

ft 

0 

> 

c 

V 

U- 

CO 

>- 

o 

x  r 

* 

<. 

h  *-* 

ft 

X 

z 

a 

2 

V. 

< 

tfl 

L3 

13 

a  - 

X 

ft 

►-» 

ft 

e> 

ft 

a- 

~ 

• 

X 

u. 

ft- 

O 

2  0. 

sJ 

2  • 

ft 

2 

♦ 

0 

CO 

> 

« 

ft 

u 

Ui 

ft 

_J 

<■> 

X 

M  « 

ft-* 

ft 

33  X 

ft 

►H 

ft 

u 

ft-  ft- 

ft 

3 

> 

<J 

0 

cr 

< 

o 

X  _J 

O 

ft 

ft 

Z  2 

X 

3; 

3 

ft- 

ML- 

ft- 

•  O 

« 

% 

ft 

(LJ 

ft 

♦H 

ft 

0 

►—  »-* 

T 

LO 

L. 

Jt 

rv 

7~ 

r 

“  . 

- 

' 

•^1 

O 

Ui 

ft- 

•K 

Ui 

ft 

X  • 

ft 

X 

ft 

CO 

u  a 

0? 

X 

*— 

k>. 

►- 

M 

O 

• 

z 

-»i r 

2 

% 

O 

ft 

ft 

Uf 

2  2 

w 

T. 

1 

CO 

ui 

X 

ft- 

*— 

K' 

►— » 

n 

ft 

Vi 

ft 

ft 

It 

* 

a 

0 

a. 

Ui 

Ui 

O 

o 

•  a 

ft 

ft 

V.  ”V 

1 

CO 

X 

u< 

_j 

• 

13 

X 

a 

n 

K  • 

« 

ft 

ft 

ft 

> 

> 

2 

> 

■*«. 

«» 

2 

►- 

ft- 

cr 

ui 

ir  m 

C>  O' 

ft 

ft 

2 

X  > 

X 

« 

1 

*1 

z 

tr 

»-4 

r> 

7- 

K'  *-• 

w 

ft 

O' 

ft 

ft 

a 

7  r 

7 

> 

r* 

UJ 

•-* 

F. 

* 

•— 

ft- 

o 

c 

_j 

X 

O' 

ft- 

u 

-  0 

ft 

ft 

_i 

r>  13 

z 

X 

a. 

2 

• 

a 

2 

ft- 

w 

-j 

»— 

UJ 

-« 

ft 

•—  » 

ft 

♦ 

cr  cr 

CO 

) 

V- 

> 

u 

v  - 

n 

<J 

O 

U- 

ft- 

Ui  9— 

X 

V 

UI 

O' 

ft 

* 

»— 

ZJ 

ft- 

X 

u 

ft- 

o 

2 

a 

a 

»-> 

— ’ 

ft 

O  ^ 

ft 

ft 

4 

N  II 

cr 

M 

It 

•« 

11 

z> 

O 

>•  a 

0 

H 

2 

cr 

ft 

ft 

ft- 

w 

w 

_J 

Oj 

Q 

l_i 

►■*  X 

u. 

2 

*-• 

ft 

u* 

ft 

ft 

U’ 

2  / 

•— 

c 

«: 

T  • 

O. 

< 

ft- 

O 

►- 

Ui  »- 

ft 

ft 

3D 

<  <1 

■< 

X 

LJ 

■*» 

X 

13 

Li 

2 

ft 

a 

ft 

« 

Ui  U 

ft- 

4* 

cr 

a 

►  - 

a 

O 

O 

v  0 

ft 

ft 

UJ 

2  T 

10 

U  1 

__i 

*J 

K1 

LJ 

ft 

K  -« 

ft 

ft 

H 

X  X 

a 

X 

Q 

LJ 

a 

<» 

ft 

ft 

r> 

lj 

tj 

CJ 

L' 

ft 

ft 

ft 

a 

O 

♦ 

• 

ft 

* 

• 

• 

rv 

ft 

ft 

2 

ftO 

3 

X 

X 

r 

V 

> 

F 

ft 

ft 

ft 

L 

—t 

C 

a 

r? 

c 

r  :> 

If' 

ft 

ft 

Li 

*1 

r< 

LJ 

Li 

LJ 

L1 

»  i 

L' 

F 

X 

If- 

F 

ft 

l 

ft 

ft 

« 

ft 

* 

• 

♦ 

ft 

r» 

O' 

e- 

U 

rj  0 

ft 

ft  ft  ft 

ft 

a 

C 

r 

o 

u 

L- 

L.' 

LJ 

L>  C» 

K' 

0  0 

c 

c 

r\ 

c> 

LJ 

LJ  L'  tJ  LJ 

U 

U  LI  U 

LJ  C_' 

t  * 

LJ 

O 

O' 

►3 

•** 

IT 

X) 

K  rc 

0 

r . 

a 

*r* 

a 

ft- 

r- 

f- 

►  ► 

►- 

r 

'-r 

<— 

»r 

rt 

o 

(Y 

fV 

fV' 

rv 

rv  cv 

rv 

( 1 

<v 

<V 

a 

r\ 

fc 

r. 

& 

e> 

«. 

CL  F 

cs. 

r. 

p- 

C- 

F 

c; 

Mi 


HD-A130  648  MULT I -TEMPORAL  ANALYSIS  OF  LANDSfiT  IMAGERY  FOR 

BATHVMETRY(U)  ENVIRONMENTAL  RESEARCH  INST  OF  MICHIGAN 
ANN  ARBOR  APPLICATIONS  DIV  F  J  TANIS  ET  AL  MAY  82 
UNCLASSIFIED  ERIM-155500-2-F  N00014-81-C-2334  F/G  8/10 


v^-'VwiN 


',SIX»  1*iN0ZlaUM  3Hi  NO  33NV«  »i»0. ' iZJ i »WdO J  d£ U  ft(C( 

U'Udfc3biS4a‘  t2)wndn*  fci  tti 

(«U312W  «•;  <2i'.  SI  NUUniOS3«  .  *Bl*.  *  »*»1  I 

*,  SiXV  "*3Iid3A  3hi  NO  3!)N  •  a  » J.VU* * *2; i»nscd  litf  liii 

l(i)iasaa)  »  ie'  (UKlIdn*  Cl;  wnoT UtS '£;3i.iM  ViiiJ 


ju«r 


2pl 


Af 

•x. 

ft 

«* 

u 

v« 

• 

S 

l» 

• 

r**  ^ 

ft 

ft 

< 

ft 

G  « 

•• 

ft 

A 

Ai 

g  r> 

ft 

A 

A 

w  » 

ft 

AJ 

•  u  ^ 

• 

ft 

z 

ft 

v«  wm 

♦ 

ft 

m 

■  Ifcl  4 

ft 

% 

* 

•>  JG 

X 

a 

• 

*M  G 

w-t 

o 

ft  •• 

• 

ft 

• 

— 

a 

•j  *■ 

tv 

:  r* 

Al 

a 

o 

^  Z  • 

ft  % 

•  ft 

ft 

H 

o 

ft  o 

Z  X 

ft  A 

«h 

• 

ftl  <9 

ft  ft 

•*  ft 

m 

• 

•% 

r  v  •• 

.J  • 

It  m 

ft 

3  _»•» 

ft  Al 

•ft 

t 

ft 

e 

Al 

wo  »l 

•  ft 

u 

t— 

b»  ^ 

w 

a  ii 

A*  •« 

♦  At 

UI 

►  ► 

ft 

•  * 

»  * 

o 

z 

A  H 

U) 

»-*  >-  4 

♦ 

>  G 

t 

© 

c  • 

A 

ft  ft  A 

ft  m 

ft 

a 

H 

•i  % 

-J 

in  U 

ft 

• 

•  ft 

•* 

A 

ft 

o 

-o  Ui 

•  a 

ft 

A 

3 

•A 

ft 

—  UJ  -i 

o 

•  •  X 

_» 

•)  ft 

• 

•  Z  ft 

*■* 

AJ 

•  ft 

O 

A  UJ 

A  ft  -»•-»  UJ 

ft  •• 

*• 

ft  ft 

ft 

A  Z  A 

O  ft 

•‘ft  ft 

Ui  ft  UI 

•  ft 

c* 

fti 

<H  J 

Z  «0  ^«u 

or  z 

a:  ai 

If 

ft 

-JO 

• 

•A  O  Z 

»H 

% 

•<  #• 

M 

K  ft  3 

•  a 

•  c 

ft  -J 

» 

a  cc 

•* 

• 

3  O 

OH  *>w  »- 

z 

•  ft 

A/ 

• 

Oft 

z 

ft  Z 

-J 

•  « 

<r 

*©  ft 

4  ou«o 

-JO  ft 

•  a 

a 

ft  c 

•  L* 

r  cz 

■c  •  - 

~ 

•* 

• 

-J* 

#■* 

«  ft 

•  ft 

C 

x  ft 

% 

<U  B 

Z  ft  © 

ft  ft  M 

ft 

K  IT  W 

w 

^  z  z 

A  ft'  ft 

ft 

«o  tn  m  h 
w  a*  ©ft 
>  >  ^«l  3 
-J  -J  ft  %  • 

cc  ir  *  • 
a  o  ^  » 

•  •  w 

I-  »-  r*- 
0  0^*1- 


r  c  at  e  #• 
*-»u/  cc 

Kft  ^ 

ft  *■  •«.  u> 


r  b* 

AJ  ft 

Z  Z  A'  ft 

A*  ft 

ft 

C* 

f  •• 

A  Z 

•  *h  r 

H  JT 

Al 

K  J 

a*  a 

h  a 

•-  w  *— Of 

H  A 

B  JC 

e  a 

ot  o 

a  a  or  o 

Or  O 

ft 

«  Z 

>  X 

ft  A 

HA  J  A 

X  A 

ft 

A  CJ  A 

ft 

O 

ft 

A 

3 

V* 

Al 

A  ft  ft 

O 

-1 

ft 

Z 

A 

• 

«t 

Ui 

« 

ft 

«  ft 

ft  A 

ft 

>  Z 

ft 

ft  ft 

ft  ft 

Al 

►*  A 

«J  IT 

d  O  Cl  ft 

UU  ft  4) 

4J  Cl 

Cl 

u. 

Z  • 

ft  ►- 

ft  > 

H  ft 

ft  •* 

A  ft  ft 

ft  X 

l‘  ft 

ft  «i 

ft  ft 

ft  O  ft 

Kt 

rt  o 

Cl  ft 

A  ft 

ft  ft  ft 

ft.  ft 

ft;  ft' 

A  U* 

ft  ft 

(5  ft  ft 

O'.  Cv 

A  A 

a 


$ 

sS 


??1 


aa 


2gx 


ftfftl  0  4 

*s 

9 

»  ft  4 

CG 

4WKy, 

cv 

L 

G 

9  9 

9/ 

ft* 

—  ft  — 

C-. 

&, 

r. 

•  j  ft 

LJ 

a 

►«  KN 

r. 

C* 

•S' 

fc-  ft 

O 

rv 

—  Ai  — 

t 

P 

T.- 

?>  rv 

C: 

<5> 

ft  SR 

ft 

ft 

t’i 

rv  c 

c. 

• 

lit 

• 

1 

1 

i  • 

«t 

O 

'9  9  9 

9 

ft 

9  9 

Ui 

a. 

•-  W  ftft' 

rv  a  9 

<t>  rv 

>- 

9 

*  *  « 

9 

9 

9 

«  « 

ft 

_J 

H  MM 

H 

Mi  Qt 

ft  -i 

Ui 

sc 

T 

LJ 

z 

X 

9 

V 

Z 

9 

o 

z 

X  X 

■4 

-J 

2 

LJ 

v« 

X 

Z'  >- 

2 

<r  ►<  »-*  _i 

7 

2 

a 

40  40 

4GOQ 

R 

X 

R» 

x  ru  rv 

«o 

ft' 

fi*  x»  ft* 

Irf* 

X 

9 

9  e 

•-• 

r> 

ft' 

p  ft  — 

4: 

ft 

— 

M 

s 

UJ 

B 

r*  ft  ft 

w 

ft 

ft 

«  f 

G 

ft 

& 

s?  rv  — * 

« 

rv 

«4 

r  6 

ft 

c 

e 

G  sir 

P 

c  e 

t 

e* 

t 

III 

i 

1 

• 

V  1 

» 

9 

«r 

9  9  9 

9 

9 

9 

9  9 

9 

a.* 

a 

0  W  ft  ftj  ft-ft  ft 

IV  « 

9 

V- 

* 

•  9  9 

* 

* 

9 

«  9 

* 

ft 

ft 

HHM 

M 

HHHJJ  ft 

_J 

ft> 

o 

1-  ft 

2 

UJ 

9 

o 

a 

ft  ft  9 

X 

ft 

-» 

X 

9  X 

UJ 

9 

ft'OM 

ft  r* 

2 

2 

ft 

X  MX 

2 

/ 

Pw 

w  «/: 

V 

rjx&o  x  rv  *r  x  rj  « 

c  (u  ium« 

ONBh^  Pr.C*G6 
GKQKKC  N  R  6B» 
OM6^(Vt-.PT%K 
GGG6>GPCwKK.e; 


ft 

«  IWILMIL  MIV4  99 

► 

ft 

B  JHH  k« 

*h  »-  rr  c?  a 

2  ft 

1- 

►•  -J  ft 

3  0  Z 

IV 

VJ 

9  -J  ft  ft. 

V  ft  M*  ►  9 

9 

J  J-J  J  _>  J 

r 

*->►■  2 

K  ►■  M  J  lU 

• 

UUUUCK 

9 

ft  9  U  H 

h  3  VDX 

O 

•0 

J  J  „•  J  If  o 

2 

«r  x  r 

2  O  ft  CO  ft 

ft 

ft 

o 

ft 

2  Z  it  2  ft  ft 

• 

a 

orcc  >  > 

a 

« 

U  U  V  u  o  o 

H  R'  4  «  ft 

ft  a  4  rv  ft 

— 

94 

*  •  4  4  •  « 

«o 

—  iQ  rv 

IA  ft  ft  ft  — 

IK 

MOO  o  oc 

m 

K  •<  f*  G  ® 

—  —  X  ft  — 

► 

u» 

CKKR  ft 

ft  ft  PA  ft  ft 

ft 

2  V  ft  2  2  2 

ft 

CD  ft  —ft 

BG  R  OG 

9 

o 

k  ec  e  k  e. 

r  r  ft  ft  o 

#■ 

o 

c  •  •  •  • 

■r 

9 

9  9  4I>  9 

•» 

a 

ft 

m 

ft 

9  IV  ft  IV  ft 

ft 

> 

ft 

i 

i 

i 

R 

R 

•  f  * 
—  9  — 


iv«  (v«  e  «v 

t  —  *  —  C,  — 

L  9 

•  ^  ft  t  f 


n  *  ►-  >  r 

<  •>  »?«*■► 

X  ^  M  m  j  V  y 
ftXftft»tOft_*X 

<  u  z>rrr* 

2  c  a  r  y 


«r  9KIV  ftft  f  M# 
•••  •  r  4  ►  e  it  x  rj 
i  y  n  h  f  r  c  •- 

r.  / 

fHV  ft  ^  tv  rv  6  9 


r,  pa  rv- 

2  ^  ■  ft  IO 

i u  ««ft‘  u  —  • 

2  •  4  *r*f“  *  <V 

H  00ftMMfV4  ^4>-« 

t  w  w  W  W  W  %•  WWW 


r*  fVl 

^  •* 
x  »in 

r*  r*  ^  4  6  *  < 

-4»8  4  -  »V- 


fti  6  4  «  ««*M44ft6fW< 
K  f»  »• 

X-  *• 


f  ^  ft  k  «.  e 
ftift'f  «  r  ft  f  r  p  ft  r  r  * 

««►  C  fwf  \  f-  r«t  e  t  f  r  ft  t 

s  c  c  r  s  *i  —  u  r?  e  c?  r  i.  f  f  r 

s  9  t  c  f  f  t  t  e  c  ft  e.  e  e.  c  e; 

lift'  0  IV  r  4  4  (k  ft  •  f  f  9  4  0  —  9 

•a  IT  ft  If  —  t*  y.  ay-yk-t.  Rftw- 
•a  IT»  ft  —  C.  ft  y  y  S  K  ^  ft 

p*  y-  41  c  *  k  iftrimfty'vaft 

®  C  —  r  c>  r  \  «  »  r  *  m:  «;  nti 


91 

r> 

• 

•  ••• 

• 

« 

•  f  • 

o 

•  *  • 

• 

• 

• 

1 

• 

• 

c 

a 

• 

a 

1 

1 

t 

m 

2 

« 

9 

9  ft  9  9 

9 

9 

9  9  9 

9 

X  X  9 

9 

X  X 

9 

X 

X-  4 

X  X 

X 

X 

0 

X 

3 

O 

•-  9 

•0 

-i 

•- 

ftp 

ft*  ft 

cr  u 

s 

ft 

ft 

X 

C  9 

ft  a 

ft  o 

uj 

w 

• 

u 

9 

c  c 

9  ft  r  u 

ft 

9 

~  a  a  a 

•  •  9  « 

a 

a 

a  c  9 

ft 

9  9A 

a 

a 

a 

a. 

a 

A. 

a 

a 

a 

9 

•» 

a* 

>  2 

UJ 

2 

U  M 

>  ft- 

>  X 

► 

9 

9 

• 

9  9  9 

ft 

9  9  9 

9 

9 

9 

r 

« 

9 

9 

9 

9 

9 

9 

• 

# 

•-»  ► 

#> 

9*  CP 

9  cr 

r  a 

9> 

a 

i 

T 

I 

T 

— 

*  - 

ft  ft  ft 

ft 

9  V  ft 

•- 

•  • 

— 

t « 

— 

— 

— 

— 

D 

Wi 

ft 

ft 

9  • 

X 

-  » 

r  a 

»»  f> 

> 

r 

r 

•• 

9  ft 

9 

ft 

ft 

9 

ft  >'  ft 

V. 

ft  UI 

“i 

«r 

a 

t 

w 

9 

at 

sc 

sc 

c: 

c: 

U 

ft  ft 

a 

ft 

9 

a1 

1 

r  ft 

>94 

ft 

UJ 

u  0 

t 

»■ 

M 

».- 

9 

« 

c 

L* 

*• 

u 

u 

ft 

ft  9 

49 

** 

— 

2 

a 

X  ft  u 

w  y  ft 

9 

X 

lil3  3 

> 

> 

•• 

(*. 

C’ 

r. 

•  4 

w 

1- 

? 

c> 

ft  i* 

c 

X 

—  ft* 

9  in 

4  ft 

ft 

9 

*  »  C  9ft 

n 

n 

c  r  c 

U 

*; 

y  *'  f. 

c 

r- 

t 

0 

;r 

> 

„  • 

.  J 

<' 

c 

c: 

O  «0 

ft 

3 

9 

2 

•t 

W  —  p^  9 

z 

z 

ft  9?  ft 

c- 

2 

a  a*  «..* 

LJ  O 

L 

4.- 

LJ 

C 

IT. 

ft 

to 

X 

•4 

— 

ft  UJ 

ft 

* 

ft 

9 

r«ffia3Kse< 


IC3VC  t*? 


2p 


c 

a 


IV>  HMV  S  9  tf  fi  IV  &  fU  «  *V  IV  «•  t  «  t,  O  o  «»  W  K  8>t  GO  «6KPf.OOi<rftj 

C  <  PsG  p  t  «.  9  9  ^  O  «MT  9  IT  *<|f  IV  IS  ft.  ^  r\B,  9  «C 

G  St  ft*  61  <5.  'VMV-  G  •*  f.  •■  P  ^  G  f  «•  B  ^  ^  ^  ^  ®  <;  c  t  ^  ^  f>  r.  -  N 

r  C  f  S  c>  f  f,<  m  ^  ^  e  «.  p  r  fs  p  «<  f  c-  ft  f .  p  g  »^  f,  f>  f»  r.  c;  ts  r  g  K  C  r  ^  r  tv  »  6  G  r  r  o  c% 

t  K  K  G  f.  6;  t>  •-  5;  B  ^  t'  f  -  P  c.  f-  K  6  t;  t  P  8?  C  (.  Pf  t'  t  C  t  p-  G  C  S  »i  e»  f!  t\  t  f  N  t  N  t> 

e  F.  6  <9  K  C  6  C  FKCffG&SFKfft  F  ft  F  fF  F  FF  F  C.'  R  FF  K  CF  l>  t  F  F  t  F  F  F  t  F 


#v>  k  «  <v  <r  &  a<\>9^Ffi«v<r9(i«C6  9  rvi^04)ti<><fc9A.fvio#irt\)s»  rv  tvi  a<  n  fv  ««  «  rt  f?  -e  n 

F/  •-/iF/^irfi*tvra-’(r>«‘«sa  f  k  to  p  p«  k*  «>  »»i  k  f  «  it  k  <\  tv  r  n  ^  <  g  f  o  9-  o  tv 

f‘>F-FS'«r^t\tir^Ki&9  It  ft  •*  F‘  ft»  ••  9>  K»  »-  IT  M  ••  K*  ft  X  K.  m  pa  ^  ft  •*  t. 

M«VF  O  A  IV  (f  9  ft  •-nfVWK'F**^  F  ft  t,  G  F  G  R  ft  R  ftft  F  F  n  •-  ^  «»  W  »r,  »«  9  jj  F:  41  ft  F  9  «1  n 

P  ®  D  F>  t'  P  RFfWe}fF«K>BFlVS»PiF6^«>»*prr-  K  K-  ft*  C  ft'  P  P'  C  F  ft  6'  ?»  Ai  -  ft  NRG 
GS  G  C-  E  GG6.FG®FFFG«FtFbFFF>FGFFFRFFFKFFFf  FC.&tt.Fbft  F 

€4>«<<^4l^^99««««9999/QqO  ^*9  9  9-^0  O  q  9  9  O  *t  O  99  9999  9  9  9  9 


iv  (V  fv  rv  rv  it  ft  i\i  (i  n  A'  ft'  ft'  ft  tt  tu  «x  *•  rv'  it  rv  rv  — *  rvi  ft-  it;  ft  rv  •<  ft'  ^  tv  tv  9  ft  ft  n  ft  »'  tv  ft  ^  -  n  •* 

T  y.*  -  .-.-  u  V  !*  tr  t  7  v  r*  ~  ** 

C  C  ►  J  J  *  ft-  k.  N.  ft  ti)  49  c  Cl  3t  li  3  *-h  ft-  *w  I-  U.'  *-  ftUinJU.OJ>  ft  9> 

w  ►-  ft  >  k.  u  •>  9  ►»  a  u  u)  a  a  «  _♦  u  •-  -®  3?  j  a  .j  _»  o  u>  tr  e  m  r  j  o  t  9  h  ►■  »  ►• 

X  ?  I)  r  i)  (9  tt*  A  7  C  Ota  ft  ft  J  f  ft  p*A'J7?  ft  KMf  7  ??  ttt  »  2  ft  l'.«  »■•  «»  9  m  j  h  Jkiti'  ft 

t»  c  c  .•  t  ft  o  vo  w  ft  ft  >  p  !)ft  ft  7  h  7  ?•  r  ?  ^  9  9  c  t.'  ►«  •■♦<•>  n  «*»••.,'«•  » '  i»  y  ft  r  p.  »n  **  v 

»*  <-.  »*  «-i  *-<  r>.  ►,HNMft*WH“)*ufjjjjjjJJlII.  ft'  ft.  T  2.  7  O  tt  »  fttt  !/)*/)•)•)►  D  r»  I*  K  ft 


95 


,%  .% 


?2a  js>s  e^a?'?? 


2pl 


M«  U 

B  J  IVAilV. 

M 

in  Cw  oj 

•■»  tt  tt  tt  ft 

tt 

•*n  — 

IA  «i  *■•  *  ft 

Ui 

fit  •■«  «  rv 

* 

If1  *  tt  ♦  tt>  — *fV 

a 

i;  &  «  & 

♦ 

c  »  b  <  be  t. 

c* 

OB  s 

K  ®  SCb 

o 

tt  • 

1  t  til 

tt 

ww  •< 

•-  «*  m  *n 

<S<<ICOV<G 
r><  <-n  r  <  «/<  if 
<<<  f  <  <  g  a  n  <>  o 


«  ^  ^  <c  «  r.<  q  m  <i>  a 

ui  n«**f<v«K.ir#  ♦in*«p.firNTK'^ 

lit  *IR»#R  *CLI>4C«Ci^c(-'c:(9^P 

O  »  K  f-  R  t  F  t,  S  CJ^  B 


r-  w  w  v*  r<  r>  *«  FI  %ti  in 
■  •■vKMMPvKvnwiflk.  mix  ^  iu 

V  Vt  V  <0  K  O  — *  f*  VI  O  *•  If  4'  <  I-  1C  II'  •>'  M 


«5(5I«  •««NGCCV« 

cinm  r  4<ik  xWin^nt-rNG 

5«i"  *  •*»>  "►-emxsn 

•  5?  *  «t>i.*in««««(i»WGf.iw(va 
6-  S  ■  t  •6'SR|k|i|»S^(^  K  lv  C-  f»  & 

oes  c-  :-:  c.  t'  c  &  u  t»  «.  c  e»  u  c. 


-  —  ^  •  »  »V 

•  •  K»s  -I)  M  (9  «  (b  rw(T  ** 

•  iwnwinv»*oc«»  w®  «  Af  ^  ^  »  At 

■»A*  /  .r  ft  «  i»*  •* 


IT 

n  r« 

AIM 

ww 

u> 

t  8  IV 
ms  a» 

an 

« 

ft  IV  * 

•• 

Ui 

•  Cfc  •  fV  « 

tt 

«  6S  t  tt  « 

Al 

IV  « 

D 

tttt  tt 

I  »  • 

tt  fO  — 


9<l*<D<C«lSf  ^  IV  f,  «  «  fV 
*<A>«v<0fOK  m  ••  r»  Ai  kuxi  ► 
«  «<  #  q  R  l\>  ft-  <V  MV 

»-•  tt  «-•  tt  *-««  t  M  ^  (V  IV  ^  9 
K  6f»  S  B  S  n  B  RbS&S 
6GGCiKiSBSBSfiS 


tt 

r-  tt 

Ui 

tt 

•<  & 

tt 

%  « 

*s 

tt  » 

tt 

•*•  tt* 

r* 

e  s 

-I 

ivtti  m 

ar 

C  tt 

u 

Ui  c 
•f  J 

a  ~ 

an  «• 

n  •-* 

t  — 

X  tt 

IV  tt 

n/. 

»  ft 

tt  A* 

» 

tt  tt 

tt  tt>  B  tt 

ft 

•  • 

tt 

an  ar  H'  ► 

tS 

m 

«  tt 

tt 

ft  IV  tt  tt 

A* 

3 

Ui 

li  B 

-i 

Ik 

•  ft  ft  *  K  b 

ft 

a. 

c* 

ft  ft  ft  k. 

ft 

• 

IV  « 

o 

•  •  •  » 

t 

>  ar 

•  « 

tt 

an  in  -* 

•4 

•*  ft 

H  -J 

u 

Z  • 

-♦ 

4  H 

tt  Ui 

tt 

_J 

Ik  ft 

H  ft 

-J 

u> 

ft  « 

V  « 

Ui 

« 

•  • 

v^» 

ft  i* 

tt  .i 

tt 

tt 

tt  —  a r  —  <r 

•  - 

O  tt 

►  i- 

tt 

«# 

«**••«•  a*  in  ft  ft 

tt  ui 

-J 

ft 

ft  S  ft  ft 

C 

tt 

• 

tilt 

Ui 

ft 

• 

*■-*  *•“  .  *  ~— 

C‘ 

ar. 

•n 

ay 

•4 

cr 

X 

Ui 

ft* 

tt 

%  %  ft 

at. 

ft 

tt 

ft 

tt  ~  ©  ft 

u» 

-i 

tt 

tt 

&  tt  in  ft* 

<r 

o 

tt 

tt* 

tt  90  ft  •« 

ft 

tt 

«ri 

ft 

u> 

ft 

at 

ft 

» 

ft 

ft 

4* 

a- 

z 

o 

u 

r> 

►4 

ft* 

tt 

c 

* 

ft 

UI 

Or 

ft 

tt 

V 

ft 

A.  tt  tt  ft 

ft 

r. 

li 

ft 

ft 

t\i  tt  tt 

ft- 

o 

an 

•« 

A  A  IT  IV  ft 

«.*> 

-i 

• 

ft 

fl.  fi  tt  *"• 

tt 

tt 

-» 

v 

•— 

ft  n  ft  r 

ft 

c 

z 

tt 

ft 

r 

ft  ft  ft  ft 

ft 

y 

Ui 

ft 

a 

i  »  •  • 

• 

tt 

ft 

Ui 

tt 

•* 

tt  tt  —  «M 

<p* 

o 

u 

u 

«r 

tt 

«n 

ft 

u 

e- 

*  * 

ft* 

an 

ft 

•  • 

ft 

»- 

ft* 

-i 

ft 

a 

Aft  ft  tt 

ft 

u 

z 

tt 

« 

%r 

ft  f.*  ft  ft 

tt 

J7 

i- 

1  * 

tt- 

ft  cv  tt  ft 

D 

u. 

o 

V) 

ft 

ft 

Ui 

